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Abstract 



We study the stability/instability of the subsonic travelling waves of the Nonlinear Schrodinger Equa- 
tion in dimension one. Our aim is to propose several methods for showing instability (use of the Grillakis- 
Shatah-Strauss theory, proof of existence of an unstable eigenvalue via an Evans function) or stability. 
For the later, we show how to construct in a systematic way a Liapounov functional for which the trav- 
elling wave is a local minimizer. These approaches allow to give a complete stability/instability analysis 
in the energy space including the critical case of the kink solution. We also treat the case of a cusp in 
the energy-momentum diagram. 
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1 Introduction 

This paper is a continuation of our previous work [16], where we consider the one dimensional Nonlinear 
Schrodinger Equation 



This equation appears as a relevant model in condensed matter physics: Bose-Einstein condensation and 
superfluidity (see [60], [32], [37], [1]); Nonlinear Optics (see, for instance, the survey [44]). Several nonlineari- 
ties may be encountered in physical situations: f{g) = ±g gives rise to the focusing/defocusing cubic (NLS); 
f{o) = 1 — f? t° the so called Gross-Pitaevskii equation; f(g) = —g 2 (see [46] for Bose-Einstein condensates); 
more generally a pure power; the "cubic-quintic" (NLS) (see [4]), where 



and a\, a 3 and a 5 are positive constants such that / has two positive roots; and in Nonlinear Optics, we 
may take (see [44]): 



/(*) = -«<?"' -V. /(<?) = "f GlT^-n^W)' /((?) = -aff(l+7ta D h(^ 5 ^)) I (1) 



where a, /3, 7, v, a > are given constants (the second one, for instance, takes into account saturation 
effects), etc. As a consequence, as in our work [16], we shall consider a rather general nonlinearity /, with 
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i— + 9^> + tf/(|*| 2 ) = 0. 



(NLS) 



f(g) = -ai + a 3 g - a 5 g 2 
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/ of class C 2 . In the context of Bosc- Einstein condensation or Nonlinear Optics, the natural condition at 
infinity appears to be 

|*| 2 ^rQ as \x\ -» +00, 

where r > is such that /(r 2 ) = 0. 

For solutions \t of (NLS) which do not vanish, we may use the Madclung transform 

* = Acxp(i4>) 

and rewrite (NLS) as an hydrodynamical system with an additional quantum pressure 
f d t A + 2d x <j>d x A + Ad 2 x <f> = Q { d tP + 2d x { P u) = 



d 2 A 

d t <f > + (d x ^ 2 -f(A 2 )-^ r =0 



d t u + 2ud x u-d x (f(p))-d x {^^-)=0, 



(2) 



with (p,u) = (A 2 ,d x 4>). When neglecting the quantum pressure and linearizing this Euler system around 
the particular trivial solution \& = r Q (or (A, u) = (r ,0)), we obtain the free wave equation 

d t A + r o d x U = 

d t U-2r f'(rl)d x A = 
with associated speed of sound 

c s = ^/-2rg/'(rg) > 

provided / satisfies the defocusing assumption /'(rjj) < (that is the Euler system is hyperbolic in the 
region p ~ rg), which we will assume throughout the paper. Concerning the rigorous justification of the free 
wave regime for the Gross-Pitaevkii equation (in arbitrary dimension) , see [9] . The speed of sound c s enters 
in a crucial way in the question of existence of travelling waves for (NLS) with modulus tending to r at 
infinity (see, e.g., [16]). 

The Nonlinear Schrodinger equation formally preserves the energy 

Ety)= [ \d x ^\ 2 + F(\^\ 2 ) dx, 



f r ° c 2 c 2 

where F(g) = / /. Since F(g) ~ 7r%(0 — r o) 2 ~ ~^{\/q — r o) 2 when g -> r§, it follows that the natural 

J g ° r o 2 

energy space turns out to be the space 

z = {v e L°°(R), d x tPeL 2 (R), |V|-r eL 2 (M)}cC 6 (M,C), 

endowed with the distance 

dz(i>,4>) = \d x ip - d x tjj\\ L 2 (R) + I - \i>\ || L 2 (R) + V(0) - -0(0) 

The Cauchy problem has been shown to be locally well posed in the Zhidkov space {ip € L°°(R), d x ip € 
L 2 (R)} by P. Zhidkov [69] (see also the work by C. Gallo [24]). For global well-posedness results, see [25] 
and [29] . More precisely, the local well-posedness we shall use is the following. 

Theorem 1 ([69], [24]) Let £ 2. Then, there exists > and a unique solution * to (NLS) such 
that *| t=0 = * in and * - f in e C([0,T*),^ 1 (R)). Moreover, £(*(*)) does not depend on t. 
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The other quantity formally conserved by the Schrodingcr flow, due to the invariance by translation, is 
the momentum. The momentum is not easy to define in dimension one for maps that vanish somewhere 
(see [10], [11]). However, if ip does not vanish, we may lift tp = Ae 1 ^, and then the correct definition of the 
momentum is given by [45]: 

P(V) = / WW) (l-r%)dx= f (A 2 - rl)d x <i> dx, 
JR v W\ 7 Jr 

where (-|-) denotes the real scalar product in C. We define 

Z hy = {v e -Z,inf \v\ > 0}, 

which is the open subset of Z in which we have lifting and where the hydrodynamical formulation (2) of 
(NLS) is possible through the Madelung transform. It turns out that if the initial datum belongs to Zh y , 
the solution of (NLS) provided by Theorem 1 remains in Z^y for small times, and that the momentum is 
indeed conserved on this time interval (see [24]). 

1.1 The travelling waves and energy-momentum diagrams 

The travelling waves with speed of propagation c are special solutions of (NLS) of the form 

*(t,.x) = U(x - ct). 

The profile U has then to solve the ODE 

d 2 x U + Uf(\U\ 2 )=icd x U (TW C ) 

together with the condition |?7(x)| — > r as x — > ±oo. These particular solutions play an important role 
in the long time dynamics of (NLS) with nonzero condition at infinity. Possibly conjugating (TW C ), we 
see that we may assume that c > without loss of generality. Moreover, we shall restrict ourselves to 
travelling waves which belong to the energy space Z (so that \U\ — > r at ±oo by Sobolcv embedding 
F^K) =-> C (M,C) = {h e C(R,C), lim ±00 h = 0}). For travelling waves U c e Z that do not vanish in E, 
hence that may be lifted U c = Ace 1 ^", the ODE (TW C ) can be transformed (see, e.g., [16]) into the system 

0*0c=£-— 2d 2 x r lc + V' c (7 lc )=0, with Vc = A 2 c -r 2 , 

and where the function V c is related to / by the formula 

V c (0^c 2 e-4(r 2 +0F(r 2 +0- 

To a nontrivial travelling wave U c is associated (see [16]) some £ c > — r 2 , such that V c (^ c ) = 7^ V^(^ c ) and 
V c is negative between £ c and — rjj, and rj c varies between and that is there holds {infjR \ U C \, sup K \U C \} = 
{ r a, \Ao + ?c}- Moreover, the only travelling wave solution (if it exists) that vanishes somewhere is for c = 
and is called the kink: it is an odd solution (up to a space translation) and then £0 = 0. 

We have also seen in [16] that any travelling wave in Z with speed c > c s is constant; and also that any 
nonconstant travelling wave in Z of speed c* £ (0, c s ) belongs to a unique (up to the natural invariances: 
phase factor and translation) local branch c^U c defined for c close to c* . 

In [16], we have investigated the qualitative behaviours of the travelling waves for (NLS) with nonzero 
condition at infinity for a general nonlinearity /. A particular attention has been payed in [16] to the tran- 
sonic limit, where we have an asymptotic behaviour governed by the Korteweg-de Vries or the generalized 
Korteweg-de Vries equation. In order to illustrate the very different situations we may encounter when we 
allow a general nonlinearity /, we give now some energy-momentum diagrams we have obtained (one is taken 
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from the appendix in [17], where we have performed numerical simulations in dimension two for the model 
cases we have studied in [16]). 

• The Gross-Pitaevskii nonlinearity: f(g) = l — g (see figure 1). 




Figure 1: (a) Energy (dashed curve) and momentum (full curve) vs. speed; (b) (E,P) diagram 
• A cubic-quintic-septic nonlinearity (I): f(g) = — (g — 1) + Mq — l) 2 — §(f? — I) 3 ( see figure 2). 




- p 



Figure 2: (a) Energy (dashed curve) and momentum (full curve) vs. speed; (b) (E, P) diagram 

• A cubic-quintic-septic nonlinearity (II): f(g) = — 4(g — 1) — 36(g — l) 3 or f(g) = —4(g — 1) — 60(g — l) 3 . 
For these two nonlinearities, the graph of E and P vs. speed c is given in figure 3, but the (E, P) diagrams 
are respectively those in figure 4. 
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Figure 3: Energy (dashed curve) and momentum (full curve) vs. speed 




Figure 4: The two (E, P) diagrams 
• A cubic-quintic-septic nonlinearity (III): f(g) = — \{g — 1) + \ {q — l) 2 — 2{g — l) 3 (see figure 5). 



■5 



E 








^P 







Figure 5: (a) Energy (dashed curve) and momentum (full curve) vs. speed; (b) (E, P) diagram 
• A degenerate case: f(g) = -2(g - 1) + 3(g - l) 2 - i(g - l) 3 + 5(g - l) 4 - 6(g - l) 5 (see figure 6). 





c->0 



cusp 



Figure 6: (a) Energy (dashed curve) and momentum (full curve) vs. speed; (b) (E, P) diagram 

• A perturbation of the previous degenerate case: f(g) = —2(g — 1 ) + (3 — I0~ 3 )(g — l) 2 — 4(g — 1 ) 3 + 5(g - 
f) 4 - 6{g- l) 5 (see figure 7). 



G 



Figure 7: (a) Energy (*) and momentum (+) vs. speed; (b) qualitative (E,P) diagram 

A saturated (NLS): f(g) = exp ( -\ — 1 with g a = 0.4 (see figure 8). 

V On > 




c->0 



Figure 8: (a) Energy (dashed curve) and momentum (full curve) vs. speed; (b) (E, P) diagram 



Another saturated (NLS): f(p) = —(-. 

v 2 V 1 + o 



1 



(1 + Q/Po) 2 (1 + 1/qo) 



, ; , with g Q = 0.08 (see figure 9). 
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Figure 9: (a) Energy (dashed curve) and momentum (full curve), (b) (E,P) diagram 
• The cubic-quintic nonlinearity: f(g) = — (q — 1) — 3(g — l) 2 (see figure 10). 




p 



Figure 10: (a) Energy (dashed curve) and momentum (full curve), (b) (E,P) diagram 

Through the study (in [16]) of these model cases, we have shown that if the energy- momentum diagram is 
well-known for the Gross-Pitaevskii equation, the qualitative properties of the travelling waves solutions can 
not be easily deduced from the global shape of the nonlinearity /. In particular, even if we restrict ourselves 
to smooth and decreasing nonlinearities (as is the Gross-Pitaevkii one), we see that we may have a great 
variety of behaviours: multiplicity of solutions; branches with diverging energy and momentum; nonexistence 
of travelling wave for some cq € (0, C s ); branches of solutions that cross; existence of sonic travelling wave; 
transonic limit governed by the (mKdV) or more generally by the (gKdV) solitary wave equation instead of 
the usual (KdV) one; existence of cusps... 
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We investigate now the behaviour at infinity of the nontrivial travelling waves, which depend whether 
c = c s or not. We consider for m G N the following assumption: 

fU)(r 2 ) „• c 2 

(A m ) f is of class C m+3 near r\. Moreover, for 1 < j < to + 2 , we have J +1 °; r o = (- 1 ) ^ but 

— rr~?"n + 7^ (_i) m + 3 ^ (note that for 7 = 1, equality always holds by definition of the speed of 

(to + 3)! 4 

sound c s = V-^o/Vo))- 



Proposition 1 Let U c G Z be a non constant travelling wave of speed < c < c s . 

(i) If c — 0, then there exist ^ £ I such that e l ^°Uo is a real-valued function and there exist two real 
constants M 7^ (depending only on f and £0) and Xo such that, as x — > ±oo, 

e l4 "'Uo(x)Tro ~ M exp ^-c s |x-x |) i/ £0 = -To, e^C/oOc)-^ ~ M exp (-c s |ar-a;o|) if £0 7^ -To- 

(m) If < c < c s , then U c does not vanish, hence can be lifted: U c = A c e l ^ c . Furthermore, there exist four 
real constants M c , 9 C (depending only on f, c and £ c ), xq and O such that, as x — > ±00, 



2rf 
c 

and 



\U c (x)\ 2 - rg = n c (x) ~ ^a^(z) ~ M c exp ( - ^cg-c 2 ^ - z |) , 



</>(x) - 0o T B c sgn(a;) ^ exp ( - ^/c 2 - c 2 \x -x \). 

2r oV c s- c V 7 

(Hi) If c = c s then U Cs does not vanish, hence can be lifted: U Cs — A Cs e^ c = . We assume that there exists 
to G N smc/i i/iat (.4 TO ) is verified and denote 



. 4 

Am - 2(m+l) 
r 



2(m+2) 2 

(m + 3)!* (r ° j + ( lj 4J 



TTien, we have, as x — > ±00, 

2r, 



l^(,)| 2 - r 2 ^ *.(*) - - ^)( (m + 1)2|Aro| , 2 ) ^ 

and 

sgn(x) In \x\ if to = 1 



x c,sgn(C c J / 4 \^+ 
nX> ~ 2r 2 V(TO + l) 2 |A m |J ] 



jSgn(x)|a;| "•+! if to > 2, 



TO 

and ifm = 0, there exists 6 Cs G M and ^ gK smc/i f/iai 

2c 

0(a;) - 0o T B Cs - sgn(&.) ' . 

'0 1 1 3 - 

in particular, since we impose U Cs G Z, we must have to G {0, 1, 2}. 

For the Gross-Pitaevkii nonlinearity (f(g) = 1 — g), we may compute explicitly the travelling waves for 
0<c<c 5 = \/2 (see [67], [10]) 



2 - c 2 , / y/2- 



^c(x) = \ — - — tanh (x- 



2 V 2 / y/2' 



1- 
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up to the invariances of the problem: translations and multiplications by a phase factor. On this explicit 
formula, the decay of the phase and modulus can be checked. In particular, as x — > ±00, we have 



Remark 1 In the above statements, the constants 4>o and xo reflect the gauge and translation invariance. 
In the spirit of the model cases proposed in [16], for 

f{e) = -2(<? - 1) + 3( Q - 1) 2 - A( e - if + 5( e - 1) 4 - i2( e - 1) 5 , 

we obtain a smooth decreasing nonlinearity tending to —00 at +00 (thus qualitatively similar to the Gross- 
Pitaevskii nonlinearity) for which we have 7*0 = 1, c s = 2, and V Cs (£) = — 4£ 4 — 8£ 5 . For this nonlinearity /, 
there exists a nontrivial sonic travelling wave of infinite energy (corresponding to £ Cs = —1/2), since m = 3. 

The aim of this paper is to investigate the stability of the travelling waves for the one dimensional (NLS). 
We recall the definition of orbital stability in a metric space (X,dx) for which we have a local in time 
existence result. 

Definition 1 Let < c < c s and U c G Z be a nontrivial travelling wave of speed c. We say that U c is 
orbitally stable in [X,dx), where X C Z, if for any e > 0, there exists 5 > such that for any initial datum 
\I> ln € X such that dx(^ ln , U c ) < 5, any solution ^ to (NLS) with initial datum \I> m is global in X and 

sup inf d x (*(t),e ie U c (--y))<e. 

t>o yes. 

9 6 R 

In the sequel, U c will always stands for a nontrivial travelling wave, and we freeze the translation invari- 
ance by imposing that \U C \ is even. Moreover, the solutions of (NLS) we consider will always be those given 
by Theorem 1. 



1.2 Stability and instability in the case < c < c s 
1.2.1 Stability for the hydrodynamical and the energy distances 

The first stability result for the travelling waves for (NLS) with nonzero condition at infinity is due to Z. 
Lin [47]. The analysis relies on the hydrodynamical form of (NLS), which is valid for solutions that never 
vanish. The advantage is to work with a fixed functional space since (n,u) — (A 2 — r%, d x <t>) € H 1 (M.) x L 2 (R) , 
whereas the travelling waves have a limit roe ±t&c (up to a phase factor) at ±00 depending on the speed c. 
The result of Z. Lin [47] establishes rigorously the stability criterion found in [12], [3]. 

Theorem 2 ([47]) Assume that < c* < c s is such that there exists a nontrivial travelling wave U Ct . 
Then, there exists some small a > such that U Ct belongs to a locally unique continuous branch of nontrivial 
travelling waves U c defined for c* — <r < c < c* + a. 
(i) Assume 



dc | c =c» 

Then, U Ct — A^e 1 ^* is orbitally stable in the sense that for any e > 0, there exists S > such that if 
* in = A in e^' n e Z verifies 

\\A in - A4 H i m + \\d x( f> in - 8 x <P4lhr) < s, 
then the solution * to (NLS) such that ^\t=o — * m never vanishes, can be lifted ^ = Ae 1 ^ , and we have 

sup inf \\\A(t) -A t (-- y)\\ H ym + \d x <l>{t) - d x <j)*{- - y)\ L ^{K)) < £• 

t>0 V^K. { J 
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(ii) Assume 

dP(U c ) 



> 0. 



dc \c=c, 

Then, U Ct — A^e 1 ^* is orbitally unstable in the sense that there exists e > such that, for any S > 0, there 
exists * in = A in e^ ln e Z verifying 

\\A in - A4 m{m + \\8j m - 9 X ^|| L 2 (K) < 5, 

but such that if >3/ denotes the solution to (NLS) with ^\t=o = then there exists t > such that ^ does 
not vanish on the time interval [0, t] but 

inf \\\A(t) - A*(- - y)\\m(R) + \d x <f>{t) - d x <f>*(- - y)| L 2 (K) } > e. 

By one dimensional Sobolev embedding H 1 (R) Co(R), it is clear that since U Ct does not vanish in R, 
by imposing || |^ ln | — \U C , \ \\h 1 (r) — \\A ln — A4\h 1 (r) small, * m does not vanish in R and thus can be lifted. 

Remark 2 We point out that C. Gallo in [24] fills two gaps in the proof of [47]: the first one concerns the 
local in time existence for the hydrodynamical system (see (15) in section 3.3) and the second one is about 
the conservation of the energy and the momentum. Furthermore, we make two additional remarks on the 
proof of [47] in section 3.3. 

Theorem 2 is stability or instability in the open set Z^ y C Z for the hydrodynamical distance 
d hy (ib, i,) = \A-A\ mw + \d x <t> - a^|| L 2 (K) + arg(t^l) , u) = Ae^, $ = Ae i +, 

which is not the energy distance. Here, arg : C* — > (— n, +n] is the principal argument. For the stability, 

it suffices to consider the phase 9 £ R such that arg ( ... — -. r ) is zero at x = 0, where y is the 

6 \e*<>U c ,(- -y)J 

translation parameter. For the instability, no matter what is the phase 9 e R. The result of [47] is based 
on the application of the Grillakis-Shatah-Strauss theory [35] (see also [13], [66]) to the hydrodynamical 
formulation of (NLS) (see section 3.3). One difficulty is to overcome the fact that the Hamiltonian operator 
d x is not onto. 

On the energy-momentum diagrams, the stability can be checked either on the graphs of E and P with 
respect to c, either on the concavity of the curve P ^ E. Indeed, we have seen in [16] that the so called 
Hamilton group relation 

dE dE(U c ) dP(U c ) 

c = or — - = c — - 

dP ' dc dc 

holds, where the derivative is computed on the local branch. Therefore, 

d 2 E _ d_dE _ dc_ 
dP 2 ~ dPdP ~ dP' 

d 2 E 

This means that we have stability when P ^ E is concave, that is < 0, and instability if P ^ E is 

d 2 E n 
convex, i.e. — — > 0. 
dP z 

Actually, the proof of [35], [47] provides an explicit control, as shown in the following lemma. 

Lemma 1 Under the assumptions of Theorem 2 and in the case (i) of stability, we have, provided dh y (^ ln , U c , ) 
is small enough, 

sup inf \\\A(t) - A»(- - y)lm<m+ \d x <j>(t) - d x <j>*(- - y)\\L?(WL)\ 

< Ky/\E(^) - E(Uc.)\ + |P(* in ) - P(U C .)\, (3) 
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as well as the control 

sup inf d hy (^(t),e ie U Ct (--y))<Kd hy ^ in ,U c J. (4) 
*>o set 
e e R 

Remark 3 The second estimate (4) is not a simple consequence of the control (3), but relies on a comparison 
to U c for some c close to c* instead on a comparison to U c , (this idea has also been used in [68]). It follows 
that, in the definition of stability for U Cr , one can take 5 = 0{e). 

Let us lay the emphasis on the fact that Theorem 2 of Z. Lin is given in the hydrodynamical distance 
dhy, which is not the energy distance dz- As a matter of fact, the Madelung transform 

^ :(Z hy ,d z )3U^ (^jSj) eff 1 (R,R)xL 2 (R,R)xS 1 , 

where U = Ae 1 ^ , r\ = A 2 — r 2 , and u = d x (j> is not so well behaved. 

Lemma 2 (i) The mapping ^ : (Zh y ,dz) — > i? 1 (K, K) x _L 2 (M,M) x S 1 is an homeomorphism. 

(ii) There exists <f> t G C 2 (M,]R) swc/i f/iaf d x (/>* G L 2 (R) and a sequence (4> n )n>i of functions in H 1 (R,R) 

such that, when n — > +00 , 

< d hy (e^ , e**» ) but ! ^ ! -> +°° • 

Therefore, Ji~ x is not locally Lipschitz continuous in general. However, for the stability issues, we 
compare the dz and the dhy distances to some fixed travelling wave U„, which enjoys some nice decay 
properties at infinity. Let us now stress the link between the two distances dh y and dz in this case. 

Lemma 3 Let < < c s and assume that U* G Z is a non constant travelling wave with speed c* that 
does not vanish. If c* = c s , we further assume that assumption (Ao) is verified. Then, there exists some 
constants K and S > 0, depending only on £/*, such that for any ip G Z verifying dz(ip, U*) < 6, we have 

-^d hy (ip,U*) < d z (tp,U*) < Kd hy (ip,U*). 

An immediate corollary of Lemma 3 is that Theorem 2 is also a stability/instability result in the energy 
distance. If one wishes only a stability /instability result, it is sufficient to invoke the fact that the mapping 
Ji is an homeomorphism. However, the use of Lemma 3 provides a stronger explicit control similar to the 
one obtained in Remark 1 (see (3)). In particular, in the definition of stability for U c , in (Z,dz), one can 
take 6 = 0(e). 

Corollary 1 Assume that < c* < c s is such that there exists a nontrivial travelling wave U Ct . Then, there 
exists some small a > such that U Ct belongs to a locally unique continuous branch of nontrivial travelling 
waves U c defined for c* — a < c< c* + a. 

(i) U < 7 then U Cm = A»e*^* is orbitally stable in (Z,dz). Furthermore, if^(t) is the (global) 

dc \c— 

solution to (NLS) with initial datum ^ ln , then we have, for some constant K depending only on U Ct and 
provided dz(^ ln ,U Ct ) is sufficiently small, 

sup inf dz(*(t),e i9 U c „(- - y)) < K^\E(^) - E(U C J\ + \P{¥^ - P(U C ,)\, 
t>o yes. 
e r 

as well as the control 

sup inf d z (*(t),e ie U c .(--y)) < Kd z (* in ,U c ,). 
t>o yeR 
Bel 

(ii) If — ^ c ^ > 7 then U Ct — A^e 1 ^' is orbitally unstable in (Z,dz)- 

dc |c=c* 
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For the Gross-Pitaevskii nonlinearity (f(g) = 1 — g), the stability (for the energy distance dz) of the 
travelling waves with speed < c < c s was proved by F. Bethuel, P. Gravejat and J.-C. Saut in [10] through 
the variational characterization that these solutions are minimizers of the energy under the constraint of fixed 
momentum. However, in view of the energy momentum diagrams in section 1.1, this constraint minimization 
approach can not be used in the general setting we consider here. Indeed, this method provides only stability, 
but there may exist unstable travelling waves. Moreover, it follows from the proof of Theorem 2 that stable 
waves are local minimizers of the energy at fixed momentum but not necessarily global minimizers. Finally, 
we emphasize that the spectral methods allow to derive an explicit (Lipschitz) control in case of stability. 



1.2.2 Stability via a Liapounov functional 

Another way to prove the orbital stability is to find a Liapounov functional. By Liapounov functional, we 
mean a functional which is conserved by the (NLS) flow and for which the travelling wave U c is a local 
minimum (for instance, a critical point with second derivative > S Id for some 6 > 0). Such a Liapounov 

functional always exists in the Grillakis-Shatah-Strauss theory when ^P^Pc) ^ ^ ag shown by Theorem 

dc |c=c» 

A in Appendix A. Its direct application to our problem leads us to define the functional in Z hy 

M , 



x (V) = Ety) - c*P(V0 + — (p(V) - P(u c j) 



where M is some positive parameter. It turns out that ££ is such a Liapounov functional when M is 
sufficiently large. Since the proof relies on the Grillakis-Shatah-Strauss framework, we have to work in the 
hydrodynamical variables. However, by Lemma 3, we recover the case of the energy distance. 

Theorem 3 Assume that for some c* e (0, c s ) and a > small, (0, c s ) D [c* — cr, + a] 3 c i->- U c e Z is a 
continuous branch of nontrivial travelling waves m th « < 0. // 

dc |c=c» 

M> dP(U c ) >0 ' 

dc |c=c» 

there exists e > and K , depending only on U Ct , such that for any ip g Z with inf d^ Y {^,e l6 U Cf {-—y)) < e, 

v e R 



■ e R 



we have 

inf R dl y (iP,e ie U c ,(- -V)) < K(%{^)-%{U C S) 
e g r 

and analogously with rfhy replaced by dz- Consequently, U Ct — A^e 1 ^* is orbitally stable in (iJhyj^hy) o,nd in 
(Z,dz)- Furthermore, if^{t) is the (global) solution to (NLS) with initial datum \& in , then we have 

sup inf d hy (^(t) 1 e ie U c ,{--y)) < K ^ J?(^) - J?{U C ,) < Kd hy (V n ,U C ,). 

t>0 S£K 
G R 

provided dh y (^ ln ,U Ct ) is sufficiently small, and analogously with dh y replaced by dz- 

For the travelling waves for (NLS) in dimension one, this type of Liapounov functional appears for 
the first time in the paper [3] by I. Barashcnkov. However, in [3], the problem is treated directly on the 
wave function whereas the correct proof holds on the hydrodynamical variables, in particular because 
of the gauge invariance (9, M> e 101 ^. For instance, the work [3] suggests that we have stability for H 1 
perturbations, whereas it holds only for perturbations in the energy space. Finally, we fill some gaps in the 
proof of [3]. 
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1.2.3 Instability via the existence of an unstable eigenvalue 

In the Grillakis-Shatah-Strauss theory [35], the instability is not shown by proving the existence of a un- 
stable eigenvalue for the linearized (NLS) and then a nonlinear instability result (see however [36] when 
the Hamiltonian skew-adjoint operator is onto). There exists, however, some general results that prove 
the existence of unstable eigenvalues. For the instability of bound states for (NLS) (and also for nonlinear 
Klein-Gordon equation), that is solutions of the form e luJt U ul (x), the work [34] by M. Grillakis shows that the 

condition — I / |?/ w | 2 ) > is sufficient for the existence of such an unstable eigenvalue. However, 

ClUJ V J^d ' \UJ— CJ* 

the proof relies on the fact that the bound states are real valued functions (up to a phase factor) and it 
is not clear whether it extends to the case of travelling waves we are studying. Indeed, since we have to 
work in hydrodynamical variables in order to have a fixed functional space, the linearized operator does not 
have (for c ^ 0) the structure required for the application of [34]. Another general result is due to O. Lopes 
[49] but it assumes that the linearized equation can be solved using a semigroup. This is not the case for 
our problem once it is written in hydrodynamical variables (see below). Finally, the paper [48] by Z. Lin 
proposes an alternative approach for the existence of unstable eigenvalues. The method has the advantage 
of allowing pseudo-differential equations (like the Benjamin-Ono equation). However, the results are given 
for three model equations involving a scalar unknown, and it is not clear whether the proof can be extended 
to the case of systems. 

The linearization of (NLS) near the travelling wave U c , in the frame moving with speed is 

i d ^-ic,d x ^ + d 2 x ^ + ^f{\U c f) + 2{4,\U c ,)f{\U c ,\ 2 )U c , =0, (5) 

and thus, searching for exponentially growing modes ip(t,x) = e xt w(x) leads to the eigenvalue problem 

i\w - ic*d x w + d 2 x w + wf(\U c , | 2 ) + 2(w\U c ,)f'(\U Ct \ 2 )U C , = 0, (6) 

with Re(A) > and ic^O. For one dimensional problems, the linear instability is commonly shown through 
the use of Evans functions (see the classical paper [57] and also the review article [63]). For our problem, 
we look for an unstable eigenvalue for the equation written in hydrodynamical variables, namely we look for 
exponentially growing solutions (77, u) of the linear problem (written in the moving frame) 

dt-q - c„d x r\_ + 2d x ((rl + f]^)u + ryu*) = 

d t u- c*d x u + 2d x (u*u) - d x {f'{rl + r/*)r?) - dj * — =d 2 x ( t=== ) - ~f f 2^ ° t 3 ^ T } = °' 

(7) 

where (77* , ) is the reference solution. The advantage is here again to work with a fixed functional space 
in variables (t],u). Due to the term d x { — 9 2 ( = )}, this equation can not be solved using a 

semigroup, except in the trivial case where 77* is constant, hence the result of [49] does not apply. However, 
system (7) is a particular case of the Euler-Korteweg system for capillary fluids (see [5] for a survey on this 
model). We may then use a linear instability result already shown for the Euler-Korteweg system with the 
Evans function method, as in [70] by K. Zumbrun for a simplified system, and more recently in [6] by S. 
Bcnzoni-Gavage for the complete Euler-Korteweg system. 

Theorem 4 Assume that for some G (0, c s ) and a > small, (0, c s ) D [c* — a , + a] 3 c U c G Z is a 
continuous branch of nontrivial travelling waves with 

dP{Uc) >Q 
dc \ c =c, 

Then, there exists exactly one unstable eigenvalue 70 G {Re > 0} for (6) and 70 G (0, +00), that is (NLS) is 
(spectrally) linearly unstable. 



14 



Once we have shown the existence of an unstable eigenvalue for the linearized (NLS) equation (5), we can 
prove a nonlinear instability result as in [38] , [20] . Note that here, we no longer work in the hydrodynamical 
variables, where the high order derivatives involve nonlinear terms, but on the semilincar (NLS) equation. 

Corollary 2 Under the assumptions of Theorem 4, U c , is unstable in U c , + C) (endowed with the 

natural H 1 distance): there exists e such that for any S > 0, there exists ^ m e U c , + H 1 (R) such that 
||* in - U c , ||jyi(R) < 6 butif^e U c , + C([0, T*), H 1 ^)) denotes the maximal solution of (NLS), then there 
exists < t < T* such that - U c , > e - 

Since the proof is very similar to the one in [38], [20], wc omit it. We may actually prove a stronger 
instability result, since the above one is not proved by tracking the exponentially growing mode. In [21], a 
spectral mapping theorem is shown and used to show the nonlinear instability by tracking this exponentially 
growing mode, which is a natural mechanism of instability. In Appendix B, we show that this spectral 
mapping theorem holds for a wide class of Hamiltonian equation. The direct application of Corollary B.2 in 
Appendix B gives the following nonlinear instability result. 

Corollary 3 We make the assumptions of Theorem 4, so that there exists an unstable eigenmode (70,10) G 
(0, +00) x iJ 1 (R), \\w\\ H i = 1. There exists M > such that for any solution ip € C(E+,i/ 1 (E,C)) of the 
linearized equation (5), we have the growth estimate of the semigroup 

Vt > \\m\\HHR) < Me^||^(0)| H i (K) . 

Moreover, U Ct has also the following instability property: there exists K > 0, 5 > and £0 > ; such that 
for any < 5 < Sq, the solution ^(t) to (NLS) with initial datum = U Ct + Sw e U Ct + H 1 (M) exists at 
least on [0, "/q 1 ln(2e /(5)] and verifies 

|*(t) - U c , - Se^wWm^ < KS 2 e 2 ^. 

In particular, for t = 7^ ln(2eo/<5) and e = we have 

inf |||*(t)| - \U c ,\(--y)\\ LHR) > e and inf ||*(t)| - \U C ,\(- - y)\\ L ~ m > e, 



which implies 



as well as 



inf \\*(t)-e*»U c ,(--y)tm(R) >e 



inf d hy (*(t),e t0 U c ,(--y))>e and inf d z (*(t), e lf) U c , (• - y)) > e. 

y G R y G R 

8 G R 8 G R 

With the above result, we then show the nonlinear instability also in the energy space, and thus recover 
the instability result of Z. Lin but this time by tracking the unstable growing mode. 

1.2.4 Instability at a cusp 

dP 

In this section, we investigate the question of stability in the degenerate case — = 0. In [35] (see also [36]), 

a stability result for the wave of speed c* is shown when the action c S(c) = E(U C ) — cP(U c ) (on the 

d 2 S dP 

local branch) is such that, for instance, — =- = — — is positive for c, but vanishes for c = c*. In the 

dc z dc 

dP 

energy- momentum diagrams of section 1.1, the situation is different since — — changes sign at the cusps, or, 

dc 

equivalently, the action c i-> S(c) = E(U C ) — cP(U c ) (on the local branch) changes its concavity at the cusp. 
The work [19] by A. Comech and D. Pelinovsky shows that for the Nonlinear Schrodinger equation, a bound 
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state associated with a cusp in the energy-charge diagram is unstable. The proof relies on a careful analysis 
of the linearized equation, which is spectrally stable, but linearly unstable (with polynomial growth for the 
linear problem). A similar technique was used by A. Comech, S. Cuccagna and D. Pelinovsky in [18] for 
the generalized Korteweg-de Vries equation. Then, M. Ohta in [55] also proves the nonlinear instability of 
these "bound states" using a Liapounov functional as in [35]. However, in [55], it is assumed that J = T'(0) 
and is onto, which are both not true here (and there are further restrictions due to the introduction of an 
intermediate Hilbert space). Recently, M. Maeda has extended in [50] the above instability result, removing 
some assumptions in [55] . We show the instability of travelling waves associated with a cusp in the energy- 

d 2 P 

momentum diagram in the generic case where — — =/= 0. Our approach follows the lines of [50], but with 

dc 2 

some modifications since our problem does not fit exactly the general framework of this paper. In particular, 
we can not find naturally a space U Y" , and some functions appearing in the proof do not lie in the range of 
the skew-adjoint operator d x involved in the Hamiltonian formalism. We overcome this difficulty using an 
approximation argument (similar to the one used in [47]). 

Theorem 5 Assume that for some 6 (0, c s ) and a > small, (0, c s ) D [c* — a, + a] 3 c U c € Z is a 
continuous branch of nontrivial travelling waves with 

dP{U c ) = Q ^ d 2 P(U c ) 
dc \ c =c, dc 2 | c =c»' 

and assume in addition that f is of class C 2 . Then, U Ct is orbitally unstable in (Z,dz)- 



1.3 Stability in the case c = 
1.3.1 Instability for the bubbles 

When c = 0, we have two types of stationary waves: the bubbles ,when £o > — r§, are even functions (up to a 
translation) that do not vanish, and the kinks, when £o = — Tq, are odd functions (up to a translation). The 
instability of stationary bubbles has been shown by A. de Bouard [20] (and is true even in higher dimension). 
The proof in [20] relics on the proof of the existence of an unstable eigenvalue for the linearized (NLS), and 
then the proof of a nonlinear instability result. An alternative proof of the linear instability of the bubbles 
is given in [58] (Theorem 3.11 (ii)). 

Theorem 6 ([20]) Assume that there exists a bubble, that is a nontrivial stationary (c = 0) wave Uq 
which does not vanish. Then, Uq is (linearly and nonlinearly) unstable in Uq + H 1 (M) (endowed with the 
natural H 1 metric), that is there exists e such that for any S > 0, there exists <I' m G Uq + H (M) such that 
||* in - U \\m(R) < S but if * e Uq + C([0,T*),H 1 (R)) denotes the maximal solution of (NLS), then there 
exists < t < T* such that \\^(t) — foll/f^fR) > e - 

Actually, in the same way that Corollary 3 is a better instability result than Corollary 2, we have the 
following stronger instability result, which is a direct consequence of Corollary B.2 in Appendix B. 

Proposition 2 Assume that there exists a bubble, that is a nontrivial stationary (c — 0) wave Uq which does 
not vanish. Then, Uq is (nonlinearly) unstable in Uq + H 1 (R), (Z,dz) and (Zh y ,dh y ) in the same sense as 
in Corollary 3. 

Finally, we would like to emphasize that we may recover the instability result for bubbles from the proof 
of Theorem 2, relying on the hydrodynamical form of (NLS), which holds true here since bubbles do not 
vanish. Our result holds in the energy space and for the hydrodynamical distance. 

Theorem 7 Assume that there exists a bubble, that is a nontrivial stationary (c — 0) wave Uq which does 
not vanish. Then, there exists some small a > such that Uq belongs to a locally unique continuous branch 
of nontrivial travelling waves U c defined for < c < a. Then, c ^ P{U C ) has a derivative at c= 0, 



dc | c =o 
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and Uo = A^e^' is orbitally unstable for the distances dz and dh y - 

Proof. We give a proof based on the argument of Z. Lin [47], which is possible since U is a bubble hence 
does not vanish and the spectral decomposition used in [47] still holds when c = 0. Moreover, it is clear that 
the mapping c i— > (r) c ,u c ) € H 1 x L 2 is smooth up to c = 0, using the uniform exponential decay at infinity 
near c = and arguing as in [16]. Therefore, it suffices to show that dP ^ c - ) > 0. From the expression of 
the momentum given in subsection 1.2 in [16], we have, for < c < a, 

since £o > — ^o- Indeed, we are allowed to pass to the limit in the integral once it is written with the change 
of variables £ = t£ c : 

Jo r < 

since £ > — r 2 . Therefore, 



o 2 +£-V c (0 Jo r 2 +t^^-VMcy 

dp(u c ) [*> e dz 



dc \ c =o 



Jo r 2 



+tv^o~m 



> 



since ^ ^0 (Uo is not trivial). The conclusion follows then from the proof of Theorem 2. □ 

dP{U c ) 

When we know that , > 0, we may also use the Evans function as in Theorem 4 to show the 

dc | c =o 

existence of an unstable eigenmode. However, due to the fact that the kink U is real-valued, we can use the 
arguments in [20], [58]. 

1.3.2 Stability analysis for the kinks 

We now turn to the case of the kinks (£o = — r o an d Uo is odd up to a translation). Since Uo vanishes 
at the origin, the hydrodynamical form of (NLS) can not be used. The stability of the kink as attracted 
several recent works. L. Di Menza and C. Gallo [21] have investigated the linear stability through the 
Vakhitov-Kolokolov function VK, defined by 

VK(A) ee f ([-d 2 x - f(U 2 ) - X]^ 1 (d x U )j (d x U ) dx, 

where Uo is the kink, for A e (A*, 0) for some A* < 0. They show that the Vakhitov-Kolokolov function VK 

has a limit VKo when A — >• 0~. If VKo > 0, then the linearization of (NLS) around the kink has an unstable 

real positive eigenvalue. When VKo < 0, the linearization of (NLS) around the kink has a spectrum included 

in iM. (spectral stability). Note that the approach of [47] (extending [35]) does not give directly unstable 

dP dP 
eigenvalues in the case — — > 0. Recently, the link between the quantity — — and the sign of VKo has been 
dc dc 

given by D. Pelinovsky and P. Kevrekidis [58] (proof of Lemma 3.10 there, the factor \/2 coming from the 
coefficients of the (NLS) equation in [58]): 

2\/2VKo = lim^^; (8) 

c^o dc 

r 1 • • dP 

and they also prove, in a different way from 21 that we have spectral stability when lim — — < and 

c^o dc 

existence of an unstable eigenvalue (in E* ) if lim — > 0. It is shown in [58] that the limit lim c ^ dP ^ c ^ 
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does exist. Actually, they prove that the function [0, c ) 9 c i-> P(U C ) is of class C 1 and that the derivative 
at c = is also given by (see (8)) 

l im d E^HA = 2V2VK = 2V2 lim / ([-d 2 x - f(U%) - A]" 1 ^)) (d x U ) dx 
c^o dc x^o- j R V / 

= 2^2 / Im(^ W*7 dx. (9) 
J R V 9c | c =o/ 

Our next lemma gives an explicit formula of the expression (9), involving only the nonlinearity /. 

Lemma 4 Assume that Uq is a kink. Then, there exists c G (0, c s ) such that U belongs to the (locally) 
unique branch [0,c ) 3 c i->- U c G Z. Moreover, P(U C ) — > rg7r as c — » and ifte continuous extension 
[0, Co) 3c4 P{U C ) has a derivative at c = given o?y 

dP([/ c ) 8r-3 , 1 f° (e-r 2 )V 1 1 ^ 



dc |c=o 3/F(0) 2/ g 3 / 2 ^^(0) 

The advantage of the formula given in Lemma 4 compared to (9) is that it allows a direct computation of 

dP(U c ) when ^ [ s known, which docs not require to compute numerically Uq and . For instance, 

dc | c =o oc | c =o 

it is quite well adapted to the stability analysis as in [22]. Let us observe that it may happen that a kink is 
unstable (see [43], [21]). 

In the case of linear instability, [21] shows that then, nonlinear instability holds. Actually, C. Gallo and 
L. Di Mcnza prove in [21] a stronger result, where they show that the L°° norm (and not only the H 1 norm) 
does not remain small. 

Theorem 8 ([21]) Assume that there exists a kink, that is a nontrivial stationary (c = 0) wave Uq vanishing 

somewhere, and satisfying \ — — > 0. Then, Uq is (linearly and nonlinearly) unstable in the sense that 

dc | c =0 

there exists e such that for any 5 > 0, there exists \l/ ln e Uq + 7J 1 (R) such that || v E' ln — ?7o||h 1 (m) < ^ but if 
* G Uq +C([0,T*),H 1 (R)) denotes the maximal solution of (NLS), then there exists < t < T* such that 
|*(t) - U \\ L °o m > e. 

The proof in [21] relies on the tracking of the exponentially growing eigenmode. One may actually improve 
slightly the result as this was done in Corollary 3. As a matter of fact, this was the result in Theorem 8 that 
has motivated us for Corollary 3. 

We focus now on the nonlinear stability issue when there is linear (spectral) stability, that is when 

< 0. Concerning the Gross-Pitaevskii nonlinearity (f(g) = 1— £?), for which wc have - — < 0, 

dc | c =o ' dc | c =o 

we quote two papers on this question. The first one is the work of P. Gerard and Z. Zhang [30] where the sta- 
bility is shown by inverse scattering, hence in a space of functions sufficiently decaying at infinity. The analysis 
then relies on the intcgrability of the one-dimensional (GP) equation. The other work is by F. Bethuel, P. 
Gravejat, J.-C. Saut and D. Smets [11]. They prove the orbital stability of the kink of the Gross-Pitaevskii 
equation by showing that the kink is a global minimizer of the energy under the constraint that a variant 
of the momentum is fixed (recall that the definition of the momentum has to be clarified for an arbitrary 
function in the energy space), and that the corresponding minimizing sequences are compact (up to space 
translations and phase factors). In this approach, it is crucial (see [10], [11]) that -Ekink < c s Pki n k = ^ s Tq^ in 
order to prevent the dichotomy case for the minimizing sequences. However, since the energy of the kink is 
equal to 

J-rl y/~Vo(0 J-rl V r l + Z 
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whereas its momentum is always equal to r^ir, it is clear that the condition -Ekink < c s Pkink = CsTqIt does 
not hold in general, as shown in the following example. 

Example. For n > 0, consider 

f{e) = i- Q + K(i- Q )\ 

which is smooth and decrease to — oo as the Gross-Pitaevskii nonlinearity. We have r = 1, c s = \/2, 
F(g) = (1 - gf/2 + k{\ - g) A /A and 



Ell „ k . 2 £ ie . 2 £ ie > vl „ _ ^ 

for k large (the left-hand side tends to +oo), and numerical computations show that it is the case for k > 14. 
Furthermore, Lemma 4 gives 



f^ dP(U c ) 8 1 /- 1 (g-l)V F(0) \ 



Since f!p4 = — ^ K -r, it can be easily checked that the right-hand side of (10) is a decreasing 

w 2(g — iy + k(q — l) 4 

function of k tending to 

8 1 f 1 (g- l) 2 / 1 \ 

when k — > +oo (by monotone convergence). In particular, for any k > 0, we have < q that 

dc | c =o 

is the kink is always (linearly) stable. The energy-momentum diagram for this type of nonlinearity with k 
large is as in on the right of figure 4 (the left part correspond to k smaller). 

In comparison with the constraint minimization approach as in [10, 11], which allows to establish a global 
minimization result, the spectral methods as in [35, 47] allow to put forward locally minimizing properties, 
which turn out to be useful for the stability analysis in dimension 1. 

In the stability analysis of the kink, one issue is the definition of the momentum P, which was up to 
now given only for maps in -2-hy, that is for maps that never vanish, but the kink vanishes at the origin. In 
[11], the notion of momentum was extended to the whole energy space Z, hence including maps vanishing 
somewhere, as a quantity defined mod 27r, and was called "untwisted momentum" . This notion will be useful 
for our stability result. 



Lemma 5 ([11]) IfipGZ, the limit 

r+li 

R^+ca 



<P(V0 = lim [ / + (iil>\d x il>) dx - r 2 (arg(^(+i?)) - axg(^(-R))) 



exists in R/(27rrQZ). The mapping CP : Z — > R/^ttTqIj) is continuous and if xp G Z verifies infg \ip\ > (i.e. 
i> e Z hy ), then^{ip) = P(ip) mod 2^. Finally, if * e C([0,T),Z) is a solution to (NLS), t/ien <£(*(*)) 
does not depend on t. 

Proof For sake of completeness, we recall the proof of [11]. Let tp £ Z and let us verify the Cauchy criterion. 
Since \ip\ — > r a > at ±oo, we may lift %p = A ± e l<t ' ± in (— oo, — R ) and in (+Rq, +oo) for some R sufficiently 



19 



large. For R' > R > R , we thus have in M/(27rr^Z) 

' [ +R md x ^) <te-r 2 (arg(V>(+i?')) - arg(V>(-#)))~ 

■ J — TV 

j + (ii/>\d x il>) dx-r 2 ( arg(^(+i?)) - arg(^(-i?))) 
iip\d x ip) dx + {iip\d x ip} dx 

J-R' 

- r 2 ( arg(V(i?')) - arg^(iZ))) + r 2 ( arg(^(-i?')) - argty(-iJ))) 
A 2 + d x (p + dx + A 2 _d x <p- dx 

J-R' 

- r 2 (<M#) - MR)) + r 2 (<£_(-#) - 0-(-J2)) 



/ (A\ - r 2 )d x <j>+ dx + (A 2 _ - rl)d x <t>- dx. 
Jr J-r' 



The absolute value of each term is < K J ±x>±R \d x ip\ 2 + (\ip\ — r ) 2 dx thus tends to zero if R — > +00. Thus, 
ty{ip) is well-defined. The proof of the continuity follows the same lines, and allows to show that is actually 
locally Lipschitz continuous. Let ip € Z. If ip € Z and dz(ip,ip) is small enough, we have \\\ip\ — \ip\ as 
small as we want. In particular, if R a > is large enough so that \tp\ > 3r /4 for |x| > R, we have \tp\ > r /2 
for \x\ > R . As a consequence, writing tp = A±e t ^ ± and tp = A±e l( ^ ± in (—00, —Ro) and in (+Ro, +00), we 
have, in M/(27rr 2 Z) and for R > R , 

»+R 



f (iip\d x iP) dx-r 2 ( arg(^(+i?)) - arg(^(-i?))) 

J — n. 

- j + R (%$\d x $) dx-r 2 ^ aig{i>(+R)) - arg(^(- J R))) 

/+Ro 
(i(<p - fyld^) + (i4>\d x (ii> - dx 
-Ro 

+ f A 2 + d x <j>+ - A 2 + d x 4>+ dx - r 2 U+{+R) - 4>+(+R) 
JRo v 

J ° A 2 _d x cj>_ -A 2 _dJ- dx + r 2 (c/ ) _(-R)-4>_(-R)) 



+ 



f+Ro 



r 

= / (iy, - i>)\d x ii>) + Wfixty - V0> dx 

J-Ro 



+ r 



■ 2 (M+Ro) - M+Ro)) + r 2 ^(-R ) - 4>-(-Ro)) 



f R r-Ro 
- I (A 2 + - r 2 )d x cP+ - (I 2 + - r 2 )dj + dx + J (A 2 _ - r 2 )d x 4>_ - (A 2 _ - r 2 )d x <j>_ dx. 



Ro 



-It 



We now estimate all the terms. For the last line, we use Cauchy-Schwarz to get | j R ^ {A\ — ro)d x 4>+ dx\ < 
K(ip)\A + — ro||L2( R )|| J 4 + 9 2: (?!) + ||i2( R ) < K(ip)dz(ip,ip), and similarly for the other terms. Moreover, using 
that (tp - tp)(x) = (tp - tp)(0) + Jq8 x (iP - tp), we get by Cauchy-Schwarz \\ip - i>tc>([-R Q ,+R ]) < |(V> ~ 
ip)(0)\ + \/Ro\\d x tp — d x ip\ L2/m < K(Ro)dz(ip,ip)- Thus, the terms of the second line can be estimated by 
K(tp, Ro)dz(tp, tp), and for those of the first line, they can also be bounded by < K(tp, Ro)dz(tp, tp)- Passing 
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to the limit as R — > +00 then gives 

W)~ W) mod 2^rg| <K{^,R Q )d z {il),i)). 

This completes the proof for the definition of To show that ^3 is constant under the (NLS) flow, we use 
that ^ G *(0) + C([0,T),_ff 1 ) and the approximation by smoother solutions (see Proposition 1 in [11]). □ 

For the stability of the kink, we can no longer use the Grillakis-Shatah-Strauss theory applied to the 
hydrodynamical formulation of (NLS), since the kink vanishes at the origin. Therefore, it is natural to 
consider the Liapounov functional ££ introduced in section 1.2.2, which becomes in the stationary case 
c = 0: 

if(V0 - Ety) + f (p(V) - P(u )) 2 . 

Since the momentum P is not well-defined in Z, we have to replace it by the untwisted momentum which 
is defined modulo 27rr 2 . Consequently, it is natural to define the functional in Z 



JTW = EM + 2Mrt sin 2 , 



2r 2 

which is well-defined and continuous in Z since sin 2 is 7r-periodic. In addition, is conserved by the (NLS) 
flow as E and 

Theorem 9 Assume that there exists a kink, that is a nontrivial stationary (c = 0) wave Uq which is odd. 
Assume also that 

dP{Uc) <Q 
dc | c =o 

Then, there exists some small /li* > such that Uo is a local minimizer of Jtf . More precisely, denoting 

={^eZ, inf <M*}, 

we have, for any ip e Vp, \ {e i8 U (- - y), e R, y e R}, 

jry>) > Jf(U Q ) = E(U ). 

The crucial point in this result is to prove that the functional controls the infimum infjg \ip\. From 

this locally minimizing property of the kink when < 0, we infer its orbital stability, provided 

dc | c =o 

we can prove some compactness on the minimizing sequences. Our method allows to infer a control on the 
distance of the solution to (NLS) to the orbit of the kink, but it is much weaker than those obtained by 
spectral methods in Lemma 1 or Corollary 1 for instance. 

Theorem 10 Assume that there exists a kink, that is an odd nontrivial stationary (c = 0) wave Uq and that 

dP{Uc) <Q 
dc | c =o 

Then, Uq is orbitally stable in (Z,dz). Moreover, if^(t) is the (global) solution to (NLS) with initial datum 
we have the control 

sup inf d z (V(t),e te U Q (--y)) < K^JX{^) - E(U ) < Kf/dz^, U ) 

t>0 V 6 R 
8 E 1 

provided that the right-hand side is sufficiently small. 
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This result settles the nonlinear stability under the condition 



dP(U c ) 



< for a general nonlinearity 



dc | c =o 

/. In particular, it may be applied to the nonlinearity / given in the example above. It shows that the 
stability of the kink holds with S = 0(e 4 ). We do not claim that the exponent 1/8 is optimal. 

For a complete study of the stability of the travelling waves, it would remain to investigate the case of 
the sonic (c = c s ) travelling waves (when they exist). The methods we have developed do not apply directly, 
and we give in section 7 some of the difficulties associated with this critical situation. 



2 Decay at infinity (proof of Proposition 1) 



For simplicity, we shall denote 



V (0 = V Cs (0 = c^- 4(^ + 0^ + 0- 



We freeze the invariance by translation by imposing \U C \ (hence also d x (j)) even, so that we can use the 
formulas in [16]. In particular, it suffices to show the asymptotics for x — > +00: the case x — > —00 follows 
by symmetry. We start with the proof of case (Hi). Under assumption (Am) and since F' = — /, we infer 
the Taylor expansion 

V(0 = c^ 2 + 4(r 2 + 0Qf(rt)e + ■■■ + T^yJ^ ft)?* 2 + ^ tf)C** + 0(^)) 



4r, 



(to + 3)! 



^(m+2)^.2^m+3 + f(m+l) ^2^m+3 + (9(£™+ 4 ) 



(TO + 2)1 



2(m+l) 
r 



2(m+2) 
'0 



(to + 3)! / (r ° ) + ( 1} 4 



m+3 



£>(£ m+4 ) = A m £ m+d + £>(£ m+4 ) 



since when (A m ) holds, all the terms 0(£ m+2 ) cancel out. The coefficient A m is not zero by assumption. 
Note that the existence of a nontrivial sonic wave, which depends on the global behaviour of V, imposes 
that A m £ m+3 < when £ is small and has the sign of £ Cs . Therefore, from the formula (following from the 
Hamiltonian equation 2d^.r] c + V'(r] c ) = 0, see [16] for example) 

x = -sgn(£ c J / 

he 

and since there holds, as r\ — > (with the sign of £ C J, 
d£ n di rv 



V(0 - A m r +3 



f 0(1) ifm = 



2 sgn(£ c 



-If 1 

1 W-A^m+l 



1 



+ < 



0(|ln 



if to = 1 



m + 1 VV-A™r +1 ^/-A m C +1 ' 

t 0(r]-^) if to 2 
(here, we use that the last integrand is 0(£~( m+1 )/ 2 ) as £ — > 0), it follows that, as a; — > +00, 



*7c„(a0 = sgn(£ c j( 



(to+1) 2 |A„ 



+ < 



if m = 
if to = 1 
if to > 2. 
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This shows the asymptotics for the modulus, or "q Ca . The asymptotic expansion for d x (p Cs is easily deduced 
from the equation on the phase 2d x (f> Cs = c s i] c J(r 2 + rj Cs ), and the phase <f) Cs is then computed by integration, 
which completes the proof of case (Hi). 

The proof of (ii) is easier. Indeed, in this case, the function V has the expansion 



V(0 = tie - 4(r 2 + OF(r 2 o +0= 0(?), 



hence 



Vc(o = v(o - (c 2 s c 2 )e = -(d c 2 )e + o(e). 

As a consequence, the result follows from the expansion, for r\ — > 



v(0 



: Sgn(£ c )- 



■ ^ 



since the integrand for the last integral is continuous at £ = 0. This yields the desired expansion for the 
modulus: 



v(0 



i- <-xp( r x c; r ' / 

v 4fFl)[fl) + V(tF^ 

= M c exp ( - a^c 2 - c 2 ) +C cxp ( - 1x\J c 2 - 



df) +C exp ( - 1x-\J c 2 - c 2 ) 



with 



= ^ c exp ( - y 



T «)^o, 

«. N/-{ 3 v.«)[^=iy?) + x/(<I-t 2 )P 1 



and hence for the phase by similar computations to those above. The proof of case (?) is similar, separating 
the case £0 = — r o 01 the kink (even solution) from the case £0 7^ — r o °f the bubble (odd solution) and is 
omitted. □ 



3 Stability results deduced from the hydrodynamical formulation 
of (NLS) 

3.1 Proof of Lemma 2 

(i) The mapping *4t is an homeomorphism. Let tp = Ae 1 ^, (tp n = A n e l ^ n ) n € Z such that ip n — > ip for 
dhy Then, A n — A — > in H 1 , d x <p n — > d x (j) in L 2 and we may assume (possibly adding some multiple of 2n 
to <j> n , that 0„(O) -> 0(0). We write, using the embedding ff 1 (R) <^->- L°°(M) for the before last line, 

dz(*Pn,Tp) = \\d x 1p n ~ d x ^\\ L 2 + \U n \ - \i>\\\ L 2 + |^„(0) - ^(0)| 

- \\e^"d x A n + iA^d^n - j+d x A - iA^d x 4>\\ L , + \\A n - A\\ L , + |A„(0)e^°) - A(0)e^\ 

< |( e **» - e^)d x A\\ L 2 + \\e^(d x A n - d x A)\\ L , + \\{A n - Ay^d x <j> n \\ L 2 + |A(e*» - ^)d x cj>\ L 2 

+ \A^{d x ^-d x ^ n )\\ L , + |K - A\\ L 2 + \(A n (0) - A(0))e*»<°>| + |A(0)(e^°> - e*< >)| 

< |( e *» - e^)d x A\\ L 2 + K\\A n - A\\ m + \(A n - A)\ m \d x <t> n \» + \\A\\ L ~\\(e^ - j*)d x <j>\\ L * 

+ \\A\ L oo \\d x <j> - d x <p n \ L 2 + |A| L oo|e^(°) - e^°>| (11) 
= ||(e^ -j*)d x A\ L . + ||A|| L oc||(e^ -e^)a x 0|| L2 + „^ +(3O (l), 
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from the convergences we have. Now observe that (f) n {x) = (j) n (0)+ J* d x (j) n (t) dt — »■ 0(0)+ J Q X d x (j)(t) dt = 0(0) 
pointwise, hence, by the Dominated Convergence Theorem, ||(e^™ — e lc ^)d x A\\ L 2 — > and similarly for the 
other term. Therefore, dz{^ n ,^P) — > as wished. 

Let now ip = Ae^, (tp n = A n e ift,n ) n e Z such that ip n -> ip for d^. Then, A n - A = \ip n \ - \ip\ ->• in 
i 2 , 9 x Vn - ^ <9:r0 m L 2 and V'n(O) — > V'(O)- Since | • | is 1-Lipschitz continuous, we infer for the modulus 

\d x A n - d x A\ L 2 = \d x \ij) n \ - ^IV-IIUs < \\d x ip n - d x ip\\ L 2. 

Moreover, tp n (0) —¥ 0(0) and this implies arg(0„(O)/0(O)) — > 0. Therefore, it suffices to show that d x (f> n — > 
in L 2 . We use the formula A 2 ^^ = (iip\d x ip), which yields 

O x q> n O x q> - ^ - - A 2 - \ l Wn\OxWn)y A ^ ft 2 ) A ^2 ^2 ' 

hence 

< l f n{ f~} A } L ~ , k - a\\ L 2 + l^fefl^ - dM„ + — ^ 1 1 v>„ - VRVI^ • (12) 

(mt R A z )(inf R A^) mf R A^ mi R ^ 

The first two terms tend to zero as n — > +00. For the last term, we use here again the Dominated Convergence 
Theorem since tp n { x ) = 0n(O) + Jq 2 d x ip n (t) dt — > 0(0) + Jo" ®xip(t) dt = -0(0) pointwise. This concludes for (i). 

Proof of (ii). Let us define 0* : R — > R by 0*(a;) = ^(lna;) 2 l x >i. Then, straightforward computations give 
d x (t>*(x) = ^^1 X >1 € L 2 (M) and, for A" > e, by monotonicity of d x 0*, 



2 -W^.^ s <!^2. (13) 



We now consider <p n : M — > R defined by (j) n {x) = 0ifx<0ora;> 3nir, <j) n (x) = x/n if < x < nn, 
(j) n {x) = 7r if nn < x < 2mr and <fi n (%) — 3tt — x/n if 2nir < x < 3mr. Then, we easily obtain 

Moreover, 

d^(e < *-,e < *- +i *») - |d*0*e^ - (9,0, + d x <j> n )e*-+»" | L2 > 10^. (e**" - 1)|| L 2 - ||d x 0„|| i2 , 
and, by our choice of 0„ and using (13), 

id X M^ n > r* *(d x <t>*f dx > ~ 

Since > W— = d hy (e*** , e'*' ) , it follows that, as wished 
y/mr V n 

d^(e i *',e i *' +< *») > -^(l + o(l)) » J—=d hy (e i +',e i +' +i +»). 
y/mr V n 

We do not know whether the mapping ^# is locally Lipschitz, but is probably not. 



24 



3.2 Proof of Lemma 3 



Note first that since £7* does not vanish, if S is sufficiently small and dz{ip, U*) < 5, then |||^>| — |{7*|||l°° < 
(l/2)inf R |[/,|, hence \ip\ > (l/2)inf R |£7*| > in R, thus ^ does not vanish, may be lifted tp = Aexp(i0), 
and we may further assume 0(0) — 0,(0) G (— 7r, +7r]. In (11), we can easily check that the terms leading to 
the "o(l)" are indeed controlled by K (ip)dh y (tp n , V0- I n other words, we have 

<W,tf.) < |(e < *-e < *')a x A.| ia + |A.| i0 o|(e i *-e i *')a x ^U2 + iir(^)d hy (^C^), 
provided dhy("07 ^*) is small enough. In order to bound the two remaining terms, we write, for x e M, 

cj>{x) - Mx) = 0(0) - 0.(0) + / d x( f>{y) - 9 x 0*(y) <*!/, 

Jo 

which implies, using that R 3 9 i->- e 10 is 1-Lipschitz and the Cauchy-Schwarz inequality, 



I _ gi(<Mx)-<Wa:)) 



< 



0(0) -0,(0) 



+ V F F - u *Il 2 - 



(14) 



Consequently, 



and 



(e^ ~e^')d x A4 L 2 < 0(0) -0,(0) \\d x A4 L 2 + \\u - u,| L 2 1| y^c^A, || L 2 



(e* - e'**) 9 X 0*[ Z ,2 < 0(0) -0,(0) ||d x 0*||z,2 + ||w - u*\lA VF|dx0*IU 2 - 



Both terms are < K(Uf)dhy(il>, £/*). Indeed, £/, e Z is a travelling wave, hence A,, d x A*, d x 0, are bounded 
functions which decay at infinity exponentially if < c < c s (cf. Proposition 1 (i) or (ii)). If c = c s , since 
assumption (.4o) is satisfied, we invoke Proposition 1 (in), which ensures that <9 X 0* and d x A* decay at the 
rate C(|x| _2 ) (d x A* decays faster actually). Therefore, \/\x\d x (j)* e L 2 . Gathering these estimates provides 



d z (ip,U*) < K(U*)d hy (iP,U*). 
On the other hand, from (12) and the estimate \\A — A^Wh 1 < dz{ip, E/*) (see the proof of (z)), we infer 

d hy {ip,U*) < K(U*)d z (ip,U,) + ^—r^\U-UJd x U4 L 2. 

lntg A* 

Using here again the estimate \tp(x) — U*(x)\ < |0(O) — 0,(0) | + ^/pr|||d x V — ^*|l 2 , we deduce 

d hy (ip,U*) < K(U*)d z (i>,U*). 
The proof is complete. □ 

3.3 Two remarks on the proof of Theorem 2 

We would like to point out two minor points concerning the proof of Theorem 2 by Z. Lin. We recall that 
the proof of Z. Lin [47] relics on the Grillakis-Shatah-Strauss theory [35] once we have written (NLS) under 
the hydrodynamical form (2), denoting ip = Ae 1 ^, (p,u) = (\ip\ 2 = A 2 ,d x (f>): 

d t p + 2d x (pu) = 



d t u + 2ud x u - d x (f(p)) - d x (^p-) 



0, 
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or more precisely, with rj = p — Tq = \ip\ 2 — r 2 , and denoting |^ the variational derivative, 

sO-'fl)- »)■ (») 

V Su I 

The first remark is that the scalar product in the Hilbert space X = H^^xL 2 can not be ((rj, u), (rj, u))h 1 xl 2 = 
J R rjfj + uu dx as used in [47], but the natural one is ((r),u), {fj,u))H 1 xL 2 = J* R VV + d x rjd x rj + uu dx. This 
requires to make some minor changes in the proof, especially not to identify (H 1 )* with H 1 . For instance, 
a linear mapping B is associated with the momentum through the formula 

Ph y (v,u) = J rju dx = ^(B(r],u), (v,u)) h i xL 2 with * 

for the (non hilbertian) scalar product ({rj, u), (rj, u))h 1 xl 2 — dx. The correct definition is actually 

f 1 ( 

Phy(?],u) = J r/u dx = -(B(ij,u), (ri,u)) x *,x with B = (, q 

where i : H 1 =— > i 2 is the canonical injection. As already mentioned in 1.2, the two points in the proof 
of [47] that have been completed by C. Gallo in [24] are that: [47] uses a local in time existence for the 
hydrodynamical system (15) in H 1 x L 2 , and not only in {p e L°° , d x p € L 2 } x L 2 ; and that the energy 
and the momentum are indeed conserved for the local solution if the initial datum does not vanish. 

The second point is that in the proof of stability (theorem 3.5 in [35]), it is made use of the fact that if 
Mel and (U n )neN € X is a sequence such that E{U n ) — > E(U) and Ph y (<U n ) — >■ Phy(M), then there exists a 
sequence (U n ) ne jq € X such that U n ~U n ^ \rt X, E(U n ) — > -E(W) and P\yy(U n ) = P^ y {U). In the context 
of bound states, the existence of such a sequence (U n )neN € follows by simple scaling in space, since then 
the momentum or charge is simply f Rd Un dx. However, for the one dimensional travelling waves for (NLS), 
the momentum P is scaling invariant. We do not know if the existence of such a sequence holds in a general 
framework, but for the problem we are studying, we can rely on the following lemma, which is an adaptation 
of Lemma 6 in [10] (see also Lemma in [11]). 

Lemma 3.1 There exists p > and K > 0, depending only on f such that for any p E (— p ,+po) an d 
fiei with \p\ < \p\, there exists w = ae iv G ^ 1 ([0, l/(2|p|)),C) verifying 

( 1 \ yl/(2|p|) 

w (°) = w {^ )> \w(0)\ = r + fi, j (a 2 - rl)d x ip dx = p 



and 



/■1/(2|P|) 

/ \d x w\ 2 + F(\w\ 2 ) dx < K\p\ 

Jo 



Proof. If p = 0, we simply take w — r . We then assume < p < po, since the case — po < p < will follow 
by complex conjugation. We then define, for some small S to be determined later, 



w(x) = ^r 2 -5 + 2p(l - |8p.x- l|) + exp i(l - \4px - 1|) + = ae ilp . 
It is clear that w € ^([O, l/(2p)], C) and that tu(0) = w(l/(2p)) = ^Jr 2 - S, thus |to(0)| = r + p provided 
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we choose 5 = —fi 2 — 2r tt = 0(|/z|). Moreover, using that the phase ip has compact support [0, l/(2p)], 

rl/(2p) /-l/(2p) 



/ {a 2 - r 2 )d xV dx = {-6 + 2p(l-\8px-l\) + jd x (l-\4px-l\) + dx 

f 1/{2v) ( \ ( \ 

= 2pl (l-\8px-l\) + d x (l-\4px-l\) + 

/•1/(2P) . , / \ 

= 2p j (l-|8pa;-l|J + a x (l-|4i)a:-l|J + 



dx 
dx. 



For the last integral, the first factor is equal to if x > l/(4p) and the second factor is equal to 4p when 
x < l/(4p). Hence, direct computation gives 



J (a 2 ~ r 2 )d x tp dx = 2p J (l - \8px - 1|J x 4p = p. 

For the energy part, notice first that 

|a 2 - rg| = | - 5 + 2p(l - |8px - 1|)J < \S\ + 2p 

is as small as we want if \S\ and po are chosen sufficiently small. Therefore, 

F(\w\ 2 )<K(a 2 -r 2 ) 2 . 

By simple computations, we have 

,l/(2p) 



,l/(4p) 



/ 

Jo 



\d x w\ 2 + F(\w\ 2 ) dx 

,l/(2p) 



< if 



y p 2 a x (l-|8pa;-l|J + + 0^1- |4px-l|J + + ( - 6 + 2p(l - \8px - 1|J J 



< /fp 3 + Jfp + jf ^ +p2 < if p 



since 6 = 0(\fi\) = 0(p), which concludes the proof. 



dx 



□ 



We then consider a sequence U n — (r] n , u n ) E X = H 1 x L 2 and show the existence of the desired sequence 
U n = {flni^n) G X- We recall that U n (resp. IX) is associated with a mapping ip n £ 2 (resp. £/*) that does 
not vanish. We have Ph y (U n ) = P(ip n ) — > P(U*), thus for n large enough, \P(ip n ) — P(U*)\ < po- For n 
fixed, we now pick R n > large enough so that 



r+oo 

/ \d x 4> n \ 2 + (|Vn| - r ) 2 < [P(V>„) - P(U»)} 2 . 



In particular, by Sobolev embedding, 



| \i>n\(Rn) ~r \< \Un\ - r | L oo ([i ^ +co)) < \j f r Wn\\ 2 + - r ) 2 dx < |P(V„) - P(C/*)|. 

We are now in position to apply (for n large) Lemma 3.1 with (p, /z) = (P({/*) — P(ip n ), \"4> n \{Rn) — ?*o)- This 
provides the mapping u; n e -H^QO, l/(2|p|)), C). Since \ip n \{Rn) — fo — > 0, for n large enough, there exists 
6 n e K such that ijj n (R n ) — e l(>n \ip n \(R n ) = e* e "(r +/i) = e j9 "w„(0). We then consider the mapping i/>„ G Z 
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defined by 



ijj n (x) = < 



tpn{x) 

e l0n w n (x - R n ) 



1 



if x < R n 

if R n < x < R n + 

if x > R n + 



2|P«,„)-P(E/„) 



2\P(ip n )-P(U« 



2\py> n )-p(u*)\. 

From the construction of w n and the phase factor 6 n , ip n is well-defined and continuous. It is clear that 

,l/(2p) 



P(V„) = P(i>n) + 



(a 2 n - r 2 )d x tp n dx = P(ip n ) +p = P(U.) 



for every (large) n, and that 
-V(2p) 
/o 



/■V(2p) 

E(i> n )=E(ip n )+ / |a x «;| 2 + F(| W | 2 ) dx = E(U,) + o(l) + 0(\p\)=E(U.) + o(l) + 0(\P(U.)-P(il> n )\) 
Jo 



converges to P{U*) as n — > +00. Denoting M„el the hydrodynamical expression of ^>„, it remains to show 
that U n — li n — > in X = i? 1 xi 2 . We thus compute, with the definition of V>„, 



l^n Unix 



+00 



d x \ip n \ ~ d x \4> n \ + (\ip n \ - |-0„|) 2 + (u n - u n ) 2 dx 



<2 



+00 



r+oo 

L 



<4K \d x ^ n \ 2 + (\^ n \-r ) 2 dx + 2 



l/(2\P(U.)-PWn)\) 



c w n \ +(\w n \ -r Q ) dx 



<4K[P(ip n ) - P(U*)} 2 + K\P{^ n ) - P(U*)\ 0. 



For the before last inequality, we have used that for \x\ > ip n has modulus uniformly close to ro, hence 
|d x |V>ra|| 2 +u 2 l < K\d x ip n \ 2 . Note that the construction still holds for the energy distance, the computations 
being similar. 

3.4 Proof of Lemma 1 

Proof of estimate (3). Instead of concluding the stability proof as in [35], we can notice that we have 
actually the bound 



E hy (U) E hy (U c J > 1 inf \U - W c .(- - y)\\ 

IS. yGM. 



(16) 



as soon as P(U C J = P hy (Z4J = Ph y {U) and W G <^ £ = |v G X, inf^R ||V - W c » (■ - y)\x < e| for some 

small e. If ^ m does not have momentum equal to Phy(W c „), we use Lemma 3.1 to infer that there exists ^>(t), 
with momentum equal to P h y{U c „) = P{U C ,), and such that £(*(f)) - = 0(\P(V(t)) - P(U C J\) 

and dby(V(t),V(t)) < 0(y/\P(^{tj) - P(U C J\). Therefore, for t > 0, denoting * hy (t) G X and * hy (t) G X 
the hydrodynamical variables for '5 and 'J'(i), 

^ «*hy(i) - U c , (• - y)\\ < inf [||* hy (i) - U c , (• - y)| + ||* hy (t) - * hy (i)|| 



< VKy/Em))-E(u c .) + o(y/\pm)) - p(Uc.)\) 

y/\E(*{t)) - E{U C ,)\ + |P(*") - P(U C ,)\ + V\P(* in ) ~ P(Uc,] 



<K 
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which yields (3). 

The above estimate is optimal when P(* m ) = P(U Ct ) since P c „ is a critical point of the action E — c*P. 
This bound shows that, in the definition of stability, one has to take S = 0(e 2 ) in general. The estimate (3) 
shows that one can actually take 5 = 0{e). 

Proof of estimate (4). The point is to compare *S?(t) to U c with c ~ such that P{U C ) — P(\& in ) instead 

dP 

of comparing to U c . In other words, we replace *f?(t) by U c . Note first that since — — < 0, there exists, 

dc | c =c» 

by the implicit function theorem, such a c ~ c*. We then proceed as follows. Let \I> ln e Z be close to f/ c „. 

dP 

Then, there exists c = c(* m ) ~ c* such that P(U C ) = P(* m ). Moreover, since — ^ 0, it follows 

dc \c=c, 

Wc-UcJ < K\c-c,\ < K\P(U C )-P(U C J\ = K\P(*' m )-P(Uc,)\ < K\\^ y -U c J < Kd hy (¥ n ,U c J. (17) 
From (16), it comes 

E hy (U) - E hy (U c ) > 1 inf \U-U C {- - y)\ 2 

K J/6R 

as soon as Ph y (W) = Ph y (W c ). The fact that the constant K can be taken uniform with respect to c for c 
close to comes directly from the proof in [35]. Therefore, for t > 0, 

inf I'M*) -Hc.(- - y)| < mf [||* hy (t) - Wc(- - y)|| + !«c(- - !/) - Wc.(- - 

< VKy/E(V(t)) - E(U C ) + 0(|P(* in ) - P(U C ,)\). 

Using that P(\I>(i)) = Ph y (vE'hy(i)) = Ph y (W c ) and that U c is a critical point of the action E\ iy — cPh y , we 
infer £(*(*)) - E(H C ) - [P hy - cP hy ](1^) - [P hy - cP hy ](W c ) = Ofltt* - U c \\ 2 ). Consequently, 

inf |*hy(t)-«c.(--tf)| < ^(lK y -Well + |K y -Uc,\\) < Kd hy (^ in ,U c ,) + K\\U c -U c J < Kd hy (* in , U c J, 

by (17). This gives (4). 

4 Instability result for cusps: proof of Theorem 5 

In this section, we set T c = E^ y — cPh y and we assume 

d 2 P c (U c ) = dP(U c ) =Q j3 _ rf 2 P(t/ c ) ^c(P c ) 

tic 2 | c =c» cfc |c=c» dc 2 |c=c» dc 3 |c=c,' 

The approach is reminiscent to the proof of M. Maeda [50]. Several modifications are necessary since for 
the skew-adjoint operator J = d x , we can not find the required Hilbert space Y. More degenerate cases can 
probably be considered as in [50]. 

We shall denote I : X -> X* and I H i : H 1 -> (P 1 )* the Riesz isomorphisms and U = G X = 

iP^R.R) x L 2 (M,R) and P = L 2 (M,M) x L 2 (M,M), endowed with its canonical scalar product. They are 
the corresponding Hilbert spaces needed in [50]. We consider the symmetric matrix 

1 

1 

which is such that B 2 = I 2 and 2P hy (W) = (BU,U) H - 
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Our assumption ^P^JJc) _ q / will simplify a little the computations in [501 . The 

dc | c =c» dc 2 |c=c. 

functions »i and 772 used in [50] become now 

"1(7) - ^ +7 (^+ 7 ) - J- C .(W C .) - 1 d ^A^ ~ _^P. 

and 

^2(7) = ^ = -p(c/c, +7 ) + m.) - - yA. 

In order to clarify the dualities used in [50], we provide some elements of the proof adapted to our context. 



Lemma 4.1 There exists 70 > small and a : (— 70, +70) — > K with (7(7) <~ — 0||7/ ,* 2 and swc/i t/iai, /or 



7 2 P, 

any 7 e (-70, +7o), 

,+7) — Phy(M*)- 



Proof. We have 

Phy(^c»+7+0"]BWc t +7) = ^(V>U Ct+ j+aU Ct+1 ,U Cii+1 +0-BU Ct+1 )H = P\-iy(Mc,+~ l )+0-Wc+~ l \ 2 H+ ~' i ^hy (^c„ +7) • 

Since |M*0j|f ^ 0, the conclusion follows from an easy implicit function argument near a = 7 = 0. In [50], 
the linear mapping B is seen from X to X* , but here, there is no confusion to define U Cr+7 + aVU Cr+7 € 
H = L 2 x L 2 . □ 

We define, for 7 e (-70, +7o), 

W( 7 ) =Z4, +7 + o-( 7 )BW c , +7 , 
which then verifies by construction Ph y (W(7)) = Ph y (W*). 

Lemma 4.2 As 7 -> 0, tfiere fto/ds J" C »(W(7)) - T C ,{U C ,) -^-P*. 

6 

Proof. Using that F'{U Ct+1 ) = 0, Ph y (W( 7 )) = Ph y (U c ,) = -^SIA and a{j) = 0( 7 2 ), we have 

dc |c=c* 

by Taylor expansion 

•Fc. (W(7)) - ^c. ) - Tc, +1 {U c , +1 + <r( 7 )]B« c . +7 ) - To, (U c , ) + 7 P hy (W( 7 )) 

= P c , +7 (Z4, +7 ) - J- c . (Wo.) - + 0( 7 4 ) ~-^P*, 

ac |c— 

as wished. □ 
We recall that we have defined the tubular neighbourhood ^ £ = {V £ X, inf yS R ||V — U*{- — y)\x < £}■ 

Lemma 4.3 For e > small enough, there exist four C 1 mappings 7, a 7 y : ^ e — > K and $ : ^ — > X, 
satisfying, forli € 

«(■ - y(W)) = W(7(W)) + <>(«) + a(W)BW c , +7(w) 
and i/ie orthogonality relations 



{■&{U),d x U c , + ^ u) ) H - [9 c W c ]| c=Ct+7(w) ) ff = (^(W),BW Ct+7(w) ) H = 0. 

Or; 



Pma%, I^f G P 2 x P 1 and 1^ §2 e #4. 
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Proof. We consider the mapping G : X xRx (-70, +70) xM^K 3 defined by 

/ («(■ - y) - W( 7 ) - aWA c , +1 , d x Uc, +7 )H 
G(U,y,j,a) = [U{- - y) - W( 7 ) - oMA c , +1 , [d c U c ] w=Ct+7 ) H 
\ {U(--y)-W{ 1 )-cMU Cr+1 ,MU c , +1 ) H 

Then G(l{*, 0, 0, 0) =0 since W(0) = W*. In order to show that G is of class C 1 , we have to pay attention 
to the translation term U{- — y), since differentiation in y requires U e H 1 x H 1 whereas we only have 
U E X = H 1 x L 2 . It thus suffices to write 

/ {U,d x U Ct+1 {-+y)) H - (W^) + aMU c , +7 ,d x U c , +J ) H 
G(U,y,-r,a) = (W, [0 c M c ]| c=c .+ 7 (- + y)) H - (W( 7 ) + oMU c , +1 , [d c U c } lc=Ct+1 ) H 
\ {U,WA Ct+1 {-+y)) H - (W( 7 ) + aMU Ct+1 ,MU Ct+1 ) H 

to see that G is indeed of class C 1 on I x I x (~ 7o> +7o) x R since c W c is smooth. Moreover, using that 
<9 7 W| 7=0 = [<9 C Z4]| C=C „, we infer 

ar / {U*,dlU*) H -{[dcU c \ c=c „d x U*) H -{WA*,d x U*) H 

1 (w,, 0,0,0)= -(w.,ax[5 w c ]|c=c)if -|[3cW c ]|c=c.Ih -(B£4,[d c w c ], c=c j ff 



d{y,l,a) 



-{u„md x u,) H -([d c u c ] lc=c „m,) H -\\«u*\\ H 



Here, [50] uses assumption 2 (Hi), that reads for us ([d c U c }\ c = Cti ,d x U*)H — ~{d x [dMc\\c=c, , U*)h = 0. It is 
indeed the case since U c is chosen even for any c (close to c*). Furthermore, (W*, <9 2 ^4)# = — \\d x U*\\ 2 H by inte- 
gration by parts, {MA*, d x U*)H = (U* ,Md x U*) h = since Md x = J is skew-adjoint, and (BW*, [dJJc\\ c =c,)H = 
dc[P\iy{U c )\\c=c, — by hypothesis. Therefore, 

dr ( -\\d x U c ,f H 

1 (£4,0,0,0) = I -p c w c ]| C=c j 2 ff 



o{y,l,a) \ Q 








is invertible, thus the implicit function theorem provides three real-valued functions y, 7 and a, defined near 
U* (in X) and with y(U*) = j[U*) = a{U*) = 0, such that G(U, y{U), j(U), a(U)J = 0. These functions 
are extended to (3 e (for e small enough) by the formulas y(U) = y(U(- — y)) + y, *f(U) = "f{U(- — y)) and 
a(U) = a(U(- — y)) for any y € R such that U(- — y) lies in the neighbourhood of £4 where y, 7 and a arc 
defined. Consequently, the mapping 



0(H) = «(• - y(W)) - W(7(W)) - a(U)MU Cr 



is orthogonal in if to d x U Cit+ ^^u), [dJJc\\ c =c,+^(u) an d ®£4,+7(i0> as desired. Since / is assumed of class C 2 , 
we have U c G iJ 4 and the regularities IT 1 ^' <E H 2 x H 1 and I"* ^ e -ff 4 follow easily. □ 



Remark 4.1 We would like to point out that in [50] (Lemma 3), it is claimed that u w(u)" is orthogonal to 
"cL</> w+ a(«)" ( we re f er to the notations in [50]). However, since ll T(6(u)) — *(A(u))" is already orthogonal 
to "(9 w tJ+ A( M )" by construction, this is equivalent to "(i?^ w +A(„), <9 w 9Wa(u)) = 0", or u d u >[Q(<l> u i)] = 0" 
at "u/ = to + A(u)". We have not understood why this should happen since in general, for the function 
oj' 1 y Q(4>u>'), the point oj is the only local critical point. For this reason, we have added a component to 
the original mapping G in [50]. Let us observe that then, Lemma 3 in [50] uses the assumption "d"(iv) = 
0". On the other hand, the derivative of G in [50] assumes "u € D(T'(0))", for otherwise the expression 
"Gi t \(u, 9, A) = (T'(0)T(6)u, T'(0)</> w+ a)" , for instance, is meaningless. We have therefore given some details 
showing clearly the smoothness of G. 

We now prove a lemma which shows that the quadratic functional associated with T'J gives a good control 
on i)(U) thanks to the orthogonality conditions on this function. This result is in the spirit of Lemma 7 in 
[55]. 
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Lemma 4.4 There exist < 71 < 70 and K > such that if 7 e (—71, +71) and if -d e X verifies 

($,d x U Ct+7 ) H = (1?, [9 C W C ]| C=C , +7 ) H = (tf,BZ4 t+7 ) H = 0, 
tfren {F l i+ 1 {Uc.+ 1 )W)x*,x > K \\ti\ x . 
Proof. As a first step, we prove that if i9 e X verifies i9 7^ 0, 

{ti,d x U*) H = (0, [^c]|c=c,)ff = (0,BH.)ff = 0, 

then (J 7 " ti) x-.x > 0. Indeed, assume that (J 7 " (14*)$, i9)x',x < 0. Let x G 1 be a negative eigenvector 
of J 7 ". We claim that can not have ((1?, x)h, {[9J^c]\ c =c,,X)h) = (0,0). For otherwise, {°&,x)h = implies 
that d is L 2 -orthogonal to x, which is the eigenvector associated with the only negative eigenvalue — /xo of J 7 " 
seen as an unbounded operator on L 2 , thus (J 7 " (14*)$, "&)x-.x > 0, and since we assume equality, this means 
that $ belongs to the kernel of J 7 "^*), which is spanned by li* = d x U*, but the condition (1?, d x U^)n = 
then implies 1? = 0, a contradiction. Therefore, there exists (a, b) e R 2 such that (a, 6) 7^ (0, 0) and 
(a[9 c W c ]| c=Ct + bd,x)H — 0. The nonzero vector p = a[d c U c ]^ c=Ct + bd then verifies (p, \)h = and 
(p, JU*)h = a{[dMc]\ c =c t 1 JM*)h + b($, JU*)h — 0, so that (J 7 "(U*)p,p)x»,x > 0. Here, we have used once 
again that ([d c U c ]\ c=Ct , JU*)h = since the left vector is an even function and the right vector an odd 
function. However, in view of the equality (J 7 " (U*)[d c U c ]\ c = c , , (f))x*.x = (BU*,4>)h, valid for any <fi E X 
(which follows from differentiation of E^ y (U c ) = cP hy (W c ) = c{MA c , -)h at c = c*)), we have 

(T:(U«)[d c U c } lc=c „$)x*,x = (M[dMc]\c=c,J)H = 0. 

As a consequence, 

< (T';(K)p,p) x *,x =a 2 (Jt'(^)[5 c Z4]| c=Ct , [d c Uc]\c=c,)x*,x + # {J\ x;x 

= a 2 (mU«, [d c Uc]\ c =c,)H + b 2 (T':{U«)V,V) X ',x = b 2 (T':(U«)d,d) x * ,x, 

since (BW*, [dJAc]\ c =c t ) h = d c {Ph y (U c )]\ c=Ct = in our situation. We reach a contradiction since the right- 
hand side is supposed < 0. 

We now prove the lemma by contradiction, and then assume that there exists sequences (# n ) n >i € X 
and (7„) n >i € (0,7 ) such that j n -> 0, ||i?„|^ = 1 and 

{ti n ,d x U c , +Jn ) H = {ti n , [d c Uc]\c=c,+^)H = (<?„,BW c . +7b ) h - 0, (18) 

but (F' l ! r+ln (U Ct + ln )'d n ,d n ) x* .x — > 0. Possibly passing to a subsequence, we may assume the existence of 
some 1? = (C, v) € X such that i?„ = (£„, u n ) — ^ $ in X = H 1 x i 2 . We then show the lower semicontinuity 
of (J r "(U*)'&,'&)x*,x- This is roughly a verification of part of assumption (A3) in [55], used in Lemma 7 
there. By compact Sobolev embedding, we may assume ( n — > ( in Lj^ c (R). A straightforward computation 
gives 

(5 X C) 2 d x (d x n Ct+1 ( 2 {d x n Ct - 



(-r'l +1 {u c , +1 )W) x ,, x = f — j 



+ Vc+ 7 ) (r 2 +Vc+ 7 ) 2 4(r 2 +n Ct+J ) 3 
+ 2(r 2 + ijc,+~ ( )v 2 + 2(2u c>+7 - (c. + 7 )K - /'(r 2 + ?7c,+ 7 )C 2 d&. 

Since r§ + r/ c>(+7ri remains bounded away from zero uniformly and r/ C:t+ ~ (n — > ?7 C>( in W '°°(M.) n iJ 1 (]R) as 
n — > +00, the weak convergence ( n £ in if 1 implies 



/ 



(g,C) 2 _ d x (d x n Ct C 2 (dxVc.) 2 dx 
2(r 2 + ry c J (r 2 + ry c J 2 4(^+^)3 



< li- / - ( f xU) , ff d ^ 2 + dx. (19) 

+ %»+ 7 J (^o+^+7J 4(r^+?7 c , +7 J J 
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For the remaining terms, we write, for some R > to be determined later, 

/ 2(r 2 + r/ Ct+7 Jvl + 2{2u Ct+ln - (c* + 7„))u„Cn - /Vo + J ?c,+ 7 „)C^ dx 

JR 



) 1/2 v n + 

JR L 

+ / + / K 

J\x\<R J\x\>R z v 



( 2m c,+7„ - (C* +7n))Cn 



2(r 2 +?7c , +7 jV2 
I , (2u c . +7 „ - (g, + 7 „)) 2 

^0 +^c,+ 7 „ 



c?x 



-2/'(r 2 + »7c, +7 j)c'^- 



For the first integral, we may use that (( n ,v n ) — 1 m ^ 2 x ^ 2 an d the fact that (-q Cr+ln , u c „ +7n ) 

converges to (r?* , ) uniformly to deduce 



2(r 2 +'?c, +7 J 1 / 2 



2(r 2 + ^)i/2 



in L 



(20) 



hence. 



( 2\(r 2 + n U/2,. , (2«« ~ c*)C 



dx 



< lim / 2 



(?"0 + ? 7c»+ 7 „) 1/2 Wn + 



(2Uc»+ 7 „ - (C* + 7n))Cn 



dx. 



2(r 2 + r 7c , +7n )V2 
Since C™ -> C in L°°([-ii, and (u c 

,+7„i'?c,+7 n ) ~^ ( u *j?7*) uniformly, it follows that 

1/ (2 ^ -c.) 2 

,2 



(21) 



|x|<fl 



= lim 

n— >+oo 



/ 5(" 

J|x|</i V 



2/'(^+??*))c 2 da: 
If (2u c „+ 7 „ - (c* +7«)) 2 



»0 +??c,+ 7 „ 



-2/'(r 2 + r ?c , +7n ))c 2 da;. 



For the last integral, we have to use the decay at infinity of rj c _ t+1 and u c „+ 7 uniformly for I7I small. This 
gives 

-2/'(r 2 + r ?c , +7 J 



(2w c , +7n - (c* +7«)) 2 



c 2 - c 2 



»0+*7c»+ 7n '0 

as |x| — > +00, uniformly in n. Since < c* < c s , there exist some small <5 > and some R > large such 
that, for any n and any x with |x| > i?, 



(2u Ct+7n - (c* + 7„)) 2 



?o + f 7c l+7 „ 



-2/'(r 2 +?7 c , +7 J > 8. 



In particular, since Cn ^ C m i 2 , 

(2w Ct + 7 „ - (c* +7n)) 2 



l|x|>ii 



(- 



ri +?7c»+ 7 „ 



- 2/'(r 2 + ^ +7 j) V2 C„ - 1 W >*( - (2 !* 2/'(r 2 + *)) ^C, 

/ v <* 



in L 2 , thus 



1/ (2m, -c*) 2 



|x|>fl 



< lim 

J\ x \ 



I i( 

J|a;|>_R ^ V 



2/'(r 2 + r?*))c 2 dx 
1/ (2u c „ +7ii - (c, +7«)) 2 



»o + r ?c,+ 7n 



V(r 2 o+Vc +7 j)C dx. (22) 
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Combining these three lim inequalities, we deduce 

(J*'(U.)0,0)x;x< lim {F'l +ln {U c , +ln )K, K)x;x = 0- (23) 

n— >+oo 

Turning back to our sequence (i? n ,7 n ), we may pass to the limit in (18): 

{ti,d x U*) H = (0, [^c]|c=c.)ff = (tf,BZ4)ff = 0. 

Comparing with (23), we deduce from our first claim that "3 = 0. This means that we must have equality in 
all the above lim inequalities. In particular, the weak convergence (22) is actually strong, thus £„ — > £ = 
in L 2 (M) (the strong convergence in {\x\ < R} being already known since ( n — > C in Going back to 

the equality in (19) thus provides d x £„ — > 9xC = in L 2 (R), since r 2 , + n Cr+ln remains uniformly bounded 

away from zero and by weak convergence, = f 5 ^ X ^ C "\ 3 dx = lim f ^^"^^ C »+T" dx. Finally, 

the equality in (21) means that (20) is actually a strong convergence, that is v n —¥ v = in L 2 since ( n — > C 
in L 2 . The contradiction then follows: 1 = |#„|& = ICnlla + I^Cnl^ + \\v n \\ 2 L 2 0. □ 

Remark 4.2 This result is also Lemma 7 in [50], and is said to be Lemma 7 in [55]. However, the hypothesis 
of Lemma 7 in [55] are not verified, and in particular assumption (B3) there. It is natural to believe that 
this assumption is verified is most physical situations, but it is not clear whether it always holds true in the 
general framework of [50] without further hypothesis. 

The next lemma provides a control for a(U). 

Lemma 4.5 Assume e > small enough. Then, there exists K > such that for any U € G e satisfying 
P\iy{U) — Phy(U*), there holds 



*(u)\<K(f{u)\m\x+mu)fx)- 



Proof. It is the same as in [50], Lemma 8, but give it for completeness. We expand and use that B 2 = Id 2 
and the definition W(-f) = U Cr+1 + cr('y)MU Ct+1 for the second line: 

P hy (Z4) =P hy (U) = P hy (U(- - y{U))) = fty(w(7(«)) + t?(«) + a{U)MU c ^ (U )) 
= P hy (W(7(^))) + P hy (tf(W)) + a 2 {U)P hy (MU c ^ iu) ) 

+ a{U){m{U),WA c , + ^ U) ) H + {WA c , + ^ u) ,m))H + a{U){WA c , + ^ u) ,WA c , + ^ u) ) H 
+ a(7(W))(W c , + ^),^(W))H + a(7(W))a(W)(W c , + ^),BW c , + ^))H. 

Since P hy (W{j{U))) = P hy (Z4), we infer 

-a{U) \\\U4 2 H + o(l)] = <t(j(U))(U c ^ (u) ,i3(U)) h + P hy (#(U)) 
and the conclusion follows since (7(7) = 0(7 2 ) by Lemma 4.1. □ 
Now, we give a lemma useful to estimate 

Lemma 4.6 Assume e > small enough. Then, there exists K > such that for any U <G G e satisfying 
P hy (W) = P hy (W*) and F^(U) - < 0, there holds 

wm\ 2 x<K\m\ 3 - 

In particular, \a(U)\ < K\^(U)\ 3 . 
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Proof. It is the same as in [50], Lemma 9. Note first that the last assertion is a direct consequence of the first 
one and Lemma 4.5. Next, we argue by contradiction and assume that there exists a sequence U n — > IA* in X 
such that J 7 *^)-^ 7 *^*) < and |i?(W„)||x > |7(^«)| 3 - For simplicity, we denote % = j(U n ), d n = $(U n ), 

a n = a{U n ). Then, by Lemma 4.5, we have \a n \ < K(f n \& n \ x + \\K\ X ) < K(\K\T + WMx) = O(\0 n Vx)- 
Therefore, by Taylor expansion and Lemma 4.3, it holds 

F*{U n ) - F{U,) =F,{U n {- - y n )) - F{U,) = J;(W(7«)) +ti n + a n WA c ,+^) - F{U*) 

= J 7 * ( W(7n) ) - HU* ) + {K ( W(7n) ) ,#n + a n BU c , +j n )x*,x (24) 

+ \(^{min))KJn)x^X + 0{\d{U n )\ x ). 

However, by Lemma 4.2, J'c. (W(7))-J" C .(W C J = ©(M 3 ), and since Ki^il)) = K{W{®))+o{l) = + 
o(l) = o(l), there holds (J*(Wtf n )), a n WA c ,+^ n ) x* ,x = o{\$ n f x ). Furthermore, using T[ = P c ^ n 
the third orthogonality condition in Lemma 4.3 and that (7(7) = 0(7 2 ), we deduce 

{K {W(%)),0 n )x;X = (K+7A W (*rn)),#n)x*,X+%(MW(%),$ n ) H 
= (- :F c»+7„(^c»+7„) +0"(7 n )BW c<i+ ^ n ,l? n )x%X +7n< J (7„)(W Ct+ 7 ?i ,t?„) ff 

= 0(tlMx) + 0(\%\ 3 \\K\\X) = oiWnWT) = 0(\\Mx)- 

For the last line, we have used another Taylor expansion with F' c ^ + ^ n (U Ct J r ^ n ) = 0. Finally, Lemma 4.4 
yields (F'c,+* ln (Uc,+>y n ) , &n->'&n)x*,x > ^ol^nlx- Reporting these expansions in (24) yields 

T*{U n )-HU*) > ^-\\K\x + o{\K\x) > ^Wnfx 
for n sufficiently large, which contradicts our assumption. □ 

We now need to find an extension of the functional "A" and "P" used in [50] (and also in [55]). In 
these works, these functional are built on what should be here " J~ 1 d c U c = Md~ 1 d c U c " , but unfortunately, 
d c (j) c docs not have vanishing integral over R (for instance, d c T] c has constant sign). We rely instead on a 
construction of a suitable approximation of " J _1 9 C W C ". A similar construction is used in [47]. 

Lemma 4.7 For any < K < 1, there exists a C 2 mapping T K : (—71, +71) — > X such that, for any 
7 £ (—71, +7i), T K (7) e H 2 x H 1 is an odd function verifying 

||JT K (7) - [d c U c ]\ c=Cr+7 \\ x < k. 

Proof. We fix an even function 60 € such that J R 6o dx = 1. For T > to be fixed later, but 

independent of 7 and k, we set t K = T/k 2 > and 

T K ( 7 )(aO =B / f[a«Mc]|c=c, +7 W- r ©off ) / [d c U c ] lc=c , +1 (z) dz] dy. 

Jo L r « vi k ' Jm. J 

It is clear that T K ( 7 ) e C^l) and that, since J = d x M and B 2 = Id 2 , 

JT K {l) ~ [d c U c ]\ c=Ct+1 = j-@o(— ) [ [d c U c ]\ c=Cr+1 {z) dz 

JR 

In particular, 



\\jT K (^)-d c [u c ] lc=c , +7 \\ 2 x = l|eo|| 2 + ^|9 x eo||| 2 [ I [d c u c ] lc=Ct+1 (z) dz) <« 

if we choose T = T(c*,U*, Oo) > sufficiently large and 71 smaller if necessary. Moreover, T re (7) is odd since 
U c and O are even. In addition, the even function y i->- [dj^c]\c=c,+-y(y) — t^r@o(^) j R [d c U c ]\ c=Cit+7 (z) dz 



2 



,2 
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decays exponentially at infinity (since 80 has compact support and d c U c decays exponentially) , and has zero 
integral (since 60 has integral equal to one), hence 

t K {i){x) = -B \[d c Uc]\ c=c , +1 {y) - r°o(f) I \dcU c ]\ c=c , +1 {z) dz] dy 

and decays exponentially at infinity. It follows easily from these two equalities that 7 i-> T K (7) G L 2 x L 2 is 
well defined and continuous, hence also 7 M> T re (7) G i? 2 x H 1 . By the same type of arguments, 



9T K 
c^7 



( 7 )W =B / [[a 2 W c ]| c=c , +7 (y) - 1 6 (f ) / [a c 2 W c ]| c=c , +7 (^) efel 



is well-defined and is a continuous function of 7 with values into H 2 x iJ 1 , and similarly for the second 
derivative. □ 

We now define, in the tubular neighbourhood e of W*, the functional (corresponding to "A" in [50]) 

n^w) = («(■ - y(u)),r K (^(u))) H = («, r K ^(u)){- + y(u))) H 

depending on k G (0, 1), which will be determined later. The first properties of K are given below. 

Lemma 4.8 For any < k < 1, Sl K : G e — >• M is of class C 1 . In addition, there exists some bounded mapping 
Afy : e — » X such that i/^hy € C 1 ([0,T),X) is a solution to (15) i/iai remains in e , then 

^fi*(*hy(t)) - s K (* hy (t)), 

where S K : ^ e — > R is defined by 

S«(W) = -(^.^(W), { JT K (7(W))(- + + («, a 7 T„(7(W))(- + i^M^H)})^. 

Proof The fact that Q K is of class C 1 follows directly from the second expression and the fact that y and 
7 are C 1 (in [50] formula (3.11), the same remark as for the smoothness of G after Lemma 4.3 holds, since 
it requires "u G £>(T'(0))"). If * hy = G C 1 ([0,T),X) is a solution to (15) that remains in e , we 

therefore have, denoting j(t) = 7(^hy(i)) and y(t) = y($?hy(t)), 

|fi K (*hy(*)) = (0 t *hy(*), T K (7(t))(" + 

+ (* hy (i),9 x T K (7(t))(.+y(i)))if(y / (fhy(t)),9 t * hy (i))x.,x (25) 
+ (*i v (t),a 7 T K (7(t))(. + y(t))) H <7 , (* hy (t)),ft* hy (t)) x . iX . 

We now observe that the invariance of Q K by translation provides by differentiation the equality, for W € 
= ^n„(W(- - y))| B=0 - (W, S X T K (7(W))(- + = (BW, JT K (7(W))(- + y{U))) H 



= {P{ iy {U),JT^{U)){-+y{U)))x*,x- (26) 
In particular, the second term in (25) vanishes. In addition, since ^h y = (Vi u ) G C 1 ([0,T),X) a solution to 
(15) that remains in 6 S , there holds, denoting by the variational derivative, 

($*hy(t), T K (7(t))(- + y(t)))* - ( J^(* hy W), T K (7(t))(- + y(t))) ^ 

= - (^(*hy (*)), JT K (7(t))(- + y(t))) H = -(< y (*hy (*)), JT K ( 7 (t))(- + i/(t))) x . )JC 

= -<^c.-Hr(*)(*hy(*))' JT «(7W)(- + y(t)))x>,x - (c. + 7(t))(iU*hy(t)), ^ K (7(i))(' + y(t))) x , x 

= -<^ c .^ (t) (*hy(*)), JT K (7(i))(- + v(*))>x-,x> 
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5P 

by (26). In addition, since y = B\&hy and JB = d x , 

dW hy 

= (f ( ^ hy ft)),j ^gg^ (^ y (t))) x ^ x + (c, + 7(^))(7 / (^h y (t)),a^ hy ft)) x%x . 

The second term vanishes since 7 is invariant by translation (by definition, see the proof of Lemma 4.3). As 
a consequence, 

tf(^y(t)),dt^y(t))x*,X = (j%j^(*h y (i))J-V(*hy(t))) x 

= -(%^(*^(*))'- 7rl ^(*))) JC 



-( ax 57 ^ ( * hy(t)) ,a x ji-l|(* hyW) ) i2 . 

The first term is simply -(J r ^ + - (w) (* hy (t)),JI- 1 7'(* hy (t))) x%x . We then define A/" 7 : G e -> X by J\f^(U) = 

Jl- l ^'{U) - (d^JI~\^(U),0) e X = H 1 x L 2 (see the regularity shown for 7' in Lemma 4.3), so that 
integration by parts yields 

(7 / (*hyW),9t*h y (t))x % x = -(^, + ^ ) (*h y (t)),A^(* hy (t)))x.,x. 
Inserting these relations into (25) then gives 

|n K (*h y (t)) = -(^ c . + ^)(*h y (t)), 

{ JT K ( 7 (i))(. + y(t)) + (*hy(*), 9 7 T«(7(*))(- + y(t)))^(*h y (i))}) xt x , 

which is the desired equality. 

If *h y € C°([0, T), A) is just a continuous in time solution to (15) that remains in G e , then the integrated 
relation 

t 



fi K (*hy(*)) = ^«(*L y ) + / S K (* hy (r)) dr 

Jo 



holds, as can be seen by using the continuity of the flow and the approximation of such a solution by smoother 
ones (see [24]). □ 

We now compute the asymptotics of E K (W(-f)) for 7 — > and small k. 
Lemma 4.9 We have 

5 K (W( 7 )) = ~^Y~ + «(7^(o,o)(7 2 )- 

Proof. The proof follows the one of Lemma 5 in [50]. As a first step, notice that 7(W(7)) = 7, J/(W(7)) = 0, 
as it can be seen from the equality G(W(j), 0, 0, 0) =0 and the local uniqueness of the solution to G = 0. 
Therefore, since P Ct+1 (U c , +1 ) = and cr( 7 ) - - 7 2 iy (2||Z4 1^), 

(W( 7 )) = J^ +7 (W c , +7 + crfrJBZ^) = a( 7 )^', +7 (Z4, +7 )[BW Ct+7 ] + o 7 ^ (7 2 ) 



7^* 



J-:'(W,)[BW,]+o 7 ^ (7 2 ). 



37 



In addition, since U c is even and T K (7) is odd, we deduce 

(W( 7 )(- + y(W( 7 ))), 9 7 T K ( 7 )) H = (U c , +1 + a( 7 )K4, +7 , 9 7 T K ( 7 )) H = 0. 

Consequently, 

^(W(7))= 9 &(^ / (^)[B^],JT K ( 7 )) X ,, X + O7 ^o( 7 2 ), 
^11"* II j? 

where "o 7 _^ ( 7 2 )" does not depend on k. Moreover, Lemma 4.7 provides ||JT K ( 7 ) — [d c U^\\ c=Ct+1 \x < k 
independently of 7 e (—71, +71), hence 

Sk(W(7)) = 2^ {T " {U * )[m * ] > [W|c=c, +7 >**,* +O(7,k)^(0,0)( 7 2 ) 

= 2^ raW * )pBW * ] ' [Wl— )^*.^ +°(7,«)-(o,o)(7 2 )- 

Finally, using once again the equality (for e X) (J 7 * (W*)[0cW c ]| c =c» , (f>)x»,x = (BW*,0)h and that J 7 " is 
self-adjoint, we infer 

(^ , (W.)pBW.],[a c W c ]| c=c .)x.,x = (^ , (W.)([5 Wc]|c=c.),BW.) x . lX - |BW.|^ - 
and reporting this into the previous expression gives the result. □ 

We now compute the asymptotics of S K for more general functions. 

Lemma 4.10 Let e > be small enough. IfU € G e satisfies Ph y (U) = Phy(U*) and F*{U) — T*{U*) < 0, 
then, we have 

uniformly for < K < \j(U)\ 3 . 

Proof. First, we may apply Lemma 4.6 and infer that + |a(W)| = 0(\j 3 (U)\). Then, we write 

E K (U) = E K (U(--y(U))) =E K (W^(U))+^(U) + a(U)MU c ,+^ u) ) = S K (W( 7 (W)) + i?(W)) +0(| 7 (W)| 3 ) and, 
recalling the expression 

3«(W) = -(K-n(u)M> {^ K (7(W))(- + y(U)) + («> 9 7 T K (7(W))(- + m))H^(U)}) x ,y 
we wish to make a Taylor expansion. First, note that 

^c.+7(H)(«)=^(K)(M7^^ 

= ^» + ^)(W(7(W)))+^: +7W (Wc, +7 (^))[i?(W)]+0(|7^ 
hence, since J :! Cit+ ^iU Cti+ ^) = and (Lemma 4.1) (7(7) = C(7 2 ), we have W( 7 ) = W Ct+7 + C(7 2 ), thus 
S K (W)-S K (W( 7 (W)))=0(|7(W)| 3 ) 

- {K+j(u)(Uc, + ^u)W(U)}, { JT k (7(W))(- + y(W)) + (W, 9 7 T K (7(W))(- + y(W))) H ^(W)}) x> ^ . 

Now, in the bracket term, we may replace U by W( 7 (W)) + (since \${U)\ 2 X = C(| 7 (W)| 3 )). By 

the computations of Lemma 4.9 and the equalities 7(W(7)) = 7, j/(W(7)) = 0, this gives 

H K (W) -H K (W( 7 (W))) =0(|7(")| 3 ) - (^'. +7 ( M )(W C . +7 ( M ))[^(W)], JT K (m)) x , x 

=oqj(u)\ 3 ) (/::„ M (w c , + ^ ) )^(M)] 1 m}| C=c , +#) ) xt ^+ow 

= 0(| 7 (W)| 3 ) - (^^ (M )(W c .-Hr( M ))[{9cZ4}|c=c.-HT(M)]^(W)) jc . x +©(«) 



; J ,8 



using Lemma 4.7 and the self-adjointness of P" + ^/ U )- Choosing < k < ^(U)] 3 and from the equality (for 

4> e X) (J*{U c )[dJ4c],<l>)x;x = (BZ4» H , we infer 

E K (U) - S„(W(7(W))) = G(|7(^)| 3 ) - {WA C ^ (U) ,${U)) H = 0(|7(W)| 3 ), 

by the orthogonality condition in Lemma 4.4. Inserting the expansion of S K (W(7)) given in Lemma 4.9 
yields the conclusion. □ 

Proof of Theorem 5. We have to show that there exists e > such that, for any S > 0, we can choose 
an initial datum at distance < 5 from U* but that escape from G e . Since W(j) — > U* in X, we shall take 
the initial datum to be W( 7 ) for some small 7, and denote ^hy(i) the corresponding solution. In view of 
Lemma 4.2, we have P*(W(7)) — F*(U*) <~ — 7 3 P*/6, hence we can choose 7 with the sign of P* 7^ so that 

J-,(W(7))-^(W,)^-| 7 | 3 |P >t |/6<0. 

We now assume that ^hy(t) is globally defined and remains in ff e , where e is as in Lemma 4.6. By conservation 
of energy and momentum and the construction of W(7), we deduce Phy(^'hy(i)) = -fV(W( 7 )) = Ph y (£4), 
and j;(\I> hy (t)) - T*{U*) = J"*(W(7)) - T*(U*) < 0. The first step is to have a control on j(t) = 7(*h y (i))- 
We denote a(t) = a($> hy (t)), = y(* hy (i)) and «?(i) = ?9(\I> hy (f)). Applying Lemma 4.6, we obtain 
+ l a (*)l = C(l7 3 (*)D- I n addition, Lemma 4.2 and Taylor expansion gives 

•F.(*hy(*))- P*(Z4) = ^*(W(7(i) + 0(t) + a(t)MU Ct+m ) - P*(Z4) 

= K(WW)) + ^(W(7(t)),^(t))x.,x + ^(^'(W(7W)^W^W)x.,x + (| 7 3 WI) 

= - + (X{ww)),m)x;x + l(F':(wm))m,m)x*,x+o(\f(t)\), 

where we have used that J r ^(W(7(t)) = o(l) (for the terms involving a(t)) and Lemma 4.2. Furthermore, 
by the orthogonality relations in Lemma 4.3 and using that 17(7) = 0(j 2 ) and F' C {U C ) = 0, it holds 

^(W(7(^)),#)}x^x={^ +^ ( t )(W(7(^)),tf(i)}x^x+7(^)(BW(7(^)),^))^ 

= <^.-Hy( t )(«c.-Hy(t) + CT (7W))BW C .+7(t)^(*)>X-,X =0(|7 7/2 W|). 

In addition, by Lemma 4.4, the before last term is > Xo|i?(t)|^-/2. As a consequence, by conservation of 
^(^hyit)), we infer, for small 7, 

> -|7| 3 |^l/3 > ^*(W(7)) - = ^,(*hy(*)) - ^.(W.) > +o(|7 3 (;)|), 

In particular, this forces j(t) to always be of the sign of P* and to satisfy \j(t)\ > | 7 |/2 (provided e and 7 
are small enough). 

Since now, we have a good upper bound for | 7 (i)|, we can choose k — n(-f) = 7 3 /8, which is such that, 
for any t > 0, k < |7(i)| 3 - In particular, we can apply Lemma 4.10 and get 

S K (* hy (t)) = -^^ + (7W 2 ). 
With this choice k = k{j), we deduce from Lemma 4.8 that 

-o K(7)( * hyW) = s K (* hy (t)) = + omf). 

Since \j(t)\ > | 7 |/2, it follows that, when P* < (the case P* > is analogous), 

|o K(7) (* hy (i))>-°^>0, 

hence fl K ( 7 )(\&hy(£)) i s unbounded as i goes to +00. However, by definition of Q, K , we have by Cauchy- 
Schwarz |Q re ( 7 )(W)| < IW||jy|T K ( 7 )|/f < C(j) for W G We have reached a contradiction. The proof of 
Theorem 5 is complete. □ 
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5 The linear instability (0 < c* < c s ) 
5.1 Proof of Theorem 4 

Existence of at least one unstable eigenvalue. The proof of the existence of at least one unstable 
eigenvalue relies on Evans function technique, as in [70], [6]. We shall actually use Theorem 1 in [6] when 
observing (see e.g. [5]) that the Euler-Korteweg system 



( d t p + 2d x {pu) =0 

8 t u + 2ud x u - d x (.f(p)) - d x (K{p)d 2 xP + ^K'{p)(d x pf) = 0, 



(EK) 



where K : (0, +oo) — > (0, +oo) is the (smooth enough) capillarity, reduces to (2) (where, we recall, VP = Ae 1 ^, 
p = A 2 and u — d x (f>), namely 

( d t p + 2d x (pu)=0 



d t u + 2ud x u - d x (f(p)) - d x (^@-) = 0, 



for the capillarity K(g) — l/(2g), as can be shown by straightforward computations. The associated eigen- 
value problem in the moving frame is 







f \(-c*d x ( + 2d x {{rl+ri*)v + (u*) 

Xv - c*d x v + 2d x (u*v) - d x (f'(r 2 + r/*)C) - 



2 \Ao +V* X ^ V r o + V* 
The link with the original eigenvalue problem (6) is done through the formula 

w = U.(±+i£ v), 



2(^ + ^)3/2 



0. 



(27) 



(28) 



since this corresponds to \& = U Cli + ijj = U Cr + e xt w(x) = (A Ct + e At £(x)) exp(i(f> Cr + ie xt J_ v). Notice 
indeed that the second equation in (27) gives J R u dx — 0. It then follows from Theorem 1 in [6] that under 

the assumption ^P(Uc) > ^ th ere exists at least one unstable eigenvalue 70 € (0, +00). 

dc \c=c. 

Existence of at most one unstable eigenvalue. The fact that there exists at most one unstable 
eigenvalue follows from arguments as in [7] (Appendix B) and is a direct consequence of Theorem 3.1 in [57], 
that we recall now. 

Theorem 11 ([57]) Let J and C be two two operators on a real Hilbert space X, with C self-adjoint and J 
skew-symmetric. Then, the number of eigenvalues, counting algebraic multiplicities, of[JC]c in the right-half 
plane {Re > 0} is less than or equal to the number of negative eigenvalues of C, counting multiplicities. 

In order to apply this result to our problem, let us write the eigenvalue problem (27) under the form 



A 



1 

1 



C 



where M is the self-adjoint Sturm-Liouville operator 

1 



M 



d 2 x {M+V*) 

S^T^T^y^T^y ' 2(r 2 +^)3/2 



'Or? +„. ) 
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(which is bounded from below) on W = L 2 x L 2 and with 



L = 



M 2m* - 
2u*-c* 2(r§ + 77*) 



We are in the setting of Theorem 11 with J — —d x y ^ ^ J skew-symmetric and C self-adjoint. We thus 

show that L has at most one negative eigenvalue. Since rg + r/* remains bounded away from zero, it is clear 
that, for a < and ((, v) given, £((, v) 1 — a((, v) f if and only if 

M^C- ( o C ;; 2 "*) 2 ■ a C = < with M^M- {C 2~ 2U * ] 1 , (29) 

since we may express w in terms of Q with the second equation. We observe that the translation invariance 
shows that d x (r]*, u*)* belongs to the kernel of C, that is, using once again the relation 2u c = 2d x 4> c = c ^ 2 , 

M^d x r]* = 0. Furthermore, has the same continuous spectrum as its constant coefficient limit as 
x — »■ ±00, namely 

1 c 2 - r 2 

2r 2 ^ + 2r 2 ' 

that is a ess {M^) = [c 2 — c 2 , +00) C (0, +00), since < c* < c s . Since 9 x r/* has exactly one zero (at x = 0), 
it follows from standard Sturm-Liouville theory that Ai^ has precisely one negative eigenvalue n < and 
that the second eigenvalue is 0. Taking the scalar product with (29) yields 



2 



Now, for s < 0, we consider the self-adjoint operator 



2(r 2 + ^) 2(r 2 +^)- S " 

Clearly, M' s=0 = M\ <r css (Ml) C [c 2 - c 2 ,+oo) C (0,+oo), and R_ 3 s i-> .Mj is decreasing. Let 
us assume now that the self-adjoint operator £ has at least two negative eigenvalues. Then, we denote 
<j\ < (72 < the two smallest eigenvalues of C (necessarily simple), and d, £ 2 two associated eigenvectors. 
Since C is self-adjoint, (Ci, C2)l 2 = 0. Furthermore, ( J =<T Ci > Ci ) £ 2 = cil^l^ < 0, hence, by monotonicity, 
(MIC2, C2) l 2 < for any a\ < s < 0. Therefore, .MJ has at least one negative eigenvalue for o\ < s < 0. 
We denote A m ; n (s) the smallest eigenvalue of Then, A m ; n (s = 0) = [i < and A m i n decreases in 

[<7i,0]. Moreover, we may choose a positive eigenvector £ s for the eigenvalue Ai(s), with £i = Since 
a{M)) n R + = {yu, 0}, it follows from the monotonicity that for any a x < s < 0, we have a(M\) H R_ = 
{A m i n (s)}. When s = er 2 € (01, 0), we then have 2 <G a(A4\ =a2 ) flK_, and thus <t 2 = A m j n (cr 2 ), which 
implies that we may choose C2 > without loss of generality. Similarly, if s = 02, we see that we may choose 
C2 > 0. We obtain a contradiction since then (Ci>C2)l 2 > and thus £i and £2 cannot be orthogonal in L 2 . 
We have thus shown that C has at most one negative eigenvalue, and then Theorem 11 shows that JC has 
at most one eigenvalue in {Re > 0}, as wished. 

5.2 Resolvent and semigroup estimates (proof of Corollary 3) 

In this section, we drop the "*" for the travelling wave we are considering. When linearizing the (NLS) 
equation in the moving frame with speed c, we obtain 

- icd x ^ + d 2 x iP + V/(|C/| 2 ) + 2(1>, U)f'(\U\ 2 )U = 0, (30) 
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or 

cd *- 2/'(|C/| 2 )^i(72 -dl f(\U\ 2 ) V'(\U\ 2 )Ui \ ( ^ 

at \<h ) ~ V d t + f(\ u \ 2 ) + V'(\u\ 2 )u? cd x + 2f(\u\ 2 )u 1 u 2 

/ l\f-d 2 -f(\U\ 2 )-2f(\U\ 2 )U 2 -cd x -2f(\U\ 2 )U 1 U 2 WVi 
-^-10^ cd x - 2f{\U\ 2 )U 1 U 2 -dl f{\U\ 2 ) 2f(\U\ 2 )U 2 J 1, V2 

We wish to show that this linear equation can be solved using a continuous semigroup. In order to handle 
later the nonlinear terms, we work in H 1 (R, C 2 ) instead of L 2 (R, C 2 ). Therefore, we consider the unbounded 
operator A : D(A) = H 3 (R, C 2 ) C if^R.C 2 ) -> ^(R.C 2 ) on ^(R.C 2 ) defined by 

A= ( cd x - 2f(\U\ 2 )U 1 U 2 -dl f(\U\ 2 ) - 2f{\U\ 2 )U 2 \ 
A -{dl + f(\U\ 2 ) + 2f'(\U\ 2 )U 2 cd x + 2f(\U\ 2 )U 1 U 2 J ■ 

It follows easily that for V = ^ ^ ^ e ^(R.C 2 ), 

Re^VMHW 2 )) =Rc((-2/'(|(7| 2 )[/ 1 [/ 2 ^ 1 ,^ 1 )h 1 ( R ,c) + (-2/'(|C/| 2 )C/ 1 C/ 2 V2>2) ffl(R ,c) 

+ ([/(|t/| 2 ) + 2/(|C/| 2 )[/ 1 2 ]^ 1 ^2> ffl (K,c) - ([/(l^l 2 ) + 2f(\U\ 2 )U 2 ]iP 2 , Vi) H i (R ,c)) 

<-^l^llffi(M,C2)- 

Moreover, the spectrum of .4 is included in the half-space {Re < cto}, hence A generates a continuous 
semigroup e tA on H 1 (R, C 2 ). 

In order to estimate the growth of the semigroup e tA on C 2 ), we could try to use the same approach 

as [21], which relics on the proof of the spectral mapping theorem in [31]. However, our situation is slightly 
different since in these studies, the reference solutions is real-valued (it is a bound state in [31] and the kink 
in [21]). Therefore, U 2 = and A has no diagonal term, and the system is much more decoupled than in our 
situation. As a matter of fact, it is not very clear whether the arguments of [31] carry over to our problem. We 
thus have chosen to use the general approach given in Appendix B. We thus verify the assumptions of Theorem 
B.4 (see also Corollary B.2) there, which are easy: A generates a semigroup in H 1 (R,C 2 ) and the spectrum 
of A is of the form iR U {— 7o,+7o}> where iR is the essential spectrum and ±70 two simple eigenvalues. 
Moreover, the eigenvector associated with 70 belongs to H 3 (R,C 2 ) = D(J). Therefore, Theorem B.4 in 
Appendix B applies and the growth estimate for the linearized problem follows. For the nonlinear instability 
result, we argue as for Corollary B.2 in Appendix B, since the manifold 9Jt = {|f/*|(- — y), y e R} is transverse 
to the curve u i-> \U*+crw\ in r a + H 1 (R). Indeed, it follows from (28) that \U* +crw\ = A + o"C + Off 1 (c 2 )- 
Assume that ( = ad x \U*\, with a£l. Then, integration of the first equation of (27) provides 

\{\UJ- r ) - c*d x \U*\ + 2((r 2 + r]*)v + u*c) = 0, 
hence, using that \U*\ = \Jr 2 + 77* and the equality 2u* — cr\„j(r 2 +r]*), we infer 

, + a{A^^ + ^ 2 C ; r ° W2 ^}=0. 
I r 2 + 77* 4(r 2 + r?*) 3 / 2 J 

Since J R v — and \U t \ — r has constant sign in R, integrating over R then implies a = 0, which in turn yields 
( = v = and = 0, a contradiction. Consequently, ( £ M.d x \U* | and the manifold DJl = {|£7* | ( — y), i/Gi} 
is indeed transverse to the curve a i-> [£/* + uw\ in r + if 1 (M). 
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6 Stability analysis for the kink (c = 0) 
6.1 Proof of Lemma 4 

Let us recall that the momentum P(U C ), for c > 0, has the expression 



P{U C ) 



since sgn(£ c ) = —1. Therefore, we decompose P(U C ) with two integrals: 



£2 



1 



Using the change of variables £ = f£ c , the second integral in (31) is equal to 

1 



di. (31) 



Ji rl + t£ c \ y=v 



i rg + t^c V-Vc(^c) V-£ c V£(£ c )(i-l) 



) 



eft 



2n3 



(-r 2 ) 

r0 



1 '0 



.( L_ 



oV-VoK) v/-4r 2 F(0)(t-l) 



rfC + °c^o(l)- 



The passage to the limit c — > being justified by the dominated convergence theorem since the absolute value 
of the integrand is < Kt for < t < 1/2 for small c and for 1/2 < t < 1, since £ c > — r 2 , r 2 

+ *& > 4(1 -t) 

and hence is equal to 

t 2 



V c (t&) - < 


;c^c)(*-i) 


V / -Vc(^c) V / -^(e c )(t-l) 


y-v c (te c ) + v/-^(e c )(t-i) 



if 



€^((1/2,1)). 



(1 - t)VT=ty/T=ty/T=t y/\ - t 

Furthermore, letting £ = £ c + (r% + ^ c )t 2 , i > 0, the first integral in (31) is equal to 
1 (& + (r 2 + £ c )i 2 ) 2 2di 



Jo 



+ (r 2 +t c f 
by direct computation. Since £ c ~ 

„2„4 

= -rl + 

thus 



1+t 2 



^{^[|-arctan(^)]-2, 2 (, 2 + U 



4+tc 



+ - 



n + 3 v r 2 + e c 



]}■ 



r 2 , is a simple zero of V c (£) = c 2 £ 2 — 4(t-q + i)F(rl + £), we have 



o 4 + c 4 r 6 ro 2 /(0) 



-2) + o c ^ (c 4 ) = -r 2 + 



4F(0) 4F(0) V F(0) 

-^(e c )=4F(0)+O c ^o(c 2 ) 
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and 

2 2 
As a consequence, the first integral in (31) is equal to 

c I ^WW) v^W 3</F(6) J c 3v /F(0) 

Gathering these two relations, we obtain 

as wished. 

6.2 Proof of Theorem 9 

Since we have a kink solution U~o for c = 0, this implies that Vo(£) = — 4(r§ + 0-^'( r o + £) ' s negative in 

2 

(-rg, 0) and that -r 2 is a simple zero of V , that is F(0) > 0. Then, F > in [0, rjj) and F(g) ~ ^-(g-r 2 ) 2 
for — > r 2 , it follows that there exists K > such that 

FG>)>-^-r 2 ) 2 

AO 

We consider for /x > the quantity 

^nin(M) = ml w e Z, inf|u|=/i|. 

The study of J^; n (0) is easy. 
Proposition 6.1 TTiere ZioZds 

^nin(O) =E(U Q ). 

More precisely, for any U E Z, 

pro pro 

E(U) > 4 / y/F(s 2 ) ds and E(U ) = 4 / y/F(s 2 ) ds. 

Jini K \U\ JO 

Finally, if U e Z, infR |£7| = and Jf(U) = E(Uq), then there exists y e R and 9 e M suc/i i/iat 
U = e i9 U (--y). 

Proof. Taking J7 as a comparison map, we see that J^i n (0) < E(Uo)- Moreover, if U e Z and inf R |[/| = 
/U > 0, we may assume, up to a translation, that \i — |J7|(0). Then, denoting 



G(r) = 2 f v/^V) ds, 



we have the inequalities 

r>+oo /*+oo Z'+OO 



/—too r-\-oo p-\-oo 

/ |9 x ;7| 2 + F(|;7| 2 )dx> / |^|<7|| 2 + F(|£/| 2 ) dx>2 / |/FW)^|C/||rfx 
Jo Jo Jo 

;-+oo /•-{-CO 

= \d x [G(\U\)]\dx> / da; = G(\U\(+<x>)) - G(\U\(0)) 

Jo Jo 

pro 

G(r ) - G{n) = 2 / T^V) ds. 

J (J, 
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Arguing similarly in (— oo,0), we get 

pro 

E(U) > 4 / y/F(s 2 ) ds. 
J fi 

For the kink Uo, which is real-valued, we have the first integral |<9a;E/o| 2 = F(Uq), hence, using the change of 
variables s — Uq(x), 

i> + oo 



E(U ) = 4 / F(Ui) dx = A VF(s 2 ) ds. 
Jo Jo 



If (j, = 0, we have then E(U) > E(U ), hence JfT(U) > E(U) > E(U Q ) as wished. 

Assume finally that U € Z verifies infR \ U\ — and Jf(U) — E{Uq). Then fi = and all the above 
inequalities are equalities. In particular, we must have \d x U\ — \d x \U\ | and equality in | 2 + F(\U\ 2 ) > 

2\y/ F(\U\ 2 )d x \U\ |, which means that \d x \U\\ = V / ^WF)- Combining this ode with the condition \U\(0) = 0, 
we see that \U\ = \Uo\, since \Uo\ solves d x Uo = y/F(U^). Finally, the fact that \d x U\ = \d x \U\\ implies 
that the phase is constant in (— oo,0) and in (0,+oo): there exist two constants 6± e R satisfying U(x) = 
e i9 ±\U \(x) for ±x > 0. Therefore, <£([/) = r 2 {9+ - 6L) mod 2irr 2 , and then 

E(U ) = Jt{u) = E(U ) + 2Mr 4 sin 2 ( ? + ~ 9 ~~' K ) 
implies 9 + — 9_ = n mod 2ir, that is U = e ie +U in R, which is the desired result. □ 

We recall the expansion P(U S ) = r 2 ir + sP a + o(s) as s — » 0, where P = ^^ s ^ . From the Hamilton 

as \s=o 

group relation j^i^j = g dP(U s ) ^ wg a j gQ - m f eJ . integration E(U S ) — E(U ) + — Po + (s 2 )- As a first 

ds ds 2 

step, we define the small parameter yu* > 0. The key point is to prove the following result. 

Proposition 6.2 There exist some constant K > and a small /i* > such that, for any < jj, < /j,*, 



^„i„(M) = inf {jf(U), U eZ, inf |u| = /i j > E(U Q ) + 



K' 



Proof. Notice first that for c > small, there exists U c travelling wave of speed c and that infjj |I7 C | = yV 2 + £ c 

2 4 2 

with £ c a smooth function in c such that £ c = — r 2 H -A- + O(o ), hence inf R |t/ c | = — + C(c 2 ) and 

4i? (°) 2 V F (°) 

is smooth. Therefore, there exists, for < \i < ^ small, a unique er^, with = 2 — + 0(fi 2 ), such 

that /U = infjR |f7<T |. In particular, taking U a as a comparison map in J£ni n (/u), we have 



X.inM <^(^J = S(I^J + 2Mr 4 sin 2 ^ ) 

= E(U ) + f Po + o(a 2 ) + 2M^sin 2 ( ^^M ) 



^o) + ^(p + Mi b o 2 )+«(^)- 



In particular, it follows that, for some positive constant K and for /i* small enough, 

J&un(/i) < S(C^,) + Kfi 2 < ^E(U Q ). (32) 
Consider now c small, a bounded open interval (x-,x+) and r\ a solution to the Newton equation 

2a 2 r? + v'M = o 
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in (x-,x+), with d x r)(x+) < < d x r q{x^) 1 r/(x + ) < — Tq + /i 2 and rj(x-) < — r§ + /i 2 . As c — > 0, V c converges 
to Vo in C 1 ([— r 2 ,0]). Moreover, Vo is negative in (— r 2 ,0) and has a simple zero at — r 2 . Therefore, if c and 
[i* > arc sufficiently small, we must have F(tq + r]) dx > \ J M F(U 2 ) dx. Consequently, if v = Ae lv 

solves (TW C ) on a bounded interval x + ), satisfies 2d x (p = —r {q = A 2 — r 2 .) and if Id is < /li* at x + 

and at x_, with d x \v\{x + ) < < d x \v\(x-), then 



J X- 



+ \8 x v\ 2 + F(\v\ 2 ) dx> \e{U ). (33) 



Here, we use that the Newton equation on the modulus \V\ actually holds true in {x-, x + ). Since F > in 
[0, r 2 ) and F(g) ~ rl(g - r 2 ) 2 when g -> r 2 , there exists if > and /c > such that F(g) > (g - r 2 ,) 2 /K for 
< g < Tq(1 + k) 2 . Hence, if inf R \v\ > fi > 0, then 

\P(v)\ < -E(v). (34) 

Moreover, arguing as in the proof of Proposition 6.1, we show that there exists x > such that if U £ Z 
and | U | takes values < /U* and > r (l + k), then 

E(U) > E(U )(l + x). 

In particular, since J^m n (M) < E(U~o) + 0(/i 2 ), we may choose /i* sufficiently small so that if U E Z and 
J?T(i7) < ^ni n (/x) + /it*, then |{/| < r (l + ft). This means that for the mappings we are considering, 
F(g)>(g-r 2 ) 2 /K. 

Step 1: Construction of a suitable minimizing sequence. There exists a sequence (V n ) n >o in Z such 
that inf K \V n \ = m = |K|(0), V n = A n e^, P{V n ) e [0,7rrg], 

2A 2 n d x K = c n {A 2 n - r 2 ), c n ee Mrl sin ~ r g (K) ) > 

and 

lim JT(K) = ^nin(A*)- 
n— >+oo 

Since /U > 0, the maps V we consider may be lifted V = Ae 1 ^ . Therefore (with u — d x (j)), 
J^M = inf { / (^ A) 2 + F(A 2 ) dx 

+ inf { J A 2 u 2 dx + 2Mr 4 sin 2 ~ r l>^ ~ , u e ^(R, R)}, 

A e r + iJ 1 (M,M),inf A = jit). 
The infimum in u may be written 

inf i°f { / dx + 2Mr 4 sin 2 ( P ~ , u € L 2 (R, R) a.t. f (A 2 - r 2 )u dx = p}. 

For each p € R, we minimize in u a quadratic functional on an affine hyperplane, with minimizer given by 

/ f (A 2 -r 2 ) 2 N-i^-rg 

Up=p {L— ^ dx ) 



(35) 
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As a consequence, the infimum in u in (35) is 

r (A 2 - r 2 ) 2 A \-i , 4,, , /p - 7rr 2 



inf 

PGR 



/ A 2 U 2 dx + 2M7>sin 2 = inf [p 2 ( / (A J q) dx) * + 2Mr 4 sin : 



2rg 



It is clear that this last infimum is achieved only for p inside [— irr 2 ,, +irr 2 ]. Indeed, the second term is 

2.r 2 -pcriodic and if p > ,r 2 , then p- 2nr 2 is a better competitor. Moreover, the function p * sin 2 (^f) 

is continuous and even, hence we may consider some p € [0, 7rr 2 ] (depending on A), which is a minimizer for 
this last infimum. The corresponding u p is then a minimizer for the infimum in u in (35). Writing that 



dp 

we deduce the relations 



2A 2 u = c(A 2 -r 2 ), c= Mr 2 sin 



2r 2 



We conclude by considering a minimizing sequence (A n ) in (35), and translating in space so that inf^An = 

m = KI(o). 



Since F > in R + , we have 



/ 



d x V n \ 2 dx<JfT{V n )<^E{U ) 



for n large. Therefore, by compact Sobolev embedding H 1 ([—R, +R]) i?, +-R]), we may assume, up 

to a possible subsequence, that there exists V € ifj^R) such that for any R > 0, 14 — 1 V in iJ x ([— _R, +i?]) 
and — > V' uniformly on [—R,+R]. Moreover, by lower scmicontinuity and Fatou's lemma, E(V) < 
lim„^. +oc E(V n ). Since \V n \ > fi > in E, |V| > /x > in R and we may lift V = Ae'^. Furthermore, 
inf R A„ = fi = |V^|(0), hence inf R ^4 = /i = |V|(0). We also know that P (V n ) £ [0,r 2 7r] for all n, hence we 
may assume, up to another subsequence, that P{V n ) converges to some Poo € [0, r^n]. We also set 

c = JR. Cn = Mrl sin (^27p) ■ 

In view of Step 1, and the convergence A n — > A uniformly on any compact interval [-R, +R], it follows that 
2A 2 d x <t> = c(A 2 - r 2 ) and d x <j> n ^ d x (j> in (36) 

Note that 



/ 



+ ^(|F| 2 - r 2 ) 2 < £(V) < +oo 



hence |V"| — > ro at ±oo. In particular, there exist — oo < i?_ < < J?+ < +oo such that |V| > /i in 
(—oo, and in +oo) and |V|(i2±) = /i. 

Step 2. There exist — oo < z_ < < z + < +oo such that 

A(x) = A c (x - R + + z + ) for x > R + and A(x) = A c (x - i?_ + zJ) for x < 

We work for x > the other case being similar. We consider \ € C^((i? + , +oo), C), i e K small 

such that V* = v n + t\ verifies \V*\ > \i in (i? + ,+oo). This is possible since inf Supp ( x ) \y n \ > yu. Then, 
\V*\ > \x in R and |V^|(0) = /j,, hence V£ is then a comparison map for J^i n (/x), and in view of the equality 
P{V*) = P(V n ) + 2t f£°(id x V n \x) dx + 0(t 2 ), it follows that 

r+oo r+co 

<XnM<^(V, t l )=<je m M + o n ^ +00 (l) + 2t (d x V n ,d xX ) dx + t 2 \d xX \ 2 dx 



f + °° / P(V ) - -Kr 2 \ f + °° 

-2t f(\V n \ 2 )(V n , X ) dar + Mt sin (-^2 °) (id x V n , X ) dx + O 

Jr+ v r o ' Jr+ 



t^o(t 2 ). 
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Letting n — > +00 and using the weak and strong convergences for V n , we infer 

»+oo ,.+00 . 2 p . r+co 

(d x V,d xX ) dx-2t f(\V\ 2 )(V, X ) dx-Mtsin( 00 j / (i^y, X ) <fc + O t ^{t 2 ). 



0<2t 



Dividing by t 7^ and letting t — > + and then t — > 0~, we deduce that V solves (TW C ) in (-R+, +00) and 
V has finite energy. Moreover, |V|(i? + ) = /1 is small, thus V — e ie +U c (- — R + + z + ) in (R + ,+oo) for some 
constants z + and 9+, and the speed c is such that ini^Ac — V'o+Cc < Mi hence c < er(/i) < if^i. Since 
|V| has finite energy in R and solves (TW C ) in (i? + ,+oo), V is C 1 in [i? + ,+oo). Moreover, \V\ reaches a 
minimum at x = R + , thus we must have d x \V\(R + ) > 0, which imposes z + > 0. Note that A c being even, 
it is possible to translate V so that R = R + = — J?_ and z = z + = — z_. Observe that \x = A c {z) > A (z), 
hence z < K^i. This yields 

/ \d x V\ 2 + F(\V\ 2 )dx= [ \d x lI c \ 2 + F(\U c \ 2 ) dx>E(U )-Kn. (37) 

J\x\>R J\x\>z 

In particular, we deduce from (32) 

2i?F( M 2 ) < I \d x V\ 2 + F(\V\ 2 ) dx < Kfi, 

J\x\<R 

hence R<K[i for fi small (F(0) > 0). 

Step 3. We prove that A = (i in (R_, R + ) = (-R, +R). 

Indeed, if it is not the case, there exists a bounded interval (x-,x + ) such that A = \V\ > [i in (x-,x + ) 
and I V|(ar±) = fi, with 9 x |F|(a; + ) < < <9 x |F|(a;_). Therefore, we can make perturbations of the amplitude 
A n localized in (x_,x + ). Hence, arguing as in Step 2, we see that then, V solves (TW C ) in (x_,x + ), with 
2A 2 d x (f> = c(A 2 - rl) and |V|(ar±) = fj,, d x \V\(x+) < < d x \V\(x-). We then are in position to apply (33), 
yielding 

d x V\ 2 + F(\V\ 2 ) dx>'-E{U ), 



L 



2 

but the combination with (37) provides 

^E(U ) > JT min ( M ) > [ X+ \d x V\ 2 + F(\V\ 2 ) dx+ f \d x V\ 2 + F(\V\ 2 ) dx 

iU Jx- J\x\>R 

> ^E(Uo) + E(U ) - K(i* = \e{U ) - K^, 

which is not possible if jU* is sufficiently small. 

Step 4. We have R = or (z = and c = a^). 

Indeed, assume R > 0, and consider C € C*((0, +00), R), C > 0, t > and V* = (A n + tQe^, so that 
\V*\ = A n + t( > fj, in R. Since R > 0, we actually have infjR \ V%\ = fi and V* is a comparison map for 
J^min(^)- Arguing as before, we thus have 

r+00 r+00 

Xninifi) < ^(V*) =Jtr m M + 0„^ +(x) (l) + 2t / O x A n d x ( dx + t 2 (d x () 2 dx 

Jo Jo 

/*+oo /"+C30 /"+oo 

+ 2t / A n ad x <t> n f dx + t 2 C 2 {d x <t> n ) 2 dx-2t f(A 2 n )A n ( dx 
Jo Jo Jo 

+ Mr 2 tsin( P(K) ~ r ° 27r ) \ R 2A n C,dA n dx + O^i 2 ). 
v r o 'Jo 
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By (36), we may pass to the limit as n — >• +00 in all the terms and deduce 



/■+00 /" + 00 /" + 00 /"+oo 

< 2i / d^d^C da; + 2i / ^C(^0) 2 dx - 2t / /(A 2 )AC dx - 2ct / ACd x </> dx + O t ^ 
Jo Jo Jo Jo 



(i 2 ). 



At this stage, we see the relevance of taking a minimizing sequence as chosen in Step 1, since it allows to 
pass to the limit in the nonlinear terms involving d x (f> n . As a consequence, using (36), 



-d 2 x A-Af(A 2 ) 



c 2 (A 2 - r 2 ) 2 



> 



c 2 rjp _ r 2\2 

in the distributional sense in (0, +oo). The term — Af(A 2 ) + — — is continuous in R. However, 

since A{x) = \x for < x < R and A(x) = A c (x — R + z) for x > R, we infer —d 2 A = —d x A c (z)5 x= R plus a 
piecewise continuous function in the distributional sense in (0, +oo). Since d x A c (z) > (recall that z > 0), 
this forces to have d x A c (z) — 0, that is z = 0. Consequently, fi — |V|(i2) = A(R) = A c (z) = A C {Q) and then 



In the next step, we take into account the loss in the weak convergence V n — V. 



Step 5. There exists K > such that 

\>. 1.KTO h t = 

n— >+oo 

Let e > be fixed but small, and pick some A > large so that 



E t >^, where E t = lim E{V n ) - E{V) > 0, P s = lim P(V„) - = - 



E(V)- f \d x V\ 2 + F(\V\ 2 ) dx <e P(V)- f (A 2 - r 2 )u dx 

J\x\<X J\x\<X 



< e. 



We claim that there exists some small fi > 0, independent of e such that \V n \ > /i for |x| > A and n 
large. Indeed, otherwise, we may argue as in Step 3 and show, as in the beginning of the proof there, that 
S\x\>x \ d x V n\ 2 + F (\ v n\ 2 ) dx > \E(U Q ). This is not possible since 



12 
IT 



W ) > lim E(V n ) > \E{U ) + f \d x V\ 2 + F(\V\ 2 ) dx > \e{U ) + E(V) - e, 

I 7i->+oo £ J\x\<X 1 

and E(V) is close to E(U ) as fi -> 0. Therefore, as for (34), 

/ (A 2 n - r 2 )u n dx < * [ \8 x V n \ 2 + F(\V n \ 2 ) dx. 

J\x\>X H J\x\>X 

Consequently, 

E(V n ) - E(V) > f \d x V n \ 2 + F(\V n \ 2 ) dx- f \d x V\ 2 + F(\V\ 2 ) dx + f \d x V n \ 2 + F(\V n \ 2 ) dx - e 

J\x\<X J\x\<X J\x\>X 

>[ \d x V n \ 2 + F(\V n \ 2 ) dx - f \ dxV \ 2 + F(\V\ 2 ) dx + f (A 2 n -r 2 )u n dx - e. 

J\x\<X J\x\<X K J\x\>X 

Passing to the liminf and using the weak convergence in [—A, +A], we infer 

lim E(V n ) - E(V) > £ lim P(V n ) - f (A 2 n - rl)u n dx 



— 6. 
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However, (36) implies 



/ (A 2 n - rl)u n dx^ {A 2 - rl)u dx, 

J\x\<X J\x\<X 



\x\<X 



so that 



E % > 



K 



p _ 

1 DO 



\x\<X 



(A 2 r 2 )u dx\-e> ^\ Poo _ P{V) \ - (l + |)e = ||P„| - (l + |)< 



Letting e — > 0, the conclusion follows. 

Step 6. There exists K > such that, if R > 0, then 



We recall the expansion P(U S ) = r 2 i\ + sP + o(s) as s — » 0, where P = . From the Hamilton 

ds |s=o 



group relation 



dE(C/,) dP(J7 a ) 



ds 

hand, by definition of c„ , 



, we also infer by integration E(U S ) = E(Uq) + ^Pq + o(s 2 ). On the other 



2Mr, 



■ Mri 



1 — cos 



, P(V n )-rfr s 



= Mr* 



1- A 1 



M 2 



for n large. Here, we have used that Mc n = sin((r 2 7r — P(V^))/r 2 ) — > Mc e [0,-ff/i*] (c/. Step 2), thus 
cos((r 2 7r — P{V n ))/r1) > 0, for otherwise, we would have, by Proposition 6.1 



Jfr(V n ) =E(V n ) + 2M4sm 2 ( 



2r 2 



) > S(J7 ) — K[i + Mr^ 



1 + 1 1 



M 2 



> E(U ) - + 2Mr A + 0{£), 



but this contradicts (32) if /i* is sufficiently small. 

We assume R > 0, so that, by Step 4, z = and c — a^. We recall the expansion ct m 
By definition of E$ , one has 

'P(V n )-r 2 n 



^■2 



f; 8 +£;(V)+Mr-g 



1-1/1- 



M 2 



< 



lim P(K)+ lim 2Mr^sin 2 (- 

n^+oo rw+oo V 



2r 2 



lim JXf(V n ) = Jf m in(n) 

n—t+oo 



since (V^) is minimizing for J(f m \ n {y). Moreover, from the expression of V, it holds (for R > 0) 



E(V)=E(U a )+2R 



°l{r 2 -^ 2 ) 2 



+ F(n 2 ) and P(V) = P(t/ CT( J + Ra^ 



M 2 ' 



Furthermore, Pj = Poo - -POO and c = Mr 2 sin((r 2 7r-P 0O )/r 2 ) with P^ e [0,r 2 7r] and cos((r 2 7r-P 00 )/r 2 ) > 
0, thus 



P S = Pco - P(V) = r 2 ^ - r 2 arcsin ( 



Mr 2 



P(U a J - Ra,- 
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Combining this with the expansion of E(U<j) and P{U a ) gives 

2 

> S(CT ) + Pj + ^P + o(<7 2 ) + M 









+ 2R 





P + 



1 i 

M 



+ o(p 2 )+4RF(0) 



> E(U ) + — r Q arcsin(cr A1 



Po + 



M 



o{n 2 )+4RF(0) 



>E(U ) + - 
>E(U ) + ^ 



M 



+ O^Pq + Ro 



(rg-M 2 ) 



2\2 



(^o 2 -M 2 ) 2 
M 2 



Pn 



o(// 2 ) + 4PP(0) 







4 


fa + si 





The right-hand side is a continuous piecewise affine function of R (the "o" does not depend on R). Since 



2^2 



hand side is a function of i? which is decreasing in [0, -Rq(m)] an d increasing in [i? (/i), +oo), with 



- > 4P(0) and P + < (since M > -P 1 by hypothesis), it follows that the right- 



4/i 2 F(0) 

Therefore, using once again that a 2 ~ j . 



^min(M) >-E(^0) + Y 



P + 



:S(C/ ) + 



P + 



M 



2 M 2 F(0) 



4i? ( M )F(0) + o(^ 2 



M 



Pi 



1 

M 



4^ 2 P(0) 



+ o(^ 2 



= p(c/ )-m ; 



2P(0) 



M 



+ ( M 2 ). 



In view of our hypothesis P + -h < 0, we infer that 



^min(M) > P(^o) + 



K 



for sufficiently small and some positive constant K, as wished. If the assumption Po + jj < is not 
verified, but if Po + > for instance, then the function of R above is increasing in [0, +oo), with minimum 
value achieved at R — and equal to 



Pn 



M 



o(v 2 ) = E(U ) + 



2/z 2 P(0) 



1 i 

M 



+ o(n 2 )>E(U ) + ^. 



We then would have concluded a stronger estimate, which is actually in contradiction with (32), hence we 
are necessarily in the case R > 0. The assumption P + jj < is however crucial for the last step. 

Step 7. We assume Po + jj < 0. Then, for ^ sufficiently small, the case R = does not occur. 

We argue in a similar way, but since R = 0, the expressions for E(V) and P(V) are given by 

^ ,„ and T>c\r\ _ vnr ^ oT c ( r o - ^?) 2 



E(V) = E{U C ) - 4 / F(|£/ c | 2 ) cfe 



P{V) = P{U c )-2 A 2 C ' <!■>■■ 
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Here, we have used that |9a;?7 c | 2 = F(\U C \ 2 ) since U c solves (TW C ). Combining this here again with the 
expansion of E(U C ) and P(U C ) gives, using that < c < Kfj,, 



^„i„(M) > E(U ) + E t + -P + o(c 2 ) + M 









1 - 




- 4 / 






Jo 



4 F(\U c n dx 



>E(U ) + ^ + C - 
> E(Uq) + arcsin(c 



2 r . ] 

^0 + 



M 



- 4zF(0) + o(p 2 



r z (r 2 - A 2 ) 2 
/M) + cP -c 1 c) dx + o(c) 

JO A c 



+ 



Mi 



Following the lines of the proof of Lemma 4, we have 

dx = 2arctan (y ~ 



Z {r 2 -Al) 2 



Al 



o( P 2 ). 



+ o(/i 2 ) -4zF(0). 



(38) 



Indeed, noticing that A c = in [0, z] with z < Kfj,, we write, expanding the square, 

r* {r 2 - A 2 ) 2 r> ^ r> ^ 



Then, using the change of variable £ = T) c (x), 



I 



*(r 2 -A 2 ) 2 



Al 



dx = 



L 
I, 

+ 



2 2 



2 2 



2 2 , 

r - r ° rj ( 

L (r 2 +£)\,r 



- arctan 
c 



(39) 



by computations similar to those for the proof of Lemma 4. This proves (38). Therefore, 



Xnin{n) > E(U ) + 









c 


[ P ° + m\ 


— 2 arctan (\l 



+ ■ 



Po + 



M 



+ o(p 2 )-AzF(0). (40) 



By (32), the left-hand side is < E(U ) + Kfi 2 . Since P + jj < 0, c < Kfj,, z < Kfj, and F(0) > 0, this 
implies 

1 



Po + 



M 



2 arctan ; 



<Kfi, 



thus 



arctan 



r 2 o + 



L) < Kfi 



and finally, for fj,* small enough, 



< 



// 



1 < Kfj 2 . 
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Combining this with the equality r 2 + £ c = — ° + C(c 4 ) seen during the proof of Lemma 4, we infer 

4_r (0) 

In particular, going back to (39) and since for < x < z, 

rl+Z c = Al{Q)<Al{x)<Al{z)=v\ 

this implies 



\<[%<\ -tan WjjL - l) + O M < ^ + *M < 



which provides (since c ~ yu) 

Inserting this into (40) and keeping in mind that the left-hand side is < E(Uq) + K^j, 2 , we deduce 



2 arctan ( ^ j 1 ) = o{jj). 



However, since arctan ( \j 2 | 1 ) > 0, this gives 



o{p) < c 



2 M ^F(0) 



yielding a contradiction for small \i since we have Pq + jj < by assumption. Therefore, the case R = 
does not occur for sufficiently small /x*. If we had Po + > 0, we would not have been able to show that 
J£nin(/-0 gi ves a control on /j,. 

The proof of Proposition 6.2 is complete. □ 

Proof of Theorem 9. Let U G f^. If ^ = inf K |£/| > 0, then Proposition 6.2 gives JT(C/) > 
P(£/ ) + fJ?/K > E(U ) = -T(U ). If inf R |[/| = 0, we deduce from Proposition 6.1 that x\v) > 
E(U ) + 2Mr 4 sin 2 ((<£([/) - 7rr 2 )r 2 ). Hence JT(C/) > E(U ) except if JT(t/) = E(U Q ). From the study 
of the equality case in Proposition 6.1, it follows that U G {e ie Uo(- — y),y £l,9 G ffi}, as claimed. 

6.3 Proof of Theorem 10 

As a first step, we shall need a quantified version of Proposition 6.1. 

Proposition 6.3 There exists e > and K > 0, depending only on f such that, for any U G Z verifying 

X(XJ) - E(U ) < e and inf \U\ < e , 

IR 

there holds 



inf d z {U, e w U (- - y)) < k(x(U) - E(U ) + inf \U\) 



1/4 



2/ 6 

e 
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Proof. First, we translate the problem in space so that fx = infjR \ U\ — |J7|(0) and shall choose the phase 
factor later. We follow the lines of the proof of Proposition 6.1 and actually get (writing U = Ae 1 ^ locally 
in {\U\ > 0}) 

/ \d x U\ 2 + F(\U\ 2 ) dx = / l lul>0 A 2 (d x <j>) 2 dx+ \d x \U\\ 2 + F(\U\ 2 ) dx 

Jo Jo Jo 



+00 



l\u\>oA (d x <f>) dx 

JO 

+00 



+00 



v*)-IWI 



dx 



+ : 



> 



2 r r ° 
dx + 2 



pro 

/ VfW) 



Jo 



s 2 ) ds. (41) 



r-roo 

2 / \VFWF)d x \U\\ dx 
Jo 

r+co r+ca 

/ i\ m>0 A 2 (d x <t>) 2 dx+ / v^£f)-IW 

Jo Jo 1 

Arguing similarly in (— 00, 0), we get 

E(U)>E(U )+ f l m>0 A 2 (d x cf>) 2 dx+ f \^F(\U\z)-\d x \U\\} 2 dx - 4 
Jr Jr l j 

The gradient of the phase is controlled using (41). We shall now estimate the modulus part. Let us denote 
A= |U"|and 

h = d x A - VfW), 

for which we have, by (41), 

Nk M) < E{U) - E(U ) + 4 f V^X^) ds < E(U) - E(U ) + A>. (42) 

Jo 

Recall that U verifies {d x U ) 2 = F(U$) in R, hence S x C/ = \/F(U§) in R + . Setting 9 = A - \U \, we infer 

d x Q = yffjA^) - ^FiU 2 ) + h in R+. 



We set, for x > 0, 



Since £/ verifies ^f/ + U (x)f(U§(x)) = and d x U = y/F(U$(x)) in R+, it follows that 

G{x,9) - y/F((Uo(x) + ey) ^F{U 2 {x)) - f^o- 

Moreover, by Taylor expansion, we infer the existence of K > and 60 > such that, for \6\ < 0o, x e R + , 

|G(a;,0)| < X6» 2 . 

The estimate is clearly uniform in view of the exponential decay of d x Uo at infinity. Therefore, 



o x e 



d 2 x Uo(x) 
d x U (x) 



e + G(x,e) + h(x). 



(43) 



We view this ode as a linear ode with source term G(x,Q(x)) + h(x). Since d x Uo solves the homogeneous 
equation, we infer from Duhamel's formula and the fact that 0(0) = A(0) — Uo(0) = \U(0)\ = \i that for 
x > 0, 



Q(x) = fi + d x U {x) / dz. 

Jo OxUo(z) 



(44) 
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We shall prove that this equation implies that if /x and |/i|l 2 (b + ) are sufficiently small, then 



|e| L2(R+) < K(ih\\ L 2 {R+) 



/i 



(45) 



We assume fi < 9q/2. Note that since Uq is a kink, we have the decays given in Proposition 1. Hence, there 
exists two positive constants K\ and K 2 such that 



Vx e R+, 



< d x U (x) < K 2 e~ 



In particular, if |©(a;)| < 9 a in the interval [0,R], then (44) implies, for x G [0,R], 



\&(x)\ < A* + K.K^-^ f k||e|| L o O([0 . ii]) |e(z)| + |/i|(z) 
Jo 

K X K 2 



dz 



[1 



KK\K 2 2 



II e 1 







by Cauchy-Schwarz. We thus choose ||/i||l 2 (r + ) + A* sufficiently small so that 

K 1 K 2l 



4(, + -^|% 2(K) )<M m in{^^}. 



Then, we consider the set "ft of all R > such that \Q(x)\ < 9 in the interval [0,R]. Since 6 e ff^R.C) 
is continuous by Sobolev embedding and |0(O)| = [i < 9 0} 1Z 7^ and is closed in R?j_. Moreover, the above 
estimate shows that for R GlZ, 



,, , KK x K 2l% „, 

«e|Uoc ([0 , fl]) < M + / lie" 2 



which gives 



and then 



II en 



l°°([o,ji])^1 ||oIl~([o,r])J < + -^HNIl 2 (r + ) 



£°°([0,fl]) + fey- I"Ii 2 (R+)l 



||@IU<*>([0,.R]) < 2 



, M2 . 



<i<% 



(46) 



Consequently, 1Z is open in By connexity, 1Z = proving (45). In what follows, we assume ||/i||l 2 (m + ) + 
\i sufficiently small so that ||©||i,°o < 9 , thus \G(x, 9)| < KQ 2 . In particular, 



|e(x)| < a* + k,k 2 f e - e '<*-*> f^||e|u=o ([0 , fl]) |e(^| + |/i|^)" 

Jo L J 



dz. 



For i? > to be determined later, we then deduce from classical convolution estimates that 

||0|U 2 ([o,ii]) < M\ / ^ + ^3||e|| L oo (K+) ||e|| L 2 ([0;fl]) +K 3 \\h\\ L 2 { 
Imposing |/i|i,2(r + ) + \i smaller if necessary, we may assume that 

1 



so that we get 



^3||e|| L oc (K+) < K 3 K(\h\ L ^ +) + M ) < 2 

II©|l 2 ([o,_r]) < K i ^v / R+ \\h\\ L 2 {m+) y 
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Reporting this into (43) provides 



ll^ef 



L 2 ([0,fl]) 



< K, 



(ji 2 R + 



li 2 (R + ) 



Arguing similarly in [— R, 0] and using (42), we obtain an H 1 estimate for G in [—R,R]: 



lie I 



HH[-R,R]) 



<K e (E(U)-E(Uo) + n 2 R + v)- 



(47) 



We now turn to the estimate in {\x\ > R}. For that purpose, we write 



\x\>R 



< :! r\U\f + ^(\U\ 2 - rlf dx < E(U) - E(U ) + 



\x\>R 



(d x U f + F(U 2 ) dx 



-f (d x \U\\) 2 + F(\U\ 2 ) dx + f (d x U ) 2 + F(U 2 ) dx. (48) 

J\x\<R J\x\<R 

Since Uq decays exponentially (see Proposition 1), it follows that 

\d x Uo\ 2 + F(U 2 ) dx<Ke~ c ° R . 



L 



\x\>R 



Furthermore, by integration by parts, 

- f (8 X \U\) 2 + F(\U\ 2 ) dx+ f (d x U ) 2 + F(U 2 ) dx 

J\x\<R J\x\<R 

= - I 2d x U d x Q - 2U f(U 2 )Q dx- f (d x e) 2 + F([U + 6] 2 ) - F(U 2 ) - 2U F'{U 2 )Q dx 

J\x\<R J\x\<R 

< f 2Q[d 2 x U + U f(U 2 )} dx - 2Q(+R)d x U (+R) + 20(- R)d x U (-R) + K\\e\\ 2 H1([ _ H +R]) 

J\x\<R 

< Ke~ c ° R + k(e(U) - E(U ) + p?R + m) . 



For the before last line, we have used that 9 ^ F([U + 0] 2 ) - F([/ 2 ) - 2U Q F' {U%)6 is O{0 2 ) as 9 and 
for the last line, that Uq solves 2 Uq + Uo/(Uq) = 0, the exponential decay of d x lIo and the uniform bound 
on 9. Reporting these estimates into (48) provides 



l lffi({|*|>.R}) 



|x|>fl 



(^|(7|-^|C/ |) 2 + (|[/|-|(7 |) 2 dx 



< 



2 f (d x \U\f + (d x \U \) 2 + (\U\ - r f + (\U \ - r ) 2 dx 

J\x\>R. 



<K 



!\x\ 

E(U) - E(U ) + e~ CsR + (i 2 R + p 



Combining this with (47), we deduce that for any R > 0, there holds 



liei 2 ^ 



< K 



E(U) - E(U ) + e'^ H + n 2 R + /i 



We then choose R — ^i 1 if > or i? ^ +oo if \i = 0, and get 

|eUi(R) < K^E(U)-E(U )+fi. 

Notice that if /' < everywhere, then we may give a quick proof of the above estimate, since, using here 
again integration by parts and that d x lIo + Uof{Uo) = 0, we may deduce that 

r+oo p 

E{U) - E(U ) > -4nd x U (0) + / (d x S) 2 dx + / F((U + 6) 2 ) - F(U 2 ) - 2U QF'(U 2 ) dx, 

Jo Jr 
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and since /' < 0, F((U + 6) 2 ) - F(U§) - 2U 6F'(U§) > 2 /K by Taylor expansion, providing the desired 
H 1 bound on 0. 
Observe now that 

Jf(U) - E(U ) > E(U) - E(U Q ) > [ l\u\>o[Adx<f>? dx, 

JR 

hence 



\\0 X U - d x U \\ L 2 m = \d x (\U \ + O)e i H m>0 + il lu>0 Ad x( f>e^ - d x U \\ L 2 (m 

< l|e^l|!7|>o^|^o| - d x U a \\ L 2 (R) + \\l\u\ >a Ad x (j)\\ L 2 {M) + ||8| L 2 (R) 

r 1V2 

< \e+l\ u>0 d x \Uo\ - d x U \\ L 2 (R) + K JfHJJ) - E(U ) + ii . 



(49) 



We distinguish now the cases [i = and \i > 0, and begin with the assumption \i > 0. Then, we have a 
global lifting U = Ae 1 * and 

dz(U, U ) = \\d x U - d x U \\ L 2 m + \\U\ - \U \\\ L 2 m + \U(0) - U o (0)\ 
= \\d x U - d x U \\ L 2(R) + |6|l 2 (r) + A 4 

r i V 2 

< \e*d x \U Q \ - d x U \\ L 2 m + K X{U) - E(U Q ) + n . 



Now, we notice that 



l^dM-dMl^ = 2 f [(d x U ) 2 - d x U d x \U \ cos(0) dx 

= 2^ + (^(7o) 2 (l-cos(0)) dx + 2j (a,C/o) 2 (l + cos((/))) dx (50) 



and that 



J(f{U)-E{U ) > 2M^sin 2 > ±-(P(U) - r 2 n mod2nr 2 ) 2 . 



2r 2 



We define (5 = pT(t/) - £(t/ ) + pt) 1/4 - By Cauchy-Schwarz, we have 

K 



f (A 2 rl)d x(l 

J\x\>5 



dx 



Inserting this into (51) gives 



< 



inf|x|><5^1 v J\x\>s 



([ (A 2 -r 2 ) 2 dx) 1/2 ( f (Ad x cfifdx) 

K J\x\>8 ' K J\x\>S ' 



1/2 



< 



K 



inf|a;|><5 A 



(E{U) - E(U ) + M ) 



1/2 



(51) 



\x\<5 



(A 2 - rl)d x 4> dx-rliiYCLO<12irrl < K \{JfT{U) - E(U )) 1/2 + t 



in£\x\>s A 



(e(U) - E{U ) + fi 



1/2- 



< 



K 



mi\ x \> s A 



(jfT{U) - E(U ) + 



1/2 



In addition, by Cauchy-Schwarz, 

A 2 d x (j> dx 

\x\<8 



< 



y/2S( sup A) (jf(U) - E{U ) + ft) 



1/2 



' \x\<5 
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Consequently, 



$(+8) - <j>(-8) - 7T mod 2tt < / (A 2 - rl)d x (j) dx - r^n mod 2nrl 

J\x\<8 

+ V2S ( sup A) (JT{U) - E{U ) + n) 1 



• \x\<S 



< 



K 



inf \ x \>s A 



+ V^( sup A)] (j{f(U) - E(U )+ v) 1/2 ■ (52) 



' kl<<5 



From our choice 8 = (JXT(U) - E(U ) + < 1 and since ||e|| L ~ (R) < K{X \U) - E(U ) + ^) 1/2 - 0(8 2 ), 
we infer inf| x |>^ A > inf| x |>5 \Uq\ — ||0|l°°(r) > 8/K. Similarly, sup| x |< 6 ^4 < sup| x |< 5 \Uq\ + || ©|jx,=~ (m) 

< K8. 

Reporting this into (52) yields 



<t>(+8) - 4>(-8) - tt mod2-7r 



< Kb. 



We now freeze the gauge invariance by imposing (j>(+8) = 0. Note that then <fi(—S) = n + 0(S). Furthermore, 
since <fi(+S) = 0, 



J\x\>6 



dx < 



j. f K A)2 f A 2 {d x <f>) 2 dx < §6* = KS 2 , 

(mf|a;|>5 A) 2 J\ x \>5 S 2 



< KSy^ 



which implies, for x > 5, 

\1 - cos{(f>{x))\ < |l-cos(0(O))|+ J d x <j>sm((j)) 

and similarly, since cos((f>(—8)) = cos(7r + 0{8)) = — 1 + 0{8 2 ) 1 for x < —5, 

|1 + cos(0(x))| < K6y/\x\. 

We turn back to (50) and infer 

\\e^d x \U \ -d x U \\ 2 L2(m <KS + 2^ {d x U ) 2 (l - cos(^)) dx + 2 J {d x U Q ) 2 (l + cos($) dx 

<KS + KS [ (d x U ) 2 = if (5. 

Inserting these estimates in (49), it follows that 

d z (U,U ) <KS. 

We now turn to the case /i = 0. Without loss of generality, we may assume that \U\ > in (— oo,0) 
(since \U\ -» r > at ±oo), and let £ > be such that \U\(£) = and |{/| > in (£, +oo). We first estimate 
^ by writing that 

1/2 



|c/ |(£) = |t/| w + e(£) = Q(£) < ||e|| L oc (K) < x (V(c/) - £(c/ ) + m) = ks 2 , 



thus 



£ < KS 2 . 



Moreover, we have two local liftings U = Ae l ^+ in [£, +oo) and U = Ae 1 ^- in (— oo,0). Going back to (49), 
we then deduce 



dz(U,U ) < \\e^-d x \U \ - d x Uolm-oo,o) + te^+d x \U \-d x Uo\\me,+oc) + KS + K Jf{U)-E{U ) 



1/2 
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Arguing as for the case fi > 0, we obtain \U\ = A > S/K in [£ + S, +00) and in (—00, —5). By definition of 
*P, we have 



(iu\d x u) + 



+00 i' — b 

{A 2 - rl)d x <t>+ dx + r 2 M£ + 6)+ (A 2 - r 2 )0 x ^ dx - r 2 ^ + {-5) 

J —00 



£+S 



in R/(27rr 2 Z), hence the same arguments as in the case \i > provide 



(/)+(£ + 5) - <p+(-6) - it mod27r 



< K5, 



since the integral J_ g (iU\d x U) is bounded by K\fb by Cauchy-Schwarz. Imposing 4>+{l + 5) for the gauge 

invariance, we infer 1 — cos{4> + {(, + 5)) = and <f>+(— 5) = tt + 0(VS) mod27r, hence l + cos(</>_(— 5)) = 0{5). 
Therefore, we conclude as before that 

d z {U,U Q ) <K5, 

which finishes the proof of the Proposition. □ 
In order to prove Theorem 10, we use Proposition 6.2, which provides 

JT{U)>E(U ) + ± (infill) 2 , 

thus 

H = inf \U\ < Ky/j(r(U) - E(U ). 
Inserting this bound in Proposition 6.3 then gives 

d z (U, U Q ) < K[jf(U) - E(U ) + K^Jtr(U)-E(Uo)] ^ < K t]X{U) - E(U ), 
and the proof is complete. 



7 About the stability analysis for the sonic waves (c = c s ) 

We have left aside in our study the case of the sonic waves (c = c s ), but would like to say a few words on 
the difficulties associated with this critical case. 

We note that if there exists a sonic nontrivial travelling wave, it does not vanish, hence we may use the 
hydrodynamical formulation (15) of (NLS) as in [47]. The point is that the Sturm-Liouvillc operator (see 
section 4 in [47]) 

L = -TAWl^)d~x) +q{x) > 

with 

a(x) - {dxV)2 9 ( 9xV ) 1 f (r 2 71) c2r ° 
91 J " 4(rg - 77)3 5iB U(rg-rj)V 2 /l V > 4(r 2 - 

has, by Weyl's theorem, essential spectrum a ess (L) = [0, +00) when c = c s . Indeed, we know from Proposition 
1 that rj Cs and its derivatives tend to zero at infinity, hence, as x — >• ±00, q(x) — > — \f'{r 2 ) — = since 
c 2 = c 2 = — 2rQ,/'(r 2 ). Therefore, there does not exist 8 > such that (Hp,p) > S\\p\\ 2 for any p orthogonal to 
the subspace spanned by the negative and the zero eigenvalue, and thus the Grillakis-Shatah-Strauss theory 
does not apply. 

dP 

In the case — — < 0, where it is natural to expect stability, a natural thing would be to try to work 

dc |c=c s 

with the functional 

jSfty) = EM) - c s P(tP) + f (Pty) - P(U Cs )) 2 
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and to follow the lines of the proof of Theorem 9. Indeed, the spectral analysis shall not give positive 
definiteness of the Hessian due to presence of essential spectrum down to 0. Therefore, we may study at 

dP 

fixed /i = infjR \tp\ close to infg \U C I. When < < c s and — — ^ 0, the infimum of \U C \ contains a 

dc \c=c, 

neighborhood of infjg \ U C , | for c close to c*. For = c s , this is no longer the case: we have only a one sided 
neighborhood of iniR \ U Cg \. It is plausible that the study for fi in this one sided neighborhood of infm \ U Cs \ 
can be done as in the proof of Theorem 9, but for the remaining values of fi, we have to find a sharp ansatz, 
which is not very easy to find. 

dP 

Furthermore, for the linear instability which is expected if — — > 0, let us mention the following 

dc \c=c s 

point. For the eigenvalue problem studied in [6], the characteristic equation for the constant coefficient limit 
at infinity, namely 

r 4 -(c 2 -c 2 )r 2 -2 Cst Ar + A 2 =0. 

becomes, when = c s , 

r 4 -2c s Ar + A 2 = 0. (53) 

The behaviour of the roots for small A is then different from the case < < c s . Indeed, there exists a 
root ~ A/(2c s ) for A — »■ 0, and for the three other roots, we use the variable r = \f\~z, which transforms 
r 4 - 2c s Ar + A 2 = into z A — 2c s z + A 2 / 3 = 0. This last equation has, for A -> 0, three roots ~ j k ^/2c~ s , 
where j = e 2t7r / 3 and k = 0, 1, 2. In particular, (53) has three roots ~ j k y / 2c s X, k = 0, 1, 2. The value 
A = is then a branching point, and we shall have a smooth problem not in A but in s/A. Since analyticity 
is not necessary for our purpose, we may define an Evans function D in R + , smooth, and such that, for 
A > 0, fl^s/A) = if and only if A is an unstable eigenvalue for (27). another difficulty comes from the fact 
that it will be difficult to find an analytic extension of the Evans function D near since, by Proposition 1, 
for = c s , m* and 77* decay only at an algebraic rate and not an exponential rate. Consequently, we can 
not use the Gap Lemma of [30] and [40]. Finally, as a straightforward computation shows, the stable and 
unstable subspaces for the eigenvalue problem are transverse for A > but their continuous extensions at 
A = have a nontrivial intersection. Therefore, both stability and instability requires some further analysis, 
and the situation is then much more delicate than the one studied in subsection 5.1. 



Appendix A: construction of a Liapounov functional in the stable 
case in the Grillakis-Shatah-Strauss framework 

We work with the notations of [35], and recall them briefly. We consider a Hamiltonian equation in a real 
Hilbert space X, with scalar product (•, -)x, under the form 

% = JE'(u), (H) 

where J : X* — > X is a closed linear operator with dense domain and skew-symmetric. Assume that T 
is a Co-group of unitary operators in X generated by T'(0), which is skew-adjoint and with dense domain, 
and that E is invariant by T, that is E(T(s)u) = E(u) for any s e K, u e X. Assume moreover that 
T(s)J = JT(-s)* for any set and that there exists B : X -> X* linear and bounded such that B* = B 
and JB is an extension of T"(0). We then set 

Q( u ) = ^{Bu,u)x%x- 
The basic assumptions of [35] are the following ones. 

Assumption 1 (existence of solutions) : For any r > there exists t* > 0, depending only on r, such that 
for any u m e X, there exists a u e C((— X) with u(0) = u m solution of (H) in the sense that for any 
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tpe D(J) c x* , 

^(u(t),tp) x *,x = -(E'(u(t)),J(p)x',x in V((-t„ t,)), 

and verifying E(u(t)) = E(u in ) and Q(u(tj) = Q(u in ) for t G (-<*,*»). 

Assumption 2 (existence of "bound states"): There exists an interval SlcK, not reduced to a singleton, and 
a mapping !] 9 ui H ^ £ of class C 1 such that, for any uj e f2, 

= wQ'(&,)> iefl(T'(0) 3 )nD(J/r'(0) 2 ), T'(o)<^o. 

Assumption 3 (spectral decomposition): For each w e 0, the operator = E"(<j> u ) — toQ"(<p u ) : X — > <Y* 
has its kernel spanned by T'(0)(/> u , has one negative simple eigenvalue and the rest of its spectrum is positive 
and bounded away from zero. 

Under assumption 2, we consider some € and the associated bound state and then define, for 
M > 0, the functional 

Mi \ 2 
X. («) = E(u) - u*Q{u) + — [Q{u) - Q(0„. ) J . 

It is clear that </> w „ is a critical point of JS^,,: JSf' (</> w ») = E'(<f> Utt ) — uj*Q'(4>u,) = 0. We denote by 

A = S%^<f> u .) = H U , + M<Q , (0 u .),->*.,*Q , (0 ta ,.) 

its second derivative, which is a self-adjoint operator. The main result of this appendix is the following. 

Theorem A We make assumptions 2 and 3 and suppose moreover that the operator (Q l (<j> u ,), •) x* ,xQ' \<i>u,) 
is a compact perturbation of H^, ■ If ^^"^ < and 

dU! |ui=u>» 

M > 



dQ(4>gj) 

did \u>=ui, 

there exists S > such that 

VveX, s.t. (v,T'(0)cf> u J x =0, (Av,v) > S\\v\\ 2 . 

In particular, for any u g X with inf sS R \u — T{s)(j)^ r || 2 < e, we have 

inf |u - T(s)4> u , f < - (jSf («) - 2 (<f> u , )) . 

Therefore, when assumption 1 is moreover verified, the (global) solution u(t) to (H) with initial datum u ln 
verifies 

sup inf \u{t) - T{s)<t> u ,f < \ (^{u m ) - &(4>u,j) < K\\u in - <^J 2 , 
tewt seR d \ / 

provided the right-hand side is sufficiently small. 

We point out that the condition that the operator (Q'(<j>w»), •) x- ,xQ '(<t>u,) 1S a compact perturbation of 
H Ut is not very restrictive, since in many cases coming from PDE's, it involves less derivatives than ff Ut 
and Q'{(j>^,) tends to zero at spatial infinity. 

This type of Liapounov functional has been used by I. Barashenkov in [3] to prove that the travelling 

dP 

waves of (NLS) in dimension one are stable when — — < 0. The proof follows basically the one in [3], but 

dc 
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some points have to be clarified. The interest of this type of Liapounov functional is that the saddle point 
4>u, is now a non degenerate local minimum for ■ This is a great advantage for numerical simulation 
of the "bound states", since a gradient flow method on Jzf^, can be used. This approach has been used, 
with a very similar functional, in [56] by N. Papanicolaou and P. Spathis for the numerical simulation of 
the travelling waves for a planar ferromagnets model. In the same spirit, in [17], we also use a gradient 
flow method on this type of functional for the numerical simulation of the travelling waves for (NLS) in two 
dimensions. 

Proof of Theorem A. Recall that the spectrum of H Utt is, by assumption 3, such that — A^ e a{Hu,), 
G a(H Ut ) and cr(iJ w J \ {— A^, 0} C [S, +oo) for some 5 > 0. Since we assume that (Q'{<t> u ,), ')Q'{<t>ui,) is a 
compact perturbation of H^, , the essential spectrum of A is the same as the one of H u , , hence is included 
in [5, +oo). Furthermore, € cr(H u ,) and ker(i? tlJi( ) = RT"(O)0 U „ by assumption 3. Since Q'{4>w,) = 
B^u, and JB is an extension of T'(0), we have (Q' '(</>„„), T' '(O)^,) x* ,x = (B<j> Ut , JB<j> u ,)x*,x = 0, hence 
A(T"(0)(/> w ,) = 0. Noticing that (Q'(<f>w,), •) x* ,xQ i s a nonnegative operator, we infer that ker(A) = 
ker(H LOt ) = KT'(O)^, is one-dimensional. Therefore, it suffices to show that A has no eigenvalues in (-co, 0). 
As we have seen that (Q'(<fiui,), - )x* ,xQ' (<Pu>,) is a nonnegative operator, we deduce that u(A) C [—A*, +oo). 
Let us first show that — A^ ^ f(A) by contradiction. If — A^ is an eigenvalue of A, then there exists »el, 
ti^O, such that = (A + Xl)v = (H + Xl)v + M{Q'((f> Ut ), v) x *,xQ'{<t>u t )- Taking the duality product with v 
yields = + X%)v, v)x*mx + M(Q'(<j) Ur ), v) 2 x , x . Since the two terms in the sum are nonnegative, this 

implies {Q' v)x»,x = and {{H Wli + X%)v, v)x*,x — 0, which in turn implies v € ker(H UJt + A*) = R\ 
(here, x is a negative eigenvector of H Ut for the eigenvalue — A^ < 0. As a consequence, we must have 
(Q'{<f>u,),x)x'.x = 0. On the other hand, differentiating the equality £7'(<^ w ) — uQ'(<t> u ) = at w = w, 
yields Q'(0 w J = H u ,<(>', where <f>' = ^ {u=Um - Thus we must have = (H u ,<j>',x)x*,x = (H u ,x, <j>')x*,x = 
— A^(x, 4>')- Therefore, <f)' is orthogonal to x an d this gives (H^^tp' , <fi')x* ,x > 0. However, this is not possible 

dQ(4>u) . dQ{<t> u ) i jj ,/ ii\ \ . f , ■ .. 

it — < U, since we nave — = —{H^^tp , q> )x*,x- As a consequence, it A is a negative 

dui \uj= oj* diO 

element of the spectrum of A, then — A^ < A < and A is an eigenvalue: there exists v e X such that v ^ 
and 

Xv = Av = H^v + M(Q'(^J, v)x*,xQ'(^,)- 
Since — A^ < A < 0, we then infer 

v = -M(Q'(<j) Ut ,),v)x*,x(Hu, - X.y 1 Q / (<f> ul „). (A.l) 

Since v ^ 0, we can not have (<3'(<^ w ,), v)x*,x — 0. Then, taking the scalar product of (A.l) with l~ 1 Q'(<fru.) 
(here, I : X — > X* is the usual Riesz isomorphism) gives 

ff(A)=0, where g(t) = 1 + M((H U , - t)- 1 Q'(<f> u ,),I- 1 Q'(<t> u ,)) x , -\l < t < 0. 

It is clear that g is smooth in (— A*, 0) and that 

g'(t) = M((JT U . - irVO^J^Wu,))* - - ^WoOl* > 0. 

We now study the limit of g at 0~. Let us recall that H u ^<f>' = Q'(</> w „) and that we have already seen that 
{Q'(<i>u.),T'(0)<l> u .)x*,x = 0, i.e. l^Q'^J is orthogonal to ker(i? w J. Therefore, as t -> 0", 

- trWuj.i-Wu,.)) -+ (^ / ,i- 1 Q / (^j) = {Q'{<t>M') X '. X - 

and it thus comes 

g(jt) ^ 1 + M MM as t _>o- 

Since ^^^"^ < by hypothesis, it follows that if M > — 1 / ( - 1 . _ ) > 0, the function g increases 

in (— Aj,0) and tends to some negative limit at 0~. In particular, g is negative, hence we can not have 
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g(X) = with A £ (— A^,0). We have therefore shown that the spectrum of A consists in a simple eigenvalue 
with eigenspace spanned by T"(0)</> u , and the rest of the spectrum is positive and bounded away from 0. 
This concludes the proof. □ 



We would like to point out that in the proof of [3], the fact that —\ 2 ^ cr(A) was not shown, the kernel 
of A was not studied and the essential spectrum was not considered. Moreover, the functional spaces are not 
given, hence we do not know for which perturbations stability holds. 

Appendix B: from linear to nonlinear instability 

We still consider in this appendix an abstract Hamiltonian equation in the framework of [35] 

^ = JE'(u) (H) 

on the real Hilbert space X, with scalar product (-,-)x- Here E : X — > R is of class C 2 and J : X* — > X is a 
closed linear operator with dense domain and skew-symmetric in the sense that (u, Jw)x = — (w,Ju)x»,x 
for u e X, w e X*. 

We assume that there exists a Co-group T of unitary operators in X generated by T'(0), which is skew- 
adjoint and with dense domain, and that E is invariant by T, that is E(T(co)u) = E{u) for any oj e R, 
u E X. Assume moreover that T(oj)J = JT(-ui)* for any uj e K and that there exists B : X — > X* linear 
and bounded such that B* = B and JB is an extension of T'(0). We then set 

Q( u ) = \{Bu,u) x *,x, 

which is invariant by the flow (H) (see [35]). By "bound state", we mean a particular solution U of (H) of 
the form U(t) = T(u!t)<p for some w€l and where <j> £ X, <fi ^ 0. In other words, E'(<p) = ajQ'(<f>). 

There exists an open interval S!cR, not reduced to a singleton, and a mapping ft 9 w ^ ^ u e X of 
class C 1 such that, for any to € ft, 

E'(M=wQ'(^), <f> u £D(T'(0) 3 )nD(JIT'(0) 2 ), T'(0)<^ ± 0. 

The solution U(t) = T(uit)<f> is said to be stable in X if for any e > 0, there exists S > such that 
any solution to {%) with initial datum u™ € Bx(4>,$) is global in time and remains in Bx(<f>,e) for t > 0. 
Otherwise, it is said unstable. This supposes some knowledge on the Cauchy problem for {%) (at least 
existence of solutions). If we are given some Banach space y D X with continuous imbedding X ^ y, 
we may also say that the solution U(t) = T(ojt)<f> is said to be stable from X to y if for any e > 0, there 
exists S > such that any solution to (H) with initial datum u m £ Bx(4>,S) remains in By ((f), e) for t > 0. 
Clearly, a solution stable in X is precisely a solution able from X to X, and is also stable from X to y, hence 
instability from X to ^ is a stronger statement that instability in X. 

In our framework, the notion of orbital stability is more relevant. Let us consider G a group and 
T:MxG^ <S J£ C (X) a unitary representation of M. x G on X, extending T : R — > X and leaving E and Q 
invariant. Then, U(t) = T(ujt)4> is said to be orbitally stable in X (for the group G) if for any e > 0, there 
exists S > such that any solution to (H) with initial datum u ln £ B((f>,5) is global in time and remains 
m U(w, ff )eRxG-E"(Tr(u;, g)<f), e) for t > 0. We may also define orbital instability from X to y D X in a natural 
way. 

In [35, 36], a general framework for the stability analysis for the "bound state" has been given. In 
particular, the nonlinear orbital instability in proved in [35] through the construction of a Liapounov type 
functional. However, this method does not give a clear understanding neither of how we get farther from 
the "bound state", nor on which timescale it occurs. 

The need for allowing an additional group of invariances G can be seen in the case of bound state solutions, 
that is U(t) = e 2 "*^ to the Nonlinear Schrodinger equation 

id t ^ + A* + */(|*| 2 ) = 0, (NLS) 
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or the Nonlinear Klein-Gordon equation in M. d 

<9 2 * = A* + tf/(|tf| 2 ), (NLKG) 

since then, the invariance by translation in space must be taken into account in the definition of orbital 
stability, and we are in a case where G = K d acts naturally by translation. The translations are taken into 
account in [15]. In [35] and [36], the notion of orbital stability is for G trivial. It is clear from the definition 
that orbital stability for G = {0} implies orbital stability for arbitrary G. For the instability in the Nonlinear 
Schrodinger equation or the Nonlinear Klein-Gordon equation, [35] and [65] work with radial H 1 functions. 
The fact that this also implies the orbital instability with the action of G = M. d by translations follows 
immediately from the fact that for any 9 € [0, 2n} the manifold DJle = {e ld (p(- — y), ye R d } is orthogonal to 

For the stability analysis of a "bound state" U(t) — T(w*i)(</> Wi J, it is natural to consider the linearization 
of {%) near <f>. More precisely, we linearize according to the ansatz u(t) = T(ui lf t)(<j) i j r + v(tj), so that the 
"bound state" becomes stationary. The linearized problem then becomes 

— = J{E"{4>) - uQ"(<f>))v = JCv, (Hlin) 

where, I : X — > X* denoting the Riesz isomorphism, J = JI : X — > X is skew-adjoint. 

The purpose of this appendix is to give a general result, for Hamiltonian equations, showing that linear 
instability implies nonlinear (orbital) instability. By linear instability, we mean that the complexification of 
[JC]c has at least one eigenvalue in the right-half space {Re > 0}. The argument follows ideas from the 
works of F. Rousset and N. Tzvetkov [61, 62]. 

Showing the existence of an unstable eigenvalue can be done through various techniques: [36] (in the 
framework of [35] when J is onto), [34] (assuming a special structure of the Hamiltonian equation); for 
uses of the Vakhitov-Kolokolov function, see [20], [21] or [58]. When J is not onto, we quote [49]. For one 
dimensional partial differential equations, one may also use Evans' function (see the survey [63]) as in [57], 
[26], [40], [70]. The paper [48] proposes another approach which allows to treat pseudo-differential equations, 
such as the BBM equation, the Benjamin-Ono equation, regularized Boussinesq equations, the Intermediate 
Long Wave equation... 

In order to pass from linear to nonlinear instability, the following result is standard. We refer to the 
paper [38] by D. Henry, J. Perez and W. Wreszinski. It can also be found in [34] and [64]. 

Theorem B.l ([38], [34], [64]) We assume that A generates a continuous semigroup on X and that a(A) 
meets the right-half space {Re > 0}. We assume moreover that F : X — > X is locally Lipschitz continuous 
and verifies, for some a > 0, ||-FXv)||x = 0(|w||^ a ) as v — »■ 0. Then, the solution <fi — is unstable for the 
equation d t v — Av + F(v). 

In [64], it is claimed that an orbital instability result can also be established. Theorem B.l shows 
nonlinear instability without assuming that the equation is Hamiltonian. However, if (Hu n ) can be solved 
using a semigroup, it does not give the growth of its norm. Moreover, it does not say that if the initial datum 
is in a most unstable direction, that is an eigendirection of A corresponding to an eigenvalue of maximal 
positive real part (plus the complex conjugate if necessary), then you can track the exponential growth of the 
solution. In particular, it does not explain the mechanism of instability and does not give any information 
on the timescale on which you see the instability. For instance, some strong instability results are shown 
by proving blow-up in finite time (see [8]), but the instability due to an exponentially growing mode holds 
on a much smaller timescale. We wish to provide here some results clarifying the instability mechanism by 
tracking the exponentially growing mode. 

B.l A spectral mapping theorem for linearized Hamiltonian equations 

When we want to prove a nonlinear instability result from a linear instability one, one needs some informa- 
tion on the growth of the semigroup JC, when such a semigroup e tJC exists, which we shall assume in this 
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appendix. The growth estimate on e tJC relies classically on the following spectral mapping result due to J. 
Priiss [59], which generalizes the work of L. Gearhart [27]. 

Theorem B.2 ([59]) Let X be a complex Hilbert space and A an unbounded operator on X which generates 
a continuous semigroup e tA on X. Fort £ (0, +oo), we have 



<j{e tA ) \ {0} = (e A *, either (x + —Z] n a (A) ± 0, or sup (A—X— —) = +oo). 

The following result is an immediate corollary. 

Corollary B.l Let X be a complex Hilbert space and A an unbounded operator on X which generates a 
continuous semigroup e tA on X. Assume that for any 7 £ R* , we have 



lim sup 

|t|->+oo 



(^A — 7 — ir^j 



< +00, 

se c (x) 



and that there exists do € [0, +00) such that a css (A) = {id, d £ R, > do}. Then, for any t £ (0, +00), 
the spectral mapping holds: a(e tA ) \ {0} = e ta ( A \ 

Proof. Since a css (A) = {id, d £ R, \d\ > d a }, we have S 1 £ e ta ^ £ a{e tA ). If A £ C docs not have 
modulus one, then note that when (A + ^Z) n cr(A) — 0, the supremum for k £ Z in Theorem B.2 can be 
+00 only when |fc| — > +00, and we conclude with our hypothesis. □ 

The fact that we exclude in the spectral mapping theorem just comes from the fact that we consider a 
semigroup and not a group. However, in most Hamiltonian PDE's, we have time reversibility and we have 
actually a continuous group and not only a semigroup. In most cases, we work with A : D{A) £ Y — > Y where 
Y is a real Hilbert space, thus for applying Theorem B.2 or Corollary B.l we have to consider, as usual, the 
complexified operator A c : D(A C ) = D(A)®iD(A) <ZY C = Y ®iY Y c defined by A c (u + iv) = Au + iAv. 

It seems that the first time Theorem B.2 is used to prove a growth estimate on a semigroup was by 
T. Kapitula and B. Sandstede in [40]. Later, the work [31] by F. Gesztesy et al. also uses this result for 
bound states for (NLS). The bounds on the resolvent in [40] were proved using the particular structure of 
the linearized operator. In [31], the computations arc more involved and rely on suitable kernel estimates of 
some Hilbert-Schmidt operators. The same type of estimates have also been used in [21]. 

The main objective of this appendix is to provide a generalization of these results to a wide class of 
Hamiltonian equations. Indeed, the approaches in [40] and [31] seem specific to the problem. In addition, 
it is not clear whether the computations in [31] and [21] can be extended to other types of equations. In 
particular, in [16] and in the present paper, we have a situation similar to the one studied in [21], namely 
travelling wave solutions to a Nonlinear Schrodinger equation with nonzero condition at infinity, but for 
nonzero propagation speeds, the travelling wave is not real-valued (as it is the case in [21] for stationary 
waves or for bound state solutions), and the bloc diagonal structure of the linearized Hamiltonian disappears. 
An additional difficulty is that in [16] and the present work, the limits of the travelling waves at +00 and 
—00 differ. 

The proof we give is based on ideas from [61, 62] and make very few spectral assumptions on C. 

Assumption (A): The spectrum of £ consists in a finite number (possibly zero) of nonpositive eigenvalues — /ii, 
... , — Hq in (—00, 0], each one with finite multiplicity, and the rest of the spectrum is positive and bounded 
away from 0. Furthermore that for any 1 < k < q, we have ker(£ + /ifc) £ D( J) and J'[ker(£ + fik)] C D{C). 
Finally, there exists i? € [0, +00) such that er css (J£) = {id, d £ R, \d\ > d }. 

The first hypothesis on the location of the spectrum of C is quite weak, since it is verified when C is 
bounded from below and has essential spectrum positive and bounded away from zero. Indeed, if 5 > is 
such that cr ess (£) £ [25, +00), then the eigenvalues of C in (—00, 5] are isolated, of finite multiplicity, and are 
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bounded from below by assumption. The second hypothesis ker(£ + /j,/.) C D(£J) is a regularity assumption 
on the eigenvectors. 

Let us recall that Theorem 11 ensures that the number of eigenvalues (with algebraic multiplicities) of 
J £ in the right-half space {Re > 0} is less than or equal to the number of negative eigenvalues of £, hence 
is finite under assumption (A). Let us now state our main result, the proof of which is given in section B.3.1. 

Theorem B.3 We make assumption (A) and suppose that J £ generates a continuous semigroup. Then, 
for any t G (0, +oo), the spectral mapping holds: a(e t ^ c ^ c ) \ {0} = e t(T ^ Jc ^ c \ Furthermore, denoting 

7o = sup {Rc(A), A e <r([JC]c) n {Re > 0}} G [0, +oo), 
for any f3 > 0, there exists M{(3) > such that, for any t>0, we have 

\\e tJ£ y c{x) <M(P)e^> + V t . 

Assume in addition 70 > and denote 

m = max ^algebraic multiplicity of A, A G cr([J£] c ) s.t. Re A = 70 j G N*. 

Then, there exists Mq > such that, for any t > 0, we have 

\e tJC \^ e( x)<M (l + t) m - 1 e^ t . 

In particular, Theorem B.3 provides a very simple proof of the spectral mapping theorem used in [31] 
and [21]. Indeed, the self-adjoint operator £ involved in these papers is block diagonal: 



£ = 



d 
£ 2 



and both £ 1; £ 2 have at most two nonnegative eigenvalues. More generally, if L\ and £ 2 are closed self- 
adjoint operators on X verifying assumption (A) and if N : X — > X is a linear bounded operator which is 
compact with respect to C\ and £2, then the self-adjoint operator 

£1 N 

also verifies assumption (A). Indeed, £ is bounded from below (since Af is bounded) and its essential spectrum 
is <T ess (£i) U o~ css (£2) C [S, +00) for some positive 6, since Af is compact with respect to £\ and £2)- In 
[40] (section 7.1 there) and [28] (see Proposition 10 there), a spectral mapping theorem is used for such an 
operator. In [40], the specific algebra of the problem was used, and for [28], the proof relies on the arguments 
in [31], but here again, in both cases, we may use Theorem B.3 to show the same result. 

B.2 Passing from linear to nonlinear instability 
B.2.1 Semilinear type models 

We start with a classical result for "semilinear" equations, proved in section B.3. 3. 
Theorem B.4 Let X be a real Hilbert space, and consider an evolution equation of the form 



GG 



where <j> : X — > X is a locally Lipschitz mapping satisfying = 0(|t>|^) as v — > and A is a linear 

operator which generates a semigroup. We assume that Ac ■ D(Ac) C Xc — > /ias an unstable eigenvalue 
in the right half plane {Re > 0} and a /zmie number of eigenvalues in {Re > 0}. We denote 

7o = sup {Re(/i), n G ct([.7£]c) n {Re > 0}} e (0, +oo) 

ana" fix X ^ cr(Ac) with Re(A) = 70 and an associated eigenvector wc G £)(.4c) sucft that ||Re(wc)||x = 1- 
Assume furthermore that there exists < (3 < 70 and M > snc/i f/ia£ 

|e M U w < M e^+«*. 

T/ien, is an unstable solution. More precisely, there exist K > 0, e > and 6"o > swc/i iftai, for any 
< S < S a , the solution v with initial datum v ln = SRe(wc) G D{A) exists at least on [0, ln(2e /<5)/7o] and 
verifies, for < t < ln(2e /^)/7o; 

\\v(t)-8Rc(e tx w c )\\ x <KS 2 e 2t ^ and \\v(t)\\ x > «5e* 7 ° - K5 2 e 2t ^, 

In particular, for < e < e , we see the instability for t = — In f^r") . 1/ F is a Banach space containing 

70 V 5 / 

X and wit/i continuous imbedding X <-} Y, the trivial solution is aZso unstable from X to Y . 

Let us observe that it is always possible to choose the (complex) eigenvector w so that Re(wc) 7^ since 
for any 9 e I, e t6 w is also an eigenvector. The following corollary deals with the orbital instability. We 
recall that under assumption (A), [JC]c has a finite number of eigenvalues in {Re > 0}. 

Corollary B.2 We make assumption (A) and suppose that J L generates a continuous semigroup. Let y be 
a Banach space containing X and with continuous imbedding X <— » y . Assume moreover that [j7X]c has at 
least one eigenvalue in {Re > 0} and choose X G C with 

Re(A) = 70 = max {Re(», fj, G <r([J£\c) n {Re > 0}| e (0, +00) 

and wc G D(Ac) an associated eigenvector such that ||Re(u>c)||;t = 1. We assume moreover that DJl = 
{T(w, g)4>u, 1 w G M, 5 G G} is a C 1 submanifold of X. We finally suppose that the equation (TL) is 
semilinear in the sense that there exists $ : X — > A" locally Lipschitz continuous such that <£>(v) = 0(||i>|||.) 
as«^0 and 

J(£' - w.Q'JWa,. + «) = - w.Q")(&>.)M + *(«)■ 

Then, there exists K > 0, £0 > andSg > 0, depending only onRe(wc) and9Jt 7 wit/i the following properties. 
For any < S < Sq, the solution u to (H) with initial datum u ln = (fy^, + <5Rc(?«c) G D(A) exists at least on 
[0,ln(2eo/<5)/7o] and verifies, for < t < ln(2£ /<5)/70; 

dist y (u(t),aJl) > |^e t7 ° - K5 2 e 2t ^. 

In particular, the "bound state" solution T(w*i)0 Wl is nonlinearly orbitally unstable from X to y and, for 

< e < £q/K, we see the instability for t = — In ( — — ) . 

7o \ 5 J 

In [38], a similar assertion is made for the orbital instability in the remark after Theorem 2 there, but 
with y = X. For applications to PDE's, the space X may be a Sobolev space H s , and y a space like L 2 or 
L°° for instance. The framework of [35] is the single energy space (for instance H 1 ), but an instability result 
established by tracking exponentially growing modes allows to prove instability from the regular space X 
(H 1 ) to the nonregular space y (L 2 or L°°). Here, we may obtain instability in L 2 . 
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Remark B.l In the framework of [35], where a Lyapounov type functional is used, it follows that the 
instability is seen for a time at most equal to K — , where K is some positive constant. This timescale is 

1 2K 

much larger than the natural one — In I — — ] . 

7o \ 5 J 

B.2.2 Some applications 

We may apply our result to the Nonlinear Schrodinger Equation 

id t ^ + Atf + tf/(|*| 2 ) =0, (NLS) 
or the Nonlinear Klein-Gordon equation 

d 2 * = A* + */(|*| 2 ) (NLKG) 
in R d . We shall consider a nonlinearity / at least C 1 , so that we are in the framework of [35]. 

• A bound state solutions for these two equations is a particular solution of the form U(t) = e JLJt (/> w . The 
instability is in general linked to the fact that 

d f , , ,o , „ „ ,„„ „ s d 



cku 



[ \4> u \ 2 dx<0 for (NLS), resp. ^(fjf \<p u \ 2 dx) < for (NLKG). (B.l) 

J R d duj \ J Rd ) 



The existence of at least one unstable eigenvalue has been shown under assumption (B.l) by [34] for radial 
bound states with an arbitrary number of nodes and in [36] for radial ground states. Corollary B.2 may 
be applied with X = H s (R d ), where s <G N, s > d/2 and assuming that the nonlinearity verifies / G C s+ , 
and y = L 2 (R d ) or L co (R d ). The result in [54] shows the instability of linearly unstable bound states for 

p — 1 

(NLS) (in dimension d — 2) with f(g) — g^s~ by showing the exponential growth of an unstable eigenmodc. 
Our result gives a simple proof of this result, but restricted to the sufficiently smooth cases, namely p an 
odd integer or p > 5 + 2s > 5 + d. For non smooth nonlinearities, the situation is more delicate (see [54]). 
An alternative approach is to combine Strichartz estimates with the growth estimate on the semigroup e tJC 
given in Theorem B.3, as in [28]. 

• Corollary B.2 also applies to the Discrete Nonlinear Schrodinger Equation 

VneZ, id t ->S n +e(* rl+1 -2*„ + *„_i) +*„/(|* n | 2 ) = 0, (DNLS) 

as studied in [53] with the saturated nonlinearity f(g)= ^ ^ , /3 > (existence of travelling wave solution) 

and in [23] (defocusing cubic (DNLS), i.e. f(g) = —fig for some (3 > 0). The numerical analysis in the 
paper [23] shows the existence of linearly unstable bound state solutions. The travelling wave solutions 
numerically obtained in [53] are linearly stable, but it may happen that for other nonlinearities /, some are 
linearly unstable. 

B.2.3 Quasilinear PDE's 

For quasilinear problems, we shall not make restrictions on the smoothness of the nonlinearity. The result 
relies on the strategy of E. Grenier [33] and the works [61, 62]. We consider the evolution equation 

rfn i 

- = J(L u + VF(u)) (E) 

for u : R d — >• W, where F € C 00 (IR ly ,M), with the following hypothesis. The operator J is a Fourier multiplier, 
skew-symmetric on L 2 , into and with domain containing H 1 . There exists a > such that the operator 
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Lq is a Fourier multiplier with domain containing H 2a , symmetric and having a self-adjoint realization on 
L 2 {R d ,W). Moreover, for some C > 0, L verifies: 

^\\uf H „ < (L u,u) L 2 < C\u\ 2 H ,. 

The framework proposed in [62] was for L coercive in H 1 , that is a = 1. For the examples below, we shall 
have a = 1/2 or a = 2, which requires very few modifications to the proof of [62]. We still assume that for 
some group G, there exists a unitary representation of G on X, T : G — > %?J^ C (X), leaving the equation (E) 
invariant. 

We consider a stationary solution of the evolution equation (E), that is some Q <E H°°(R d , R") such that 
LqQ + VF(Q) = 0. We are interested in the stability of this solution. We assume that the commutator 
[J, V 2 F(Q)] is bounded in L 2 , which is the case when J is bounded in L 2 or when d = 1 and J = d x . We 
suppose that for the problem 

OH 

— = J(L u + VF(u a + u) - VF(u a ) + G), 

where u a is smooth, bounded as well as its derivatives and G € C(M, H s ) for every s, we have local well- 
posedncss for s large enough: there exists a time T > and a unique solution in C([0, T],H S ). We moreover 
assume that for some continuous non-decreasing function k : R + — > R + with k(0) = 0, the tame estimate 

\(d«J{\7 2 F(w + v)lv}} 7 dZv) L2 \ < k(H^.,.+Mh.)|4 1 

with |a| < s holds true. In order to control high order derivatives, we finally require that for s large enough, 
there exists a self-adjoint operator M s and C s such that 



sU,v) L 2\ < C s \\u\\ H s\\v\\ H ^ (M s m,u) L 2 > \\u\\ Hs - C s ||u|| HS _ min(<ril) 

and 

Re(JLu, M s u) L 2 < C s |u||_y s llMll^s-mint^.i) 
(for a criterion which ensures the existence of such a multiplier, see Lemma 5.1 in [62]). 

Adapting the strategy of [61, 62], we may deduce the following result. Since the proof is very similar, we 
omit it. 

Theorem B.5 We make the above assumptions and moreover that Lq + W 2 F(Q) verifies hypothesis (A) 
in L 2 . We assume furthermore that [J{L + V 2 F(Q))]c has an unstable eigenvalue in the right half-plane 
{Re > 0}, denote 

7 o = sup {Rc(A), A e <7([J(L + V 2 F(Q))] C ) n {Re > 0}} e (0, +oo) 

and fix A e cr([J(L a +\7 2 F(Q))] c ) with Rc(A) = 7 and an associated eigenvector w c e D([J(L +V 2 F(Q))]c) 
such that ||Re(wc)||ff« = 1- There exists sq € N such that, if s > so, Q is nonlinearly unstable from H s 
to L 2 and to L°° : there exists K > 7 Eq > and S > such that, for any < S < <5 , the H s solution 
u to (E) with initial datum u ln = Q + 8Re(wc) € H s exists at least on [0, m(2e /(5)/7o] and verifies, for 
0<t< ln(2 £o /<5)/7o, 

\\u(t) -Q- SRe(e tx w c )\\ H ° < K5 2 e 2t ^, 

hence 

\\u(t) - Q\\ L 2 > Se^« - K5 2 e 2t ~>° and \\u(t) - Q\ L oo > de^" - K5 2 e 2t ^°. 

If, in addition, SOT = {T(g)Q, g <E G} is a C 1 submanifold of H s , then we also have 

dist L 2(u(i),9H) > KSe^ - KS 2 e 2t ^ and dist^oo (u(t),DJl) > KSe^° - KS 2 e 2t ^°. 

In particular, for < e < So/K, we see the nonlinear orbital instability for t = — In ( — — \ . 

7o V S J 
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B.2.4 Some applications to nonlinear dispersive wave equations 

Some model quasilinear equations are given by wave equations (in one space dimension) such as the gener- 
alized Korteweg-de Vries equation 

d t u + d x (f(u)) + d 3 x u = o, ( g Kdv) 

the generalized regularized Korteweg-de Vries equation, also called Benjamin-Bona-Mahony equation or 
Peregrine equation when f(u) = u 2 /2, 

8 t u + 8 x u + 8 x {f{u)) - d t d 2 x u = 0, (gBBM) 

the generalized regularized Boussinesq equation 

d 2 u - d\u - d 2 x (f(u)) - 8 2 t d 2 x u = 0. (grBsq) 

Each of these equations admit nontrivial solitary wave solutions u(t, x) = U c (x — ct) for c in (0, +oo), (1, +oo) 
and (—oo, —1) U (l,+oo) respectively. For these solitary wave solutions, the momentum is, respectively, 

P{U C ) = f U 2 dx = ||£/ c ||| 2 , P(U C ) = [ U 2 + (d x U c f dx, P(U C ) = c [ U 2 + {d x U c f dx. 
Jr Jr Jr 

The existence of exactly one unstable eigenvalue has been shown with the use of an Evans function by R. 

Pego and M. Wcinstcin in [57] for these three equations under the condition ^^ c ^ <- q The paper [49] by 

O. Lopes also gives a linear instability result. Equations (gBBM) and (grBsq) turn out to be semilinear due 
to the regularization effect. Indeed, they may be written 

d t u + (i - d 2 x )- l d x u + (i - dly'dMiu)) = o, d 2 u - (i - dlr^lu - (i - dlr'dluiu)) = o. 

Therefore, Corollary B.2 applies to these two models and this shows the nonlinear instability when linear 
instability holds. 

In [48] , some generalizations of the equations (gKdV) , (gBBM) and (gBBM) have been proposed that 
take into account pseudodifferential operators. These are respectively 

d t u + d x (/(«)) -d x Mu = 0, (I) 

d t u + d x u + d x (f(u)) + d t Mu = o (ii) 

and 

d 2 u~d 2 x u-d 2 x {f{u))+d 2 Mu = 0. (Ill) 

Here, Ai is a Fourier multiplier of symbol M.: A4w = Atw (here, " denotes Fourier transform). We assume 
M > (otherwise, see [48]). When M = -d 2 , these equations reduce to (gKdV), (gBBM) and (gBBM) 
respectively. The Benjamin-Ono equation (M = |^|), the Smith equation (M = \J\ + £ 2 — 1) and the 
Intermediate Long Wave (or Whitham) equation (M. = $/tanh(^_ff) — 1/H, for some constant H > 0) are 
common models of dispersive wave equations that are of type (I). We refer to [48] for references on these 
models and the existence of solitary waves. The associated momentum is 



Pi(U c ) = [ U 2 dx = \\U C \\ 2 L2 P U (U C ) = ||(1 + M) 1/2 C/ C ||i 2 Pm(U c ) = 4(1 + M) 1/2 U, 
Jr 



cilia- 



For these models, Evans' function type arguments do not work since we no longer have a differential equation 
(it is nonlocal). The paper [48] by Z. Lin proposes another approach than the Evans' function technique for 
establishing the existence of unstable eigenvalues. However, it is not completely clear whether this method 
extends easily to the case of systems such as the Euler-Korteweg system (EK) (given at the beginning of 
section 5.1). 
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Theorem B.6 ([48]) We consider one of the equations (I), (II) or (III) with f of class C 1 satisfying 
/(0) = /'(0) = and \f(u)\ \u\ for \u\ — > +00. We assume moreover that M. is even, nonnegative, 

and verifies, for some m > 1, < lim < lim < 00. Assume that c >->• <i c = U c (x — ct) is a 

C 1 branch of travelling wave solution to (I), (II) or (III) wii/i E/ c G i/ m / 2 (R) defined near c* and suppose 
that the linearized operator L has exactly one negative eigenvalue, that ker£ is spanned by d x U Cr and that 

— < 0. Then, U Ct is linearly unstable. 

dc | c =c» 

It is not easy to determine whether the hypothesis of Theorem B.6 hold true when Ai is not a (differential) 
Sturm-Liouville operator. See however [2] on this question. It is clear that if the assumptions of Theorem 
B.6 are verified, then assumption (A) is also satisfied. As for the (gBBM) and the (grBsq) equations, the 
equations (II) and (II) turn out to be semilincar, thus we may prove nonlinear orbital instability by applying 
Corollary B.2. 



The Kawahara equation (or fifth order (KdV) equation) 

d t u + d x (f{u))+adlu + pdlu = Q, (K) 

with a, f3 7^ two real constants, is another relevant dispersive model. For this equation, it may happen 
that the linearized equation around the solitary wave has more than one negative eigenvalues, in which case 
the works [35, 36, 49, 48] do not give a clear necessary and sufficient condition for stability. The paper [14] 
by T. Bridges and G. Derks gives a sufficient condition for linear instability for solitary wave solutions, but 
also for other types of travelling solutions. This condition is probably not necessary since it may happen 
that there exist at least two unstable eigenvalues, or two complex conjugate eigenvalues. 

Instead of stating a general result for nonlinear orbital instability, we shall consider several model cases 
on which we will verify the hypothesis of Theorem B.5, in particular the question of the existence of the 
multiplier M a . 



Proposition B.l We consider the equation (I), namely 

d t u + d x (f{u)) - d x Mu = 

with f of class C 1 satisfying /(0) = f'(0) = and \f(u)\ » |u| for \u\ — > +00. We assume that M is one of 
the following functions: 

-£ 2 (KdV); C 4 + < 2 (Kawahara); |f| (Benjamin - Ono); 



-!- (Intermediate Long Wave); \J\ + £ 2 — 1 (Smith). 



tanh(fff) H 

There exists so > such that, if there exists c € M such that (I) has a nontrivial solitary wave U c € L 2 
which is linearly unstable, then for any s > s , it is also nonlinearly unstable from H s to H s , to L 2 and to L°° . 

By application of Theorem B.5, we are thus able to show the nonlinear instability from H s to L 2 or 
L°° by tracking the exponentially growing mode (this question was left open in [48] and also in [49]). In 
particular, we obtain the L 2 nonlinear instability of the linearly unstable solitary waves for these models. 

Proof. All the assumptions for Theorem B.5 for these types of models are verified in section 8.1 in [62], 
except the existence of the multiplier M s . 

For (KdV), where a = 1, we shall take (for s > 2 integer) 

m s ee (-iyd 2 x s + li^-ir^rM/wr 1 -}, 
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as the computations from [62], section 8.1 show. For the Kawahara equation, with a = 2, we take (for s > 4 
integer) 

m s (-iyd 2 x s + !±^(-i)«-^-2{/'(Q)3j-2. } 

and since the computations are very similar, we omit them. For the Bcnjamin-Ono equation, we have 
= |£| and a = 1/2, and we will then have to deal with pseudo-differential operator which are Fourier 
multipliers with homogeneous symbol. For this type of operator, we shall need some commutator estimates. 
We denote by &{w) or w the Fourier transform of w, and the Fourier multiplier with symbol — zsgn(£) 
(this is the Hilbert transform). 

Lemma B.l (i) Let h G L°° (R) with ,^{M^h) G L\R) (for instance, h e H a (R) for some a > I). Then, 
there exists C > such that, for any v G 



\M*(hv) - hMzv\\ L 2 (R) < C\\v\\ L 2 



2 W ■ 



(ii) Let h G L°° (R) with &{M^h) G L X (R) (/or instance, h G # CT (R) /or some er > 2). Then, there exists 
C > sucft iftai, for any v G i? 5 (M), 



L-'( 



<c|«| 



£ 2 ( 



Cmj Le£ ft G L°°(R) toitft ^(d x M^h) G L 1 (R) (/or instance, h G £P(R) /or some a > 2j. Then, there exists 
C > smc/i £ftai, /or any u G i? 5 (R), 



c^M^/u,} _ hd x Miv - -[d x h]M*v 



<C\\v\ 



Proof. We have 



&(M*(hv)-hM*v){0= f \tfh(£-Qv(Q dC- / |Cl'^-0*(0 dC- 

V y JR JR 



Using the inequality 



-|CI * < C|C — C! = , we thus obtain 



|/iK-OI-l*(0|dC = c{|^(M»/i)|*|*IKO 



and we conclude with the classical convolution estimate L 1 * L 2 C L 2 . This argument does not provide the 
sharpest bound in ft, since it involves ||=^"(.A/f 5ft)| L i, whereas the use of paradifferential calculus will use only 
||ft|| i . However, we shall to use this refinement here. 
The starting point for the second inequality is 

|d§-|C|i-|ci»sgn(C)(€-0| <C\£-(\i. 
To prove this, note that by homogeneity £ = 8(, 9 G R, it suffices to prove 

0|i-l-!(0-l)|<C|0-l|§, 
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which is easy. Therefore, 

^(mHIiv} - hM$v - -[d x h]M* Jfvjit) 

= [ \tfkt-Qv(QdC- [ K\ih^-Ov(C)dC- /||C|»sgn(C)K-C)^-0*(C)dC 
Jr Jr Jr 

<C [ \t-(\^\h(t-()\-\v(()\d( 
Jr 

= C\&(M$h)\ * \v\, 

and we conclude as before. For the third inequality, we argue in a similar way with the estimate 

^ici* — *cici* — *§ici*c^ — o| < ci^ — ci s - 

The proof is complete. 



□ 



For the Benjamin-Ono equation, A4(£ = a = 1/2 and the index s will be half an integer: s <G N/2. 
Therefore, we set s = [s] + {s}, with [s] integer and {s} € {0; 1/2}. Let us define, for s e N/2, s > 1, 

(-l)*e +lsMh s x - 1 {f'(Q)d s x - 1 M^} if {s} = 

(-l)M^MM +7s aW{/'(Q)9W-} if { S } = I, 

for some real constant 7s to be determined later. It is clear that M s is self-adjoint on L 2 and that there 
exists C s > such that 



\(M„u,v) L *\ < C,\u\h-Mh- 



and 



(M s u,u) L 2 > \\uf H s - C s \\u\\ 



H 2 

2) 



In order to verify the assumptions for the multiplier M a , it remains to study Re(J(Lo + V f(Q))ii,M s !i)y. 
When {s} = 0, i.e. s € N, this quantity is 

Re(d x (M + c+f'(Q))u,M s u) L 2 

= Re(d x Mu,(-l) s d 2 x s u) L 2 + ls Rc(d x Mu,M l zd x ~ 1 {f(Q)d x - 1 Miu}) L 2 
+ Re(^[/(Q) U ],(-l) s ^ s u) i2+7s Re(^[/'(Q) U ],A^5^- 1 {/'(Q)^- 1 Al^}) i2 
+ cRc(d x u,M s u) L 2. (B.l) 

By skew-adjointness, the first and last scalar products are zero. By integration by parts and Leibniz formula, 
we deduce, since Q <G H°°, 

Re(d x [f(Q)u], (-l) s d 2 x s u) L 2 = Re(d s x +1 [f'(Q)u],d s x u) L 2 

<Re(/'(0)^ +1 U ,a>) L 2 + (,s + l)Rc(^[/'(Q)]a>,a>) L2 +C s || U ||^||u|| H3 -i 

< [s+ ^)Re(d x [f(QW x u,d s x u) L 2 + <7 a |uM«|ff.-i- 
Similarly, using the easy estimates ||A^ < and < C(/i)||u|| H i for h E L°° with 



"H2 



"H2 



3?(M2h) e L 1 (this is an immediate consequence of Lemma B.l) 

ls Re{d x [.f l {Q)u],M h d x - 1 {f'{Q)d s - 1 M^u}) L 2 = ls {-iy- 1 Re{M^d x [f{Q)u],f'{Q)d s - 1 M^u) L 2 
< ls (-ir- 1 Re(Mi[f , (Q)d x u}J'(Q)d x - 1 Miu) L 2+C\\u\\ 



II 2 
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Using Lemma B.l, we deduce \M^[f [Q)d s x u] - f{Q)M^d s x u\\ L 2 < C(Q)\u\ H >, thus 

7.Re(0 x [f'(Q)u] , M * d s x - 1 {f'(Q)d s x - 1 M ^u}) L , 

< ls {-l) s - 1 Re{f\Q)d x M^uJ\Q)d x - 1 M^u) L 2+C\u\ H s\\u\\ uS _ h 



= |(-i)^Re(a x [/'(Q)]ar 1 x^ )/ '(Q) 9r ixl u ) i2 + c|| w ||^lkll HS -i 

We now turn to the term 

7 a Ite(0 I A4u,A4»^-V(Q)0r 1 *^ 
Using Lemma B.l, we write 

<C7(Q)|«| £r ._i, 

which implies 

j s Re{d x Mu,M^d s x - 1 {f{Q)d s x - 1 Miu}) L 2 

< Js (-iy- 1 Re(d x u,f'(Q)d x - 1 M 2 u) L 2 + ^,(-l)- 1 Re(^ U) fl I [/'(g)] J M^{y i - 1 Miu}) 1J 

Noticing that M 2 = -<9 2 and M^J^d^M^ = d^MJ? = -d s x (since has symbol equal to 

we infer 

ls Re(d x Mu,M^d x - 1 {f(Q)d x - 1 Miu}) L 2 

< 7s (-l) s Re(a>,/'(Q)^ +1 U ) i2 + * ls (-iyRe(d x u,d x {f(Q)]d s x u) L 2 



by integration by parts. 

Reporting these estimates into (B.l), we infer 

Re(0 I (M + c + / , (Q))u,M i u) La < (s + i)Re(cU/'(Q)]c>>, 9» L2 + 7a (-l)'Re(^u, d x [f {Q)]d s x u) L 2 

+ C\\ul H s\\u\\ H3 _i. 

Therefore, the choice 

7. -(-I)- 1 (.+ 5) 

provides the desired control 

Re(d x (M + c + f'(Q))u, M s u) L 2 < C\u\\ H s \u\ H ._ 1 . 
When {s} = 1/2, the computations arc similar: (B.l) becomes now 

Re(d x (M + c+f'(Q))u,M s u) L 2 

= Re(d x Mu,(-lpd^Mu) L 2 + ls Rc(d x Mu,d^{f(Q)d^u}) L 2 

+ Rc(a x [/'(Q) u ],(-i)Wa 2 WA^ u ) L2+7s Re(a x [/'(Q) u ],a[ s H/(g)aw u }) i2 

+ cRe(d x u,M s u) L 2, (B.2) 
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and the first and last scalar product still vanish. Moreover, by integration by parts and Leibniz formula, we 
deduce, since Q e H°°, 



7, 



,Re(a x [/ / (Q)n],aW{/ / (Q)aW4) i2 = 7s (-l) [sl Re(aW +1 [/ / (Q)n],/ / (g)4 s ] U ) i2 
< 7s (-l) [sl Rc(/'(Q)aW +1 U , f{Q)d^u) L , + C\\uf H[s] 



<7a(-i) w - 1 Re(a x [/ , (Q)]aw«,/ , (Q)aw«) ia + q«|| rW 

<C\\ut 2 H[s] =C\\u\\ 



|| 2 _ /-f||„.||2 



Furthermore, 

Re(5x[/'(Q)«],(-l) w ^ w M«) La =Re(M3aW+ 1 [/'(Q)u],aWM5«) i2 

< Re(9 x 7W 3 {/' (Q)aW u } ; flW^i U ) L2 

+ WRe(a x Af^a x [/ , (Q)]5W-M,^ 1 ^»«)L»+GH ffW -il«l ff [.] + j- 

For the second scalar product, we write, by Lemma B.l, 

Re(^M s R [f'md^-'u} , 0M 

= Re(A^H^[/'(g)]4 sl " 1 «},4 sI ^^) L 2+Re(A^^a x [/ / (Q)]4 sl «},4 sl >'^) L 2 
<C|| U || ffH _i|h|| ffW+ i +Rea[/'(g)]4 s ]>(^,aWAll M ) i2+ C|| U || H[sl || W || ff[3l+ i 

<Re(a I [/ , (Q)]aWAf5u,aWA45«) i2 + q«| ff ._i|u| ff .. 

For the first scalar product, we use Lemma B.l once again: 

Re(d x M 1 Hf'(Q)d^uh dWM k >u) L 2 

<Re(f(Q)d x Mid^u,d^Miu) L 2 + ^Re(d x [f'(Q)]d^Miu, d^M^u) L , + C\u\ H ._i \u\ H - 
<Re{d x [f{Q)]d^M^u,d^M^u) L ,+C\u\ H3 _ h \u\ H s. 

As a consequence, since [s] = s — \, 

Re(d x [f'(Q)u],(-l)^d 2 x ^Mu) L 2 < ( s + ^jRc(d x [f'(Q)]d^M^ Ul d^M^u) L2 +C\u\ H ._i \u\ H .. 
We turn finally to the term 

ls Re(d x Mu,d^{f(Q)d^u}) L , = ls (-l) [s] Re(d^M^u,d x MHf'(Q)diM)L^ 
and infer, by Lemma B.l, 



7. 



: Re(d x Mu,dW{f(Q)dWu}) L 2 <j s (-l)^Re(d^M^u,f(Q)d x M^d^u) L 2 



+ ^7 s (-l) W Re(4 sl ^^ I ^[/'(Q)]Al^W u ) L2+ C'||u|| H3 _ i ||u||^ 
= 7 5 (-l) [sl Re(9 x s ^^,9x[/ / (Q)]^^W U ) i2 +C||n|| ffS _i||n||^. 



Therefore, 

Re(d x (M + c + f(Q))u,M s u) L 2 < (« + * + 7a (-l)M)Re(0WM*u, ^[/'(Q)]^^)^ + CH ff .-l H| ffs , 
hence choosing 7s = (— l)M _1 (s + |) gives the result. 
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It remains to study the cases of the Intermediate Long Wave equation and the Smith equation, for which 
M is, respectively, 

tanh(^) H ' V + ^ 

We denote Mo the operator with symbol |£| (the one of the Benjamin-Ono equation), and define M s as for 
the Benjamin-Ono case (hence with "M" = Mo)- We observe that in both cases, M = M — Mo is bounded 
on L 2 . Indeed, its symbol is continuous in K and, for £ — > ±00, 

• M «» - ssfeo - h - sg „(o + o(e-^«) - 5 - « - 1 + w ' e " 2|il ") 

and 

^t(o = a/t+^- 1 = leiv/r+r^ - 1 = Ki - 1 + Oder 1 ) 

respectively. In the quantity Re(d x (M + c+ f (Q))u,M. s u) , we then have to bound from above the extra 
term Re(d x (Mu) ,M s u) l 2 , that is (using the skew-adjointness for the higher order derivatives in M s ), 

Re(d x (Mu), ls Mld s x - 1 {f'(Q)d s x - 1 Mlu}) L , = j s (-l) s ~ 1 Re(d x M^ (Mu), f'(Q)d x ~ 1 M^ u) L 2 if { s } = 0; 

Re{d x {Mu), ls d^{f{Q)d^u}) L , = ls {-l)WRe{dW + \Mu),f (Q)d^u) L , if {s} = \. 

We then note that in both cases, one may actually split M = M — Mo = M c + Mh, where M c is the 
multiplication by —1/H (resp. —1) and Mh has a symbol which is continuous in R and 0(|£| _1 ) at infinity, 
so that Mh is bounded from H a to H a+1 if a > 0. Therefore, when {s} = 0, we easily get 

Re(d x (Mu), 7s Mld s x - 1 {f'(Q)d s x - 1 Mlu}) L 2 

= ls (-iy- 1 R C (dlMl{M c u)J'{Q)dl- 1 Mlu) L >+ ls (^^ 

< ^.{-l)'MS' x - 1 Miu,M c d x \f'(Q)]d' x - 1 Mlu) L , +C\u\ 2 H ,_ h < Cluf HS _ h 

and similarly when {s} = 1/2. Therefore, the estimate 

Re{8 x {M + c + f{Q))u,M s u) L 2 < C\u\ H ._ 1 \\u\\ H s 

remains true for the Intermediate Long Wave equation and the Smith equation. The proof of Proposition 
B.l is thus complete by applying Theorem B.6. □ 



B.3 Proofs 

B.3.1 Proof of Theorem B.3 

We shall prove the resolvent estimate required in Corollary B.l. Let us consider A = 7 + it e C with 7 / 
and the resolvent equation {JC — \)v = S, or 

(~f + ir)v = JC(v)-H. (B.3) 

By hypothesis, the essential spectrum of JC is of the form i[R \ (—-do, +#o)]- Moreover, we have seen that 
JC has a finite number of eigenvalues in the half-space {Re > 0}, hence, for |t| > t sufficiently large, we 
know that there exists a unique solution v to (B.3). By taking the scalar product with C{v), we deduce the 
conservation law 

1 {v,C{v)) x = -Re{V,C{v)) x . (B.4) 
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By our assumption, there exists a finite (possibly empty) number of eigenvalues in (— oo,0], {—Hi, fJ-q), 
each one of finite multiplicity. For any 1 < k < q, we fix an orthonormal basis {xk,t)i<i<n k of the eigcnspace 
Ker(£ + fik). By assumption (A), any eigenvector Xk,e is smooth in the sense that \k,l S D{J) and 
J X k,e e £>(£). 

We then make a spectral orthogonal decomposition 

q «fc 



fc=i £=i 

where £(xk,t) = HkXk,t and (v + ,£(v + )) x > ( 5|| w +|| 2 t f° r some positive S. In the double sum, we have a finite 
number (independent of v) of terms. Inserting this into (B.4) yields 

W\v+\% < h\S(v + ,£(v + )) x < s[\Re(^£(v)) x \+J2^\a k jf] < K[X\ x \v\ x + Jr£> M | 2 

k,£ k,l 

Using the inequality ab < ea 2 + b 2 /(4e) with a = \\v\\ x , b = K\\T,\\ X and e = |7|<5/2, the equality \\v\\ x = 
\\ v +\\x + J2k i l a Ml 2 and incorporating the term |7|<5|u+|^/2 in the left-hand side, we infer 

^Xv+^KK'Y.WA 2 +K"\n x - (B.5) 

k,e 

On the other hand, since Xk.e € D{J) and Jxk.e € D{£) by assumption (A), taking the scalar product of 
(B.l) with Xk,e provides 

(7 + ir)a k j = -(v, £Jxk,t)x - (£, Xk,t)x- 

Consequently, 

(H + M)|a M | <if M H^+ A-|EU, 

thus 

(l7l + |r|) 2 £ K,,| 2 < K \\v\\% + K\m% =K J2 + K\v + \% + K\E\\%, 

k,e k,e 

which implies, if |r| > 1 + \fK§ — \y\, 

2 
x 



Vu 12 < y IKI* + ll^ll 



Reporting this into (B.5) gives 



If |t| > 1 + v^r+4^VhP - |7|, we deduce 



and it follows that 

\\v\ x = \\v + \ x +Y J M 2 <K2\\n x , 

where K 2 does not depend on |r| (large enough) as wished. 



The proof of the first semigroup estimate follows then easily, see, for instance, Proposition 2 in [59] . 
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Proof of the semigroup estimate when 70 > 0. Here, we assume 70 > 0. As a consequence, the 
spectrum of [JC]c is of the form <r s U a u , where cr 0SS ([JX]c) C <r s C {Re < 0} and 7^ er u C {Re > 0} 
consists in a finite number of eigenvalues of finite algebraic multiplicities. Therefore, we may define (see, e.g. 
[41], [39]) the spectral Riesz projection 



1 

2^ 



([JC] C - z)- 1 dz, 



where Y is any simple (positively oriented) closed curve enclosing ct u . As a consequence, P is bounded, 
commutes with [J £} c on D([JC] C ) and verifies a([J£] c P) = cr u , a([JC] c (ld - P)) = a s . Moreover, [JC] C P 
is bounded, hence generates a continuous semigroup, e'^^ 1 ", given by the exponential series 



t n ([J£] c P) n 

n=0 



In addition, [J r £] c (Id — P) = [JC\c — [JX]cP also generates a continuous semigroup and that we have 

e t[JC]c _ e tJCF e t[JC]c(ld-P) ^ 

The semigroup generated by the bounded operator [J7X]cP is easily analyzed. We shall now apply the 
spectral mapping theorem of J. Priiss (Theorem B.2) to [JX]c(Id — P) in order to control the growth of its 
norm. By Corollary B.l, it suffices to estimate its resolvent [[JX]c(Id — P) — (7 + *7")] _1 for large |r| (note 
that <r([JX]c(Id — P)) = er s C {Re < 0}. If £ £ Ac and |t| is large, it is clear that the solution u e Ac to 
[[J£] c (Id - P) - (7 + ir)]u = S is given by 



u = [[JC} C - (7 + ^T)]- 1 (Id - P)E - 



7 + IT 



-PS, 



thus, for |r| large, 

[[^£] c (Id-P)-(7 + »r)]- 1 



< 



[[J£]c-(7 + *r)]- 1 



is bounded. Consequently, by Theorem B.2 and since cr s C {Re < 0}, ^e'l^W' 11 -')) = e t<r([JC] c (id-P)) = 
e tu s c £)(o ? l). It follows that for any e > 0, there exists K e > such that 



e t[J£lc(Id - p) |^(^)<^e et . 



V t > 0, 

Since e*[- 7£ ]« :P is given by the exponential series, we also have the optimal estimate 



V t > 0, 



s t[XC]cP|| 



by definition of m. We conclude by taking e = 70/2 for instance. 



B.3.3 Proof of Theorem B.4 

Since A generates a continuous semigroup, v is a solution to d t v = Av + $(u) if and only if it is a mild 
solution: 

v(t) = e tA v m + f e (t - T ^ A ^(v(T)) dr. 
Jo 

There exists r > such that ||$(w)||x < M||u||^ if \\ v \\x < r o- We choose v ln — <5Rew, where ||Rew||x = 1 
and w is an eigenvector for the eigenvalue A, and write the solution under the form v = e tA v ln + v = 
Re(e tA w) + v. If A E R, we can choose w € D(A) C D(Ac)- Then, 



v(t) 



f e < t - T )- A $( ( 5Rc(e* A «;) + v(t)) dr. 
Jo 
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Let us denote ri = min (r , 



70-/3 



_ i , ^ and let T > be the maximal time such that T < ln(r/(2<5))/7 and 

ll^( T )llx < r i/2 m [0, T"), where < r < r\ will be determined later. We shall work for < t < T, so that 
|<5Re(e tA w) + v(t)\\ x < 5e^° + n/2 < n < r . Then, 





e (*-r)^ 


/o 





M|<5Re(e rA w;) + w(T)|^ dr 



< 



< 



2M M ^*e(70+/3)(t-r)^2 e 2r 70 + | 5 ( r )|2^ dr 



2M M 



7o 



f oM 5 2 e 2t 7o+riMoM /" 

— P Jo 



a (70+/3)(t-T)||~ 



i5(r)[x dr, 



since /3 < 70. Applying now the Gronwall inequality to e ^ 7o+/3 ^*||v(t)||x then gives, since M Mr\ < 70 — /3, 



2M M 



2nM 2 M 2 



+ 



7o~/3 (7o-^)(7o-^-riM M) 



We now choose r = yJr\/K, so that the right-hand side is < Kr 2 /A < ri/2, and this implies that u exists 
at least on [0,ln(r/(2<5))/7 ]. In addition, for < t < T, 

\\u(t)\\ x > Se 1 '' - \u(t)\x > Se^« - K8 2 e 2t ^, 

as desired. We conclude choosing e > so small that 2e — Kef, > e . 



B.3.3 Proof of Corollary B.2 

We pick some < f3 < 70 (for instance f3 — 7o/2) in order to have the semigroup estimate required in 
Theorem B.4. The solution u(t) = T(u}*t)(<j> u , + v{t)) verihes, for < t < 7^ ln(2e /£), 

M*)l* = \\T(-ojJ)u(t) - + 5Re(e tx w))\\ x < KS 2 e 2t ^. 

Hence, T(— u*t)u(t) remains at distance < Keq from <j> Ut € SOT and therefore 

dist*(u(t),07t) > dist x (SRc(e tx w),m- 4>) - K8 2 e 2tla . 

Assume Ael. Then, we observe that the straight line R 3 9 Qw is transverse to the tangent space T^SPT 
of the manifold Wl, since w is an eigenvector of JC for A ^ 0, hence does not belong to the kernel of £. 
Therefore, &st x {6w, M - 4>) > \9\fKx for small \6\. Thus, 

dist*(u(t),0R) > 7^<5e tA - K5 2 e 2t ~<«. 

If A € C \ R, the equation [JX]c(w) = Aw splits as JC(Rew) = Re(A)Reu> — Im(A)Imu) and JC(\nvw) = 
Im(A)Rew + Re(A)Imro. Therefore, Rew and Imw do not belong to ker(£). Consequently, the surface 
C 3 9 Rc(9w) is transverse to the tangent space T^DJl of the manifold 9Jt, and we conclude as before that 

dist*(u(t),S!R) > -^6e^° - KS 2 e 2t ^°. 
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